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PAAS  USER  MANUALS 


In  this  appendix  the  user  manuals  for  the  PAAS 
software  are  presented The  user  manuals  are  arranged  in 
alphabetical  order  and  each  is  self  contained  except  for 
Table  1.  Table  1 is  a list  of  Honeywell  CCOS  file 
management  supervisor  status  codes  which  are  used  to 
indicate  errors  in  file  handling  to  the  user.  Since  disc 
files  (PRMFL's)  are  used  by  all  PAAS  software  modules  this 
table  is  printed  one  time  here  for  reference  by  all  the 
user  manuals.  All  of  the  modules  discussed  herein  except 
PLTDVR  can  be  used  via  any  remote  terminal  compatible  with 
the  RADC  H6180  CCOS/TSS . The  module  PLTDVR  is  designed  to 
transfer  plots  to  the  Dedicated  User  Interface  Subsystem 
(DUIS) . 


All  TSS  file  references  made  herein  are  to  the  user 
master  catalog  BECAGD01  in  the  RADC  H6180  CCOS  system.  In 
this  documentation  all  RUN  and  OLD  commands  are  shown 
assuming  the  user  has  previously  accessed  the  named  file  or 
has  a copy  of  the  file  under  his  user  master  catalog.  The 
procedure  for  accessing  files  that  are  stored  under 
BECACD01  is  as  follows: 

From  SYSTEM  level  enter  YFORT  and  then  enter  the 
ACCESS  subsystem  by  typing: 

ACCE 


The  HS180  will  reply 


Function? 


Enter  the  following: 


AF , BECAGD01/PLARY , R 


. •». , aclLsizSi 


TABLE  1 


HONEYWELL  GCOS  FILE  MANAGEMENT  SUPERVISOR  STATUS  CODES 


Status  codes: 


4000  NO  ERRORS 

4001  NAME  NOT  III  MASTER  CATALOG 

4002  I/O  ERROR  Oil  DEVICE  XXX  SA  = NNN NNN 

4003  PERMISSIONS  DENIED 

4004  FILE  BUSY:  TPY  LATER 

4005  INCORRECT  CAT/FILE  DESCRIPTION  AT  AAA AAA 

4006  LLINK  SPACE  EXHAUSTED,  DEVICE  XXX 

4007  UNDEFINED  DEVICE  YYY  ZZZZZZ 

4010  LINK  SPACE  EXHAUSTED,  DEVICE  XXX 

4011  NON-UNIQUE  NAME 

4012  SIZE  REQUESTED  LS  THAN  ALLOCATED 

4013  SPACE  REQUEST  GR  THAN  ALLOWED 

4014  PASSWORD  REQUIRED  AT  AAA AAA 

PASSWORD  AAA AAA  AT  AAA AAA  INCORRECT 

4015  I-D-S  FILE  IN  ABORT  STATUS 

4016  FILE  CANNOT  BE  ALLOCATED  FOR  TS  USE 

<017  SEEK  ERROR  ON  DEVICE  XXX  SA  = NXN NNN 

4020  FAILURE  IN  NAME  S CaN  (IMP • ) 

4021  UNDEFINED  DEVICE  (IMP.) 

<022  DEVICE  LINK  TABLE  CHKSUM  ERROR 
<023  INCONSISTENT  FSW  BLOCK  COUNT 

4024  INTERNAL  LINK  TABLE  CXSM  ERROR 

4025  REQUESTED  ENTRY  NOT  ON-LINE 

4026  NON-STRUCTURED  FILE  ENTRY 

4027  FILE  IN  EFFECTIVE  STATUS 

4030  ILLEGAL  PACK  TYPE 

4031  ACCESS  GRANTED  TO  1-D-S  FILE 

4033  CAT/FILE  SECURITY  LOCKED 

4034  ILLEGAL  CHARACTER  IN  CAT/FILE  NAME 

4035  ILLEGAL  CAT/FILE  LIST  REQUEST 

4036  AFT  IS  FULL 

4037  FILE  ALREADY  IN  AFT 

4040  MAXIMUM  PAT  SIZE  EXCEEDED 

4042  INVALID  FILE  CODE  OR  PAT  PTR 

4043  INVALID  CATALOG  3LOCX  ADDRESS 

4044  PERMISSION  DENIED  - SHARED  FILE 

4045  INVALID  SPACE  IDENTIFIER 

4051  CHECKSUM  ERROR  - DEVICE  XXX  SA  * NNN NNN 

4052  DEVICE  XXX  RELEASED 


WHERE:  XXX  ^DEVICE  NAME  (STl,DSl, . . . ) 

NNN NNN  -OCTAL  REPRESENTATION  OF  THE  SEEK  ADDRESS 

AAA AAA  =12  BCD  CHARACTERS  OF  THE  CATALOG  ELEMENT 

IN  ERROR 

YYYY  =TYPE  /OR/  NAME 

ZZZZZZ  -DEVICE  NAME  OR  CLASS  OF  DEVICE 


The  H6180  will  reply: 

Successful 

Function? 


Type  a carriage  return  and  control  will  return  to 
YFORT.  In  the  above  example  the  file  named  PLARY  (shown 
underlined)  was  accessed.  For  any  other  file  Just 
substitute  the  appropriate  file  name. 
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PROGRAM  FFT2DX 


I .  PURPOSE 


The  program  FFT2DX  transforms  an  antenna  aperture 
distribution  into  its  equivalent  far-field 
distribution.  FFT2DX  generates  the  far-field 
distribution  via  a two-dimensional  Fast  Fourier 
transform  of  the  aperture  illumination  function, 
far-field  is  stored  on  a PRMFL  designated  by  the 
user . 

II.  PERIPHERAL  DEVICES  REQUIRED 


The 


UJD 

NAME 

USE 

01 

Input  PRMFL 

Aperture  Input 

02 

Output  PRMFL 

Far-field  Output 

03 

Temporary  disk  file 

Intermediate  FFT 
Storage 

05 

Batch  data  deck 

Program  Input  Data 

06 

Batch  system  output 

Computer  Output 

OPERATING 

PROCEDURE 

1.  Enter  the  TSS  CARDIN  system 

2 . The  message 
OLD  OR  NEW? 

is  printed  on  the  operator's  console: 

3.  Type: 

OLD  FFT2DX 
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The  program  FFT2DX  will  be  loaded  into  the 
current  file. 

4.  Co  to  line  1460  and  change  the  PRMFL  name  that 
follows  the  characters 

* S : PRMFL : 01 , R/V , R , 

to  the  appropriate  PRMFL  name  which  contains  the 
input  aperture  which  was  created  using  PLARY, 
RNDERR  or  FILMOD. 

5.  Co  to  line  1470  and  change  the  PRMFL  name  that 
follows  the  characters 

*1 : PRMFL : 02 , R/V, R, * 

to  the  appropriate  PRMFL  name  in  which  the 
far-field  energy  distribution  is  to  be  stored. 

6.  Co  to  line  1500  and  change  the  value  of  the  input 
parameters  in  the  namelist  FFT  to  correspond  to 
the  specific  Input  aperture  and  desired  far-field 
output . 

N2  - The  power  of  2 which  determines  the  number 
of  points  on  the  side  for  the  2-d  transform 

LRJ  - The  number  of  blocks  in  the  x-direction  of 
the  input  aperture 

LRK  - The  number  of  blocks  in  the  y-directlon  of 
the  input  aperture . 

7.  Co  to  line  1510  and  change  the  values  of  the 
parameter  in  the  namelist  FILOUT  to  correspond  to 
the  desired  blocks  of  the  far-field  output  to  be 
stored  in  the  output  PRMFL. 


nw  »*■  i^/tt  |*-s.-»"«->it^’,,1‘,‘'1" M'lliT&i'* •'.' . l'". ■ — - r > ..  v,-y-  -~«i  »\ 
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LRJIN  - The  number  of  blocks  to  be  skipped 

in  the  x-direction,  starting  on  the  left, 
before  beginning  to  store  the  output. 


LRJVID  - The  width  in  blocks  of  the  desired 
far-fleld  output. 


LRKIN  - The  numne.-  of  blocks  to  be  skipped 
it.  !;h-  y-direction,  starting  at  the 
top,  before  beginning  to  store  the 
output . 


lRKWID  - The  height  .in  blocks  of  the  desired 
fur-field  output. 


?J3 


a.  hr  this  pom’,  tne  job  definition  file  is  complete 
and  the  user  ;uay  SAVE,  RESAVE  and/or  RUN  the 
current  file. 


9.  Type: 

RUN 

to  run  the  prepared  file. 

10.  The  message 
SNUMB  XXXXT 

is  printed,  where  XXXXT  is  the  job  identification 
number  and  is  used  to  learn  the  status  of  the  Job 
at  later  points  in  time. 


r..  -I  .W, . - . 
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IV.  SUBPROGRAMS  REQUIRED 

FFT2D 

V.  RESTRICTIONS.  REQUIREMENTS , AND  MISC . DATA 

1.  This  program  must  be  executed  under  Honeywell 
GCOS  TSS  CARDIN  subsystems. 

2.  N2  must  be  In  the  range  4 £ N2  £ 10. 

3.  LRJ,  LRK  must  be  exactly  the  same  value  as  was 
specified  in  loading  the  aperture. 

4.  In  picking  values  for  LRJIN,  LRJVID,  LRKIN, 

LRKWID  remember  that  the  size  of  the  output  PRMFL 
must  be  large  enough  to  hold  ail  of  the  specified 
output . 

5.  The  input  PRMFL  and  the  output  PRMFL  must  be  in 
different  files. 
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PROCRAM  FILMOD 


I .  PURPOSE 


The  program  FILMOD  allows  the  user  to  modify  a 
previously  generated  aperture  that  has  been  stored  in 
a permanent  disk  file  (PRMFL).  This  program  allows 
the  user  to:  <1)  list  element  values,  (2)  change 
element  values  (one  by  one),  or  (3)  'punch  holes' 
with  a radius  and  center  both  specified  by  the  user. 
The  modified  aperture  is  then  stored  on  a user 
designed  PRMFL  which  may  be  the  same  as  the  input 
file. 


II.  PERIPHERAL  DEVICES  REQUIRED 


LUD 

NAME 

USE 

01 

Input  PRMFL  (Optional 
output  PRMFL) 

Aperture  .Input 
(Optional  outpu 

02 

Output  PRMFL 

Aperture  Output 

09 

Time  sharing  terminal 
keyboard 

User  Input 

06 

Time  sharing  terminal 
printer 

Computer  Output 

III .  OPERATING  PROCEDURE 


1.  Enter  TSS  YFORT 


2 . Type 

RUN  FILMOD 


3 . The  message 

INPUT  FILE  NAME 


W a.  .flMiAtaA  — . *ai»  «r 
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is  printed. 

The  user  responds  with  the  name  of  an  existing 
PRMFL  which  contains  the  aperture  to  be  modified. 
The  file  name  must'  be  followed  by  a semicolon( ; ) . 


4.  If  the  PRMFL  name  is  not  acceptable  to  the 
computer  the  message 

UNSUCCESSFUL  ATTACH  ISTAT  - X 

is  printed,  where  X is  the  first  status  word 
returned  by  the  File  Management  Supervisor  (see 
the  Time  Sharing  System  Programmers  Reference 
Manual , BR39,  p.(3-39)  or  Table  1 herein)  or  will 
contain : 


1 - file  is  currently  open 

2 « teletypewriter  requested  in  batch  mode 

(illegal) 

3 - additional  memory  needed,  request  denied 

(time  sharing  user  will  be  terminated) 

4 - CATFIL  all  blanks 


NOTE:  See  Honeywell  series  600/6000  Fortran 

manual,  BJ67,  p.(6-35)  - (6-36)  for  more 
details  on  the  subroutine  ATTACH. 


After  the  message  is  printed  the  program  returns 
to  step  3 . 

5.  If  the  PRMFL  name  in  step  3 is  accepted  the 
message 

OUTPUT  FILE  NAME 
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Is  printed. 

The  user  responds  with  the  PRMFL  name  in  which  he 
wishes  the  output  to  be  stored.  The  user  may 
specify  the  input  file  and  output  file  as  the 
same  file  if  he  wishes.  The  file  name  must  be 
followed  by  a semicolon( ; ) . If  the  file  name  is 
not  acceptable  the  program  will  go  through  step  4 
and  return  to  step  5 so  that  the  user  may  try 
another  file  name.  Otherwise  the  program  will 
proceed  to  step  6. 


!?<«!,  LRK 
is  printed. 

LRJ  - The  width  of  the  input  aperture  measured  in 
blocks  . 

LRK  - The  height  of  the  input  aperture  measured 
in  blocks. 

NOTE:  The  numbers  entered  for  LRJ,  LRK  must  be 
exactly  the  same  as  the  values  specified 
in  the  program  PLARY  which  originally 
loaded  the  aperture  file. 

The  message 

MODIFY  OR  HOLE?  <0  or  1) 
is  printed. 

• 0 The  program  begins  the  question  end  answer 
sequence  for  modifying  individual  elements 
(Proceed  to  step  8) . 
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■ 1 The  program  begins  the  question  and  answer 
sequence  for  making  holes  in  the  aperture 
(Proceed  to  step  18)  . 

8 . The  message 


IBLK 


is  printed. 

The  user  must  respond  with  the  lowest  block 
number  in  which  he  wishes  to  make  any 
modifications.  Each  subsequent  request  for  IBLK 
must  be  answered  with  a larger  number  than  the 
previous  response. 

9 . The  message 

ANY  ELEMENTS  LISTED?  <Y  or  N) 
is  printed. 

- Y The  user  wishes  to  have  the  values  of  some 

of  the  elements  in  the  specified  block 
1 1 sted . 

- N No  elements  are  listed  (Proceed  to  step  12). 

1C.  The  message 

JSTRT,  KSTRT,  JSTP,  KSTP 
is  printed. 

JSTRT  - The  horizontal  coordinate  to  begin  the 
value  listing. 

KSTRT  - The  vertical  coordinate  to  begin  the 
value  listing. 


12 
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Example:  (JSTRT,  KSTRT)  - (1,  1)  if  the  user 

wishes  to  begin  in  the  upper  left  hand 
corner . 

JSTP  - The  horizontal  coordinate  to  end  the 
value  listing. 

KSTP  - The  vertical  coordinate  to  end  the  value 
listing. 

Example:  (JSTP,  KSTP)  - (16,  16)  if  the  user 

wishes  to  end  in  the  lower  right  hand 
corner . 

NOTE:  The  user  must  give  all  coordinates 

starting  from  the  upper  left  hand  corner 
with  (1,  1)  going  to  the  lower  right  hand 
corner  with  (16,  16).  As  an  example,  if 
the  user  wished  to  list  the  values,  of  the 
whole  block  (256  values)  he  would  type 

1,  1,  16,  16 

in  response  to  the  message. 

The  value  listing  is  printed  in  two  columns  with 
the  real  part  of  the  element  value  in  the  left 
column  and  the  imaginary  part  in  the  right 
column.  The  values  are  double  spaced  with  one 
element  per  line  starting  with  the  element 
(JSTRT,  KSTRT)  and  proceeding  to  (JSTP,  KSTRT). 
The  next  value  will  be  (JSTRT,  KSTRT  + 1)  and  so 
forth  until  (JSTP,  KSTP)  is  reached. 

11 . The  message 

ANY  ELEMENTS  CHANCED?  (Y  or  N) 

is  printed. 

» Y The  user  'wishes  to  change  some  element 

13 
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values  . 

(Proceed  to  step  12) 

- N The  user  does  not  wish  to  change  any  element 
value . 

(Proceed  to  step  14) 

12 . The  message 

HOW  MANY  ELEMENTS  CHANCED? 

Is  printed. 

The  user  responds  with  the  number  (up  to  100) 
that  he  wishes  to  change  in  the  specified  block. 
If  he  wishes  to  change  more  than  100  the  answer 
to  step  14  must  be  Y.  If  the  user  responds  with 
0 the  program  will  jump  to  step  15. 

13.  The  message 

IELJ,  IELK,  VREAL,  VIMC 
is  printed. 

IELJ  - The  horizontal  coordinate  of  the  element 
to  be  changed. 

IELK  - The  vertical  coordinate  of  the  element  to 
be  changed. 

NOTE:  The  coordinate  locations  are  specified 

according  to  the  explanation  given  in  step 

10. 


VREAL  - The  real  part  of  the  new  element  value. 

VIMC  - The  imaginary  part  of  the  new  element 
value  . 
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The  above  message  is  repeated  the  number  of  times 
specified  in  step  12  for  the  number  of  elements 
to  be  changed. 

14 .  The  message 

ANY  MORE  MODS  OR  LIST?  (Y  or  N‘ 

is  printed. 

- Y The  user  wishes  to  list  or  modify  more 
elements . 


Return  to  step  9 . 

« N The  user  is  finished  with  this  bio,'k. 

The  program  proceeds  to  step  15 . 

15.  The  message 

ANOTHER  BLOCK?  (Y  or  N) 
is  printed. 

• Y The  user  wishes  to  list  or  modify  the  values 
of  some  elements  in  another  block. 

Return  to  step  8 . 

- N The  user  is  finished  listing  or  modifying 
values  of  the  elements  in  the  aperture. 

Proceed  to  step  16. 

16 . The  message 

ANY  HOLES?  (Y  or  N> 
is  printed. 


E&l 
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- Y The  user  wishes  to  make  holes  in  the 
aperture . 


Proceed  to  step  17. 


- N The  user  is  finished  with  the  aperture 
changes  . 


Proceed  to  step  21. 


The  message 


ICNTJ,  ICNTK 


is  printed. 


ICNTJ  - The  coordinate  for  the  center  of  the  hole 
in  the  horizontal  directions. 


ICNTK  - The  coordinate  for  the  center  of  the  hole 
in  the  vertical  direction. 


NOTE:  Coordinates  for  the  center  of  the  holes 
are  given  with  respect  to  the  upper  left 
corner  element  increasing  to  the  right  and 
down.  The  upper  left  corner  element  has 
the  coordinates 


(ICNTJ,  ICNTK)  - (1,  1) 


Example:  If  the  center  of  the  hole  were  to  be  in 
the  lower  left  corner  of  a 64  x 64 
element  aperture  the  coordinates  would 
be 


(ICNTJ,  ICNTK)  - (1,  64) 


The  message 


KHOLE 


fa  Mg q jA..- 
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is  printed. 

KHOLE  - Th8  radius  of  the  hole  to  be  'punched*  in 
the  aperture . 

19.  The  hole  is  punched  in  the  aperture  and  the 
following  message  is  printed 

ANOTHER  HOLE?  (Y  or  N> 

The  user  responds  with  either  of  the  following: 

- Y The  user  wishes  to  have  another  hole 

punched . 

Return  to  step  17. 

- N The  user  is  finished  punching  holes  in  the 

specified  aperture. 

20.  The  message 

ANY  ELEMENTS  CHANCED?  (Y  or  N) 
is  printed. 

- Y The  user  wishes  to  list  or  modify  some 

element  values  in  the  specified  aperture. 

Return  to  step  8. 

■ N The  user  is  finished  with  the  aperture 
changes . 

21 . The  message 

ANOTHER  OUTPUT  CENERATED?  (Y  or  N) 
is  printed. 

- Y The  user% wishes  to  generate  another  output 


A_  JSsLtejo.. 
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aperture  from  the  original  input  aperture. 
Proceed  to  step  22. 

- N The  user  is  finished  with  the  original  input 
file. 

Proceed  to  step  23. 

22 . The  message 
DETACH  OUTPUT  FILE  ISTAT  - X 

is  printed,  where  the  detach  is  successful  if  X ■ 
0;  otherwise  X * 1 . 

Return  to  step  5. 

23 . The  message 
ANOTHER  FILE  MODIFIED?  (Y  or  N) 
is  printed. 

- Y The  user  wishes  to  begin  the  program  again 
with  a new  input  file. 

Return  to  step  3 . 

- N The  user  is  finished  with  the  program. 

In  both  cases  the  message 
DETACH  OUTPUT  FILE  ISTAT  - X 
DETACH  INPUT  FILE  ISTAT  - X 

is  printed,  where  if  the  detach  is  successful  X » 
0;  otherwise  X - 1 . 
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IV . SUBPROCRAM  REQUIRED 


None 


V.  RESTRICTIONS.  REQUIREMENTS , AMD  MISC.  DATA 


1.  If  the  input  file  and  the  output  file  are  the 
same,  then  an  affirmative  answer  in  step  21  to 
the  question 

ANOTHER  OUTPUT  CENERATED? 

causes  the  second  generated  output  to  be  a 
modification  of  the  first  output.  This  is  a 
result  of  the  fact  that  the  first  modification  is 
written  over  the  original  aperture. 

2.  If  the  response  to  the  question 
HOV  MANY  ELEMENTS  CHANCED? 

in  step  12  is  less  than  zero  the  program  jumps  to 
step  15.  If  the  response  is  greater  than  100  the 
program  will  return  to  the  beginning  of  step  12. 

3.  This  program  can  be  executed  only  under  Honeywell 
CCOS  TSS  YFORT  subsystem. 
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PROGRAM  PDFESTR 

1 .  PURPOSE 

The  program  PDFESTR  generates  a histogram  of  the 
radiating  elements  in  a statistically  loaded 
aperture.  The  width  of  each  annulus  may  be  varied 
and  the  center  of  the  annuli  may  be  specified. 

II.  PERIPHERAL  DEVICES  REQUIRED 


LUD 

NAME 

USE 

01 

User  PRMFL 

Input  Aperture 

05 

Time  sharing 
keyboard 

terminal 

User  Input 

06 

Time  sharing 
printer 

terminal 

Computer  Output 

III .  OPERATING  PROCEDURES 

1.  Enter  TSS  YFORT 

2.  Type: 

i RUN  PDFESTR 

3.  The  message 
LRJ,  LRK 

’ i 

is  printed. 

LRJ  - Number  of  blocks  in  the  x-direction 
(horizontally) 

LRK  - Number  of  blocks  in  the  y-direction 
(vertically) 
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4 . The  message 
OFSTJ,  OFSTK 


is  printed. 


OFSTJ  - Value  added  to  the  calculated  center  of 
the  aperture  to  offset  the  origin 
horizontally . 


OFSTK  - Value  added  to  the  calculated  center  of 
the  aperture  to  offset  the  origin 
vertically. 


5.  The  message 
INPUT  FILE 


is  printed. 


The  user  should  respond  with  an  existing  PRMFL 
which  contains  a statistically  loaded 
aperture.  Follow  the  entry  with  a semicolon(j) 


Example i 


«/SUBC ATSPSVRD/FILENAMf PSVRD ; 


If  the  file  name  is  improper  then  the  message 
ATTACH  FAILED  ISTAT  - X 


is  printed.  Where  X is  the  first  status  word 
returned  by  the  File  Management  Supervisor  (see 
the  Time  Sharing  System  Programmers  Reference 
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Manual > BR39,  p.(3-39)  or  Table  I,  herein)  or 
will  contain: 


1 - file  is  currently  open 


2 - 


teletypewriter  requested  in 
batch  mode  (illegal) 


3 - 


additional  memory  needed,  request 
denied  (time  sharing  user  will  be 
terminated) 


4 - CATFIL  all  blanks 


NOTE:  See  Honeywell  series  600/6000  Fortran 

manual,  BJ67,  p.  (G-35)  - p.  (6-36)  for 
more  details  on  the  subroutine  ATTACH. 


If  the  file  name  is  unacceptable  the  program  will 
return  to  the  beginning  of  step  5.  Otherwise, 
proceed  to  step  6 . 


6.  The  message 
RINC,  RLIM 


is  printed 


RINC  - The  incremental  radius  or  annuli  width 
used  in  generating  the  histogram  values 


RLIM  - The  maximum  radius  of  interest. 


7.  The  message 
ICON,  HDPACK 
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Is  printed. 


ICON  - Hode  flag 

- -1  The  program  halts 

- 0 Histogram  data  is  normalized  to 

give  a unity  cummulative 
distribution 

- 1 Histogram  data  is  converted  to 

probability  density  estimate  data. 
The  new  data  is  divided  by  the 
product  of  the  cell  width  and  the 
total  number  of  elements. 

- 2 Raw  histogram  data. 

NDPACK  - The  number  of  incremental  radius 

histogram  cells,  RINC,  combined  to  make 
each  output  histogram  cell. 

8.  The  program  repeat?  to  step  7. 

IV.  SUBPROGRAMS  REQUIRED 

None 

V.  RESTRICTIONS.  REQUIREMENTS.  AND  MISC . DATA 

1.  This  estimator  can  only  be  used  on  apertures  that 
have  been  statistically  loaded. 

2.  The  program  TBLS  can  be  used  to  compare  the 
histogram  values  of  an  aperture  with  similar 
values  calculated  for  the  particular  weighting 
function  used  to  load  the  aperture.  This 
comparison  gives  some  degree  of  'goodness'  for 
the  particular  load. 
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PROGRAM  PLARY 


PURPOSE 


The  program  PLARY  loads  a PRMFL  with  a user 
specified  antenna  aperture  illumination  function 


II.  PERIPHERAL  DEVICES  REQUIRED 


LUD 

NAME 

USE 

01 

User  PRMFL 

Aperture  Output 

05 

Time  sharing 
keyboard 

terminal 

User  Input 

06 

Time  sharing 
printer 

terminal 

Computer  Output 

III . OPERATING  PROCEDURES 


1.  Enter  TSS  YFORT 

2.  Type: 

RUN  PLARY 

3 . The  message 
OUTPUT  FILE  NAME 
is  printed. 


The  user  should  respond  with  an  existing  PRMFL 
name  which  is  large  enough  to  store  the  aperture. 
Follow  the  entry  with  a semicolon< ; > . 


Example : 

-/SUBCATSPSWRD/FILENAMSPSWRD ; 
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If  the  file  name  is  improper  then  the  message 


ATTACH  FAILED  ISTAT  - X 


is  printed,  where  X is  the  first  status  word 
returned  by  the  File  Management  Supervisor  (see 
the  Time  Sharing  System  Programmers  Reference 
Manual , BR39,  p.(3-33>  or  Table  1 herein)  or  will 
contain : 


- file  is  currently  open 


teletypewriter  requested  in  batch  mode 
(illegal) 


additional  memory  needed,  request  denied 
(time  sharing  user  will  be  terminated) 


- CATFIL  ail  blanks 


NOTE:  See  Honeywell  series  B00/6000  Fortran 

manual,  BJB7,  p.(G-35)  - p.(6-3B)  for  more 
details  on  the  subroutine  ATTACH. 


If  the  file  name  is  unacceptable  the  program  will 
return  to  the  beginning  of  step  3.  Otherwise 
proceed  to  step  4 . 


4.  The  message 


STATISTICAL  TAPER? 


is  printed.  The  user  responds  either  Y (yes)  or 

N (no). 


YES  The  user  wishes  to  load  the  array  using  a 
space  tapered  or  ’thinned’  loading 
technique.  Proceed  to  step  5. 


PAAS-UM-PLARY 
1 MAY  78 


NO  The  user  wishes  to  load  the  array  using 
element  to  element  amplitude  tapering. 
Proceed  to  step  6. 

5 . The  message 
XKK,  MAD1,  JRND 
is  printed. 

XKK  - This  value  is  the  probability  that  an 

element  will  occur  at  the  normalized  peak 
of  the  design  weighting  function.  (Usually 
XKK  - 1.0) . 

MAD1  - The  starting  address  for  selecting  random 
numbers  from  the  random  number  array 
<is:madi5120>. 

JRND  - Random  number  generator  initialization 
constant  ( 0 < JRND  < 236  - 1). 

6.  The  message 
IAPTFL 

is  printed. 

The  user  responds  with  a number  which  determines 
the  shape  of  the  aperture . 


- 1 circular 


- 2 elliptical 

- 3 rectangular 
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7.  If  IAPTFL  - 1,  the  message 
HEDGE,  XHOLE 
is  printed. 

HEDGE  “ Outside  radius  of  the  aperture. 


XHOLE  - Radius  of  a hole  centered  In  the  aperture 
(For  no  hole,  XHOLE  - 0.0). 


Proceed  to  step  10. 

8.  If  IAPTFL  - 2,  the  message 
NMAJOR,  NMINOR,  XHOLE 
is  printed. 


NMAJOR  - Length  of  the  semi-major  elliptical 
axis  . 


NMINOR  - Length  of  the  semi-minor  elliptical 
axis  . 


XHOLE  - Radius  of  a hole  centered  at  the 

intersection  of  the  major  and  minor  axes 
(For  no  hole,  XHOLE  « 0.0). 


Proceed  to  step  10. 

9.  If  IAPTFL  » 3,  the  message 
NWIDTH,  NHIGH 
is  printed. 

NWIDTH  - Width  of  the  rectangular  aperture, 
IJHICH  " Height  of  the  rectangular  aperture 
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10.  The  message 
IWTFLC 

is  printed. 

The  user  responds  with  a number  which  represents 
the  desired  weighting  function. 

-0  no  weighting  function 

- 1 cosine  to  a power  on  a pedestal 

- 2 Blackman  weighting  function 

- 3 Kaiser  weighting  function 

- 4 triangular  weighting  function 

- 5 Taylor  weighting  function 
• 6 Bessel  weighting  function 
■ 7 cubic  weighting  function 

- 8 Bayliss  weighting  function 

11.  If  IWTFLC  - 0,  proceed  to  step  21. 

12.  If  IAPTFL  - 1,  the  message 
VTRAD 

is  printed. 

The  user  responds  with  the  desired  radius  of  the 
weighting  function. 

Proceed  to  step  14  . 
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13.  If  IAPTFL  * 2 1 or  I APTFL  m 3,  the  message 
ZJRAD , ZKRAD 

is  printed. 

ZJRAD  - Half  the  span  of  the  weighting  function 
in  the  x-direction. 

ZKRAD  - Half  the  span  of  the  weighting  function 
in  the  y-direction . 

NOTE:  These  refer  to  the  elliptical  and 

rectangular  weighting  functions  which  are 
products  of  the  orthogonal  weighting 
functions . 

14.  If  IWTFLC  - 1,  the  message 
WTPED , NWTPOW 

is  printed. 

WTPED  - The  height  of  the  pedestal . 

NWTPOW  - The  power  of  the  cosine. 

Proceed  to  step  21. 

15.  If  IWTFLC  - 2,  proceed  to  step  19. 

1G.  If  IWTFLC  - 3,  the  message 

WKASIR 
is  printed. 

WKASIR  The  Kaiser  variable  for  the  trade-off 
between  mainlobe  width  and  sidelobe 
amplitude . 
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Proceed  to  step  21. 

17.  If  IVTFLG  - 4,  proceed  to  step  21. 

18.  If  IVTFLG  - 5,  the  message 
DB,  NBAR 

is  printed. 

DB  - Sidelobe  in  dB  with  reference  to  the  main 
lobe . 

NBAR  - Number  of  zeros  used  in  approximating  the 
Dolph-Chebyschev  weighting  distribution. 

Proceed  to  step  21. 

19.  If  IVTFLG  - 6,  the  message 
BESCAL,  BESEDC 

is  printed. 

BESCAL  - The  maximum  amplitude  at  the  center  of 
the  aperture  . 

BESEDG  - The  scale  factor  used  in  calculating  the 
argument  for  evaluating  the  Bessel 
function  for  the  actual  radial  location 
on  the  aperture. 

Proceed  to  step  21 . 

2®.  If  IVTFLG  - 7,  the  message 

CUBK,  VTRAD 

is  printed. 
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CUBK  - The  amplitude  scaling  constant. 

WTRAD  - The  half  span  of  the  weighting  function. 

21 . The  message 
NBITS 

is  printed. 

NBITS  - The  number  of  bits  used  to  control  the 
digital  phase  shifters. 

22.  The  message 

ANY  BEAM  STEERING? 
is  printed. 

The  user  responds  with  a Y(yes>  or  an  N(no). 

23.  If  the  response  in  step  22  is  Y the  message 
DELPHJ,  DELPHK 

is  printed. 

DELPHJ  - The  beam  steering  in  degrees  in  the 
x-di rection . 

DELPHK  - The  beam  steering  in  degrees  in  the 
y-direction . 

If  the  response  to  step  22  is  N proceed  to  step 

24. 

24.  The  message 
QUADRATIC  ERROR? 
is  printed. 


,-v  iWi 
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The  user  responds  Y(yes)  or  N(no) . 

25.  If  the  response  to  step  24  is  Y,  the  message 
PHERX , PHERY 


is  printed. 


PHERX  - The  maximum  phase  error  in  degrees  at  the 
edge  in  the  x-direction. 


PHERY  - The  maximum  phase  error  in  degrees  at  the 
edge  in  the  y-direction. 


If  the  response  to  step  24  is  N then  proceed  to 
step  26. 


26  . The  message 

BESSELL  ERROR? 


is  printed. 

The  user  responds  Y(yes)  or  N(no). 

27.  If  the  response  to  26  is  Y,  the  message 
BESERR,  BSCAL 


is  printed. 


BESERR  - The  maximum  phase  error  in  degrees  at 
the  center  of  the  aperture. 


BSCAL 


- The  scaling  factor  used  in  calculating 
the  argument  for  evaluating  the  Bessell 
function  from  the  actual  radial  location 
on  the  aperture. 
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28  . The  message 


LRJ,  LRK 


is  printed. 


LRJ  - The  number  of  blocks  in  the  x-Jirection. 


LRK  - The  number  of  blocks  in  the  y-direction. 


A message  is  now  printed  which  shows  the  user 
what  values  he  assigned  to  the  program  parameters 
and  will  be  used  in  calculating  the  aperture 
illumination  pattern. 


30.  The  message 


DETACH  ISTAT  - K 


Is  printed  after  the  aperture  is  loaded.  X - 0 
if  the  detach  was  successful;  otherwise  it  is  1. 


31.  Program  aperture  terminate.,  and  the  user  is 

returned  to  build  mode  under  TSS  YFORT.  This  is 
Indicated  by  an  asterisk  (*). 


IV . SUBPROCRAMS  REQUIRED 


EXPND 

BESS 

CAM 

WEIGHT 


V.  RESTRICTIONS,  REQUIREMENTS.  AMD  MISC . DATA 


1.  If  the  weighting  function  span  is  less  than  the 
aperture  span  then  those  elements  outside  the 
weighting  function  span  are  set  to  zero. 


.fta. WiBf.  * 
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2.  The  Kaiser  variable,  WKASIR,  should  be  within  the 
range 


2 £ WKASIR  £ 1© 


3. 


The  number  of  zeros,  NBAR,  in  the  Taylor 
approximation  must  be  in  the  range 


3 £ NBAR  £ 2© 


4.  LRJ  and  LRK  must  be  even  or  the  message  LRJ,LRK 
will  be  retyped  and  the  user  must  respond  with 
two  even  numbers . 


S.  It  should  be  noted  that  if  the  user  has  specified 
a rather  complex  aperture  illumination  pattern 
(Example:  a large  circular  aperture  (XEDCE  - 
40.0)  with  a Taylor  distribution  (NBAR  - 20))  the 
time  required  to  load  may  be  quite  long  (perhaps 
5 minutes) . 


G.  This  program  can  be  executed  only  under  Honeywell 
CCOS  TSS  YFORT  subsystem. 
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PROCRAM  PLTDVR 


I . PURPOSE 


The  program  PLTDVR  formats  either  an  aperture 
illumination  or  a far-f.eld  energy  distribution  that  is 
stored  in  a permanent  di  -k  file  (PRMFL)  for  making 
pseudo-3d  plots  using  tht  DUIS . The  data  is  formatted  and 
transmitted  to  the  DUIS  for  recording  and  subsequent 
production  of  3d  plots. 


II.  PERIPHERAL  DEVICES  REQUIRED 


LUD 

NAME 

01 

Input  PRMFL 

03 

Time  sharing 

keyboard 

06 

Time  sharing 

printer 

OPERATING 

PROCEDURE 

USE 


Aperture  or 
far-field  input 


Load  RJE  program  into  DUIS  .and  start  program 
execution.  See  User's  Manual  on  RJE-300  or  RJE-1200. 


1.  Enter  TSS  YFORT 

2 . Type 

RUN  PLTDVR 
3 . The  message 

INPUT  FILE  NAME? 


is  printed.  The  user  should  respond  with  the 
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name  of  an  existing  file  which  contains  the  data 
to  be  formatted  and  recorded.  The  file  name  must 
be  followed  with  a semicolonC ; ) . If  the  user 
types  STOP  the  program  will  terminate  execution 
and  the  user  is  returned  to  build  mode  under  TSS 
YFORT.  This  is  indicated  by  an  asterisk(»). 

4.  If  the  file  name  is  not  acceptable  the  message 

UNSUCCESSFUL  ATTACH  ISTAT  * X 

is  printed,  where  X is  the  first  status  word 
returned  by  the  File  Management  Supervisor  (see 
the  Time  Sharing  System  Programmers  Reference 
Manual , BR39,  p.(3-39),  or  Table  1 herein)  or 
will  oontain: 


1 - file  is  currently  open 

2 - teletypewriter  requested  in  batch  mode 

(illegal) 

3 - additional  memory  needed,  request  denied 

(time  sharing  user  will  be  terminated) 

4 - CATFIL  all  blanks 

5.  If  the  file  name  in  step  3 is  accepted  the 
message 

JVID,  JWIDSP,  NBMAX 
is  printed. 

JWID  - The  width  in  blocks  of  the  part  of  the 
array  to  be  formatted  and  transferred 

JWIDSP  - The  width  in  blocks  of  the  total  input 
array 
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NBMAX  - The  logical  record  number  of  the  last 
block  to  be  transferred. 

If  JVIDSP  is  greater  than  10  the  value  of  JWIDSP 
is  set  to  10  and  the  following  warning  message  is 
printed: 

JWIDSP  > 10... SET  TO  10 

If  JWIDSP  is  greater  than  JWID,  the  value  of 
JWIDSP  is  set  to  JWID  and  the  following  warning 
message  is  printed: 

JWIDSP  > JWID... SET  JWIDSP- JWID 

6 . The  message 

Enter  FIRST  LREC,  ISDN 
is  printed. 

FIRST  LREC  - The  logical  record  .of  the  block  in 

the  upper  left  corner  of  the  part  of 
the  array  to  be  formatted  and 
transferred.  If  the  response  is  -1 
the  program  transfers  to  step  8. 

ISGN  - - 0 data  is  processed  for  magnitude 

plotting 

- 1 real  component  is  processed  for 
bipolar  plotting. 

7.  The  data  transmission  to  the  DUIS  now  begins. 

When  all  of  the  requested  data  has  been 
transmitted  the  Tektronix  display  will  beep  once 
and  erase.  The  user  should  then  press  carriage 
return  and  the  following  message  is  printed: 


TH  - X 
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where  X is  the  maximum  value  of  all  the 
transmitted  data.  This  value  is  used  in  scaling 
the  plots  produced  using  the  DUXS  subprogram 
H3DPL . The  program  now  returns  to  step  6. 

8.  The  message 

DETACH  ISTAT  - X 

is  printed,  where  X - 0 if  the  detach  is 
successful  and  X - 1 otherwise.  The  program 
returns  to  step  3. 

SUBPROCRAMS  REQUIRED 

PLOTD 

CHOP 

RESTRICTIONS.  REQUIREMENTS , AMD  HISC . DATA 

1.  This  program  can  be  executed  only  under  Honeywell 
CCOS  TSS  YFORT  subsystem. 

2.  The  transmission  time  for  a 64  x 64  point  plot  is 
about  10-15  minutes  at  300  baud.  This  should  be 
kept  in  mind  if  larger  plo*s  are  to  be  attempted. 

3.  If  the  response  to  FIPST  LREC  is  greater  than  the 
response  to  NBMAX  the  program  will  repeat  step  6. 

4.  The  magnitude  of  the  main  lobe  determines  the 
value  of  the  least  significant  bit  (LSB) . If  the 
ratio  of  mainlobe  to  sidelobe  level  exceeds  the 
accuracy  of  the  12  bit  word  length,  errors  in 
magnitude  representat i on  result.  The  dynamic 
range  of  the  12  bit  word  is  72  dB  (8  dB/bit). 
Magnitude  errors  begin  to  appear  if  the  sidelobes 
are  more  than  80  dB  down . 
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PROCRAM  RNDERR 


PURPOSE 

The  RNDERR  program  modifies  an  antenna  aperture 
distribution  (resident  in  PRMFL)  by  adding  a phase 
error  distribution  specified  by  the  user.  The 
resulting  aperture  is  stored  in  another  PRMFL  where 
it  can  be  accessed  for  further  processing. 


PERIPHERAL 

DEVICES  REQUIRED 

LUD 

NAME 

USE 

01 

User  PRMFL 

Input  Aperture 

02 

User  PRMFL 

Output  Aperture 

05 

TSS  terminal  keyboard 

User  Input 

06 

TSS  terminal  printer 

Computer  Output 

OPERATING 

PROCEDURE 

1.  Enter  TSS  YFORT 

2 . Type 

RUN  RNDERR 

3.  The  message 

NTYPE,  MAD1 , JRND , LRJ,  LRK 
is  printed. 

NTYPE  - Determines  the  type  of  error  distribution 
to  be  added  to  the  aperture. 
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- 1 uniform  distribution 
“ 3 Gaussian  distribution 

MAD1  - The  starting  address  for  selecting  random 
numbers  from  the  random  number  array. 
(l<iMADl£128> 

JRND  - Random  number  generator  initialization 
constant.  (0<JRNDS236  - 1) 

LRJ  - The  number  of  blocks  in  the  x-directlon 
of  the  input  aperture. 

LRK  - The  number  of  blocks  in  the  y-direction 
of  the  input  aperture. 

4.  If  NTYPE  - 1,  the  message 

UMEAN,  UUEXT 

is  printed. 

UMEAN  - Mean  value  of  the  uniform  distribution  in 
degrees . 

UUEXT  - Width  of  the  uniform  distribution  in 
degrees . 

Proceed  to  step  6. 

5.  If  NTYPE  - 3,  the  message 

XMEAN , SIGMA 

is  printed. 

XMEAN  - Mean  value  of  the  Gaussian  distribution 
in  degrees. 

SIGMA  - Standard  deviation  of  the  Gaussian 

42 


$i 


* 1 
i -a® 


j 


, 'if 
i 3| 

m 


l I 


wMJ SmitX  it  JtoL  3tk£tt2i 


PAAS-UH-RNDEFiR 
6 APR  78 


distribution  in  degrees. 

S . The  message 

INPUT  FILE  NAME 

is  printed.  The  user  should  respond  with  an 
existing  PRMFL  name  which  has  the  input  aperture 
stored  in  it  followed  by  a semicolon( ; ) . If  the 
file  name  is  not  acceptable  the  following  message 
will  be  printed: 

UNSUCCESSFUL  ATTACH  ISTAT  - X 

where  X is  the  first  status  word  returned  by  the 
File  Management  Supervisor  (see  the  Time  Sharing 
System  Programmers  Reference  Manual,  BR39, 
p . (3-39 ) or  Table  1 herein)  or  will  contain: 

1 - file  is  currently  open 

2 - teletypewriter  requested  in  batch  mode 

(illegal) 

3 - additional  memory  needed,  request  denied 

(time  sharing  user  will  be  terminated) 

4 - CATFIL  ail  blanks 


NOTE:  See  Honeywell  Series  600/6000  Fortran 

manual,  BJ67,  p.(6-35)  - (6-36),  for  more 
details  on  the  subroutine  ATTACH. 

If  the  input  file  name  is  unacceptable  the 
program  returns  to  the  beginning  of  step  6. 
Otherwise,  proceed  to  step  7. 
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7.  If  the  PRMFL  name  input  in  step  B is  acceptable 
the  message 

OUTPUT  FILE  NAME 

is  printed.  The  user  should  respond  with  an 
existing  PRMFL  name  in  which  the  output  aperture 
is  to  be  stored,  followed  by  a semicolonC ; ) . If 
the  PRMFL  name  is  not  acceptable  the  same 
procedure  as  described  in  step  B applies  here  and 
the  program  will  return  to  the  beginning  of  step 
7.  Otherwise,  proceed  to  step  8. 

8.  The  antenna  aperture  distribution  is  modified  and 
stored  in  the  specified  output  PRMFL. 

9.  When  the  processing  is  finished  the  message 
DETACH  ISTAT  - X 

is  printed  twice.  The  first  time  is  for  the 
input  PRMFL  and  the  second  is  for  the  output 
PRMFL.  If  X - 0 then  the  detach  is  successful, 
otherwise  X - 1 . 

10.  Program  execution  terminates  and  the  user  is 

returned  to  build  mode  under  TSS  YFORT.  This  is 
indicated  by  an  asteriskC*). 


SUBPROGRAMS  REQUIRED 


RRAND 
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V . RESTRICTIONS,  REQUIREMENTS,  AND  MISC . DATA 

1.  If  MAD1  Is  outside  of  the  range  1 < MAD1  £ 128 
the  computer  will  print  the  message  of  step  3 
again  and  the  user  must  respond  correctly. 

2.  This  program  can  be  executed  only  under  Honeywell 
GCOS  TSS  YFORT  subsystems. 

3.  The  input  PRMFL  and  the  output  PRMFL  can  be  the 
same  file. 


45  / 46 


PAAS-UM-RTI4 
1 MAY  78 


PROCRAM  RTI4 


PURPOSE 


The  program  RTI4  is  a subroutine  which  converts 
data  amplitude  to  letters  of  the  alphabet.  The 
resulting  characters  are  organized  into  a matrix 
which  is  printed  on  a remote  terminal.  Range  is 
usually  shown  vertically  and  cross  range  is  displayed 
horizontally  with  intensity  displayed  by  the 
character  placed  in  the  cell,  e .g . 0 represents  0 dB 
with  respect  to  the  reference,  A represents  -10  dB 
with  respect  to  the  reference,  Z represents  -38  dB 
with  respect  to  the  reference.  This  subroutine 
processes  the  same  type  of  data  as  PLTDVR  but  does  so 
in  a much  mere  compact  form. 


II.  PERIPHERAL  DEVICES  REQUIRED 


LUD 


01 


05 


08 


NAME 


USE 


Input  PRMFL  (permanent  Aperture 
disk  file)  far-field  input 

Time  sharing  terminal  User  Input 
keyboard 

Time  sharing  terminal  Computer  Output 
printer 


III.  OPERATING  PROCEDURE 


1.  Enter  TSS  YFORT 


2 . Type 

RUN  RTI4 
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3 . The  message 

INPUT  DESIRED  FILE  NAME 
is  printed. 

The  user  should  respond  with  the  name  of  an 
existing  PRMFL  vhlch  contains  the  data  to  be 
displayed.  The  file  name  must  be  followed  by  a 
semieolon< ; ) . If  the  user  types  STOP  the  program 
jumps  to  step  9. 

4.  If  the  PRMFL  name  is  not  acceptable  the  following 
message  is  printed: 

UNSUCCESSFUL  ATTACH  ISTAT  - K 

where  X is  the  first  status  word  returned  by  the 
File  Management  Supervisor  (see  the  Time  Sharing 
System  Programmers  Reference  Manual , BR39, 
p.(3-39)  or  Table  1,  herein)  or  will  contain: 


1 - file  is  currently  open 

2 ■ teletypewriter  requested  in  batch  mode 

(illegal) 

3 * additional  memory  needed,  request  denied 

(time  sharing  user  will  be  terminated) 

4 - CATFIL  all  blanks 

NOTE:  See  Honeywell  Series  600/6000  Fortran 

manual,  BJ67,  p.(6-35)  - (6-36)  for  more 
details  on  the  subroutine  ATTACH.  After 
the  messege  s printed  the  program  returns 
to  step  3 and  begins  again. 
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5.  If  the  PRMFL  name  in  step  3 is  accepted  the 

message 

FLOOR,  YINC,  JWID,  NBHAX 

is  printed. 

FLOOR  - The  reference  in  dB  below  which 

everything  is  represented  in  dashes  (-) 
on  the  RTI  plot.  All  the  data  with  a 
greater  value  is  represented  by  a letter, 
number,  or  punctuation  symbol  and  the 
value  relative  to  the  FLOOR  is  calculated 
using  YINC. 

YINC  - The  increment  in  dB  between  each 

successive  letter,  number,  or  symbol  is 
determined  by  the  value  assigned  to  YINC. 
See  Table  2. 

JWID  - The  width  in  blocks  across  the  whole  side 
of  the  data  array  as  determined  by  the 
parameter  LRJWID  in  the  program  FFT2DX  or 
LRJ  in  PLARY,  RNDERR , or  FILMOD. 

NBMAX  - The  number  of  the  last  block  to  be 

displayed.  The  final  display  will  be  4 
blocks  vide  and  as  long  as  the  user 
chooses  depending  on  the  value  of  NBMAX. 

NOTE:  Two  examples  of  the  relative  dB  values  for 
a specified  FLOOR  and  YINC  are  shown  in 
Table  2. 

8 . The  message 

Enter  FIRST  LREC,  DISPLAY  WIDTH 
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is  printed.  The  user  should  respond  with  the 
number  of  the  block  that  he  wishes  to  be  placed 
in  the  upper  left  hand  corner  of  the  RTI  and  the 
width  of  the  RTI  in  samples.  The  maximum  number 
of  samples  is  128.  If  the  user  types  -1  then  the 
program  goes  to  step  8. 

7.  The  computer  now  begins  to  transmit  the  RTI 
display  to  the  DUIS  or  a time-sharing  terminal 
printer  and  will  continue  until  it  reaches  the 
specified  stopping  block  number,  NBMAH. 

8 . The  message 
DETACH  ISTAT  - X 

is  printed  where  X - 0 if  the  detach  is 
successful  and  X - 1 otherwise.  Co  to  step  3. 

9.  An  asterisk  (*)  is  printed  and  the  program  is 
finished . 

IV.  SUBPROGRAMS  REQUIRED 

None 

V.  RESTRICTIONS,  REQUIREMENTS,  AND  MISC . DATA 

1.  This  program  can  be  executed  only  under  Honeywell 
CCOS  TSS  YFORT  subsystems. 

2.  If  the  response  to 
INPUT  LREC  DESIRED 


is  greater  than  NBMAX  then  the  program  will 
repeat  step  6. 
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3.  Any  dB  level  below  the  value  of  FLOOR  will  be 
represented  by  a dash  (-)  . Any  dB  level  abovs 
the  value  calculated  for  'S'  will  be  represented 
by  a . 
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TABLE  2 

EXAMPLES  OF  RELATIVE  VALUE  OF  LETTERS,  NUMBERS,  AND  SYMBOLS 


FLOOR  - -20 .0 

20.0 

t 

30.0 

FLOOR  - 10.0 

YINC  - 1.0 

19  .0 

* 

29  .5 

YINC  - 0.3 

18.0 

* 

29  .0 

17  .0 

+ 

28.3 

16  .0 

0 

28.0 

15  .0 

1 

27  .3 

14.0 

2 

27.0 

13.0 

3 

26  .5 

12  .0 

4 

26  .0 

11.0 

5 

23.3 

10.0 

6 

23  .0 

9.0 

7 

24.3 

8.0 

8 

24.0 

7 .0 

9 

23.5 

6.0 

A 

23.0 

3.0 

B 

22.5 

4.0 

C 

22.0 

3.0 

D 

21.5 

2 .0 

E 

21 .0 

1.0 

F 

20.5 

0.0 

u 

20.0 

~1,0 

H 

19.5 

-2,0 

I 

19.0 

•*3  .0 

J 

10  .5 

dB  level 

-•4  .0 

K 

18.0 

dB  level 

-5  .0 

L 

17  .3 

-6  .0 

M 

17  .0 

-7.0 

N 

16.3 

-8  .0 

0 

16  .0 

-9  .0 

P 

15 .5 

-10.0 

Q 

15  .0 

-11  .0 

R 

14  .5 

-12.0 

S 

14  .0 

-13  .0 

T 

13  .5 

-14  .0 

U 

13  .0 

-15  .0 

V 

12.5 
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-IS  .0 

W 

12  .0 

-17  .0 

X 

11 .5 

-18.0 

Y 

11 .0 

-19.0 

Z 

10.3 

-20.0 

- 

10.0 

53/54 


■■ 


v1J*W y^jM*^ 


SkSsss* 


PROCRAM  TBLS 


I .  PURPOSE 


PAAS~Ull-TBLS 
6 APR  78 


The  program  TBLS  computes  and  tabulates  sample 
values  of  selected  weighting  functions.  The  program 
also  generates  data  which  is  used  in  checking  the 
probability  density  function  of  space  tapered  arrays, 


The  program  may  be  used  in  one  of  the  following 
three  modes: 

1.  TBLS  generates  the  value  of  the  weights  for 
specific  distribution  width  and  weighting 
function  type. 

2.  For  the  Taylor  and  Bayllss  functions,  the 
following  three  modes  are  available  (3).: 

(a)  For  a specific  distribution  width  TBLS 
generates  all  the  sample  weights  for 
sidelobe  levels  from  20  to  90  dB  in  steps 
of  5 dB  with  n ranging  from  3 to  20. 

(t)  For  a specific  distribution  width  and  dB 
level,  TBLS  generates  the  sample  weights 
with  n ranging  from  3 to  20. 

(c)  For  a specific  distribution  width,  dB 
level,  and  n,  TBLS  generates  the  sample 
weighting  function. 

This  program  produces  tables  similar  to  Hansen's 
(1)  but  with  more  flexibility,  greater  accuracy,  and 
greater  range  in  dB  and  n.  Tables  may  be  generated 
f'or  the  Bayliss  as  well  as  the  Taylor  di stri bution . 

3.  For  all  of  the  above  modes  TBLS  generates  data 
which  is  either  in  the  form  cf*  standard  weights 
or  in  a form  which  may  be  used  to  compare  with 
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the  data  generated  by  PDFESTR  to  check  the 
density  function  of  space  tapered  arrays. 


PERIPHERAL  DEVICES  REQUIRED 


NAME 


TSS  terminal  keyboard  User  Input 


TSS  terminal  printer  Computer  Output 


III.  OPERATING  PROCEDURE 


i.  Enter  TSS  YFORT 


2 . Type 


RUN  TBLS 


3 . The  message 


PDFESTR  DATA  OR  TABLES?  (0  or  1) 


is  printed. 


“ 0 The  user  wishes  to  generate  data  that  will 
be  compatible  with  the  data  from  PDFESTR. 


- 1 The  user  wishes  to  generate  tables  of 
weighting  function  values. 


The  message 


NTYPE 


is  printed. 


The  user  should  respond  with  a number  which 
determines  the  weighting  function  type. 
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- 1 cosine  on  a pedestal  to  a power 

(Proceed  to  5) 

- 2 Blackman 

(Proceed  to  6) 

- 3 Kaiser 

(Proceed  to  7) 

- 4 Bartlett  or  triangular 

(Proceed  to  6) 

- 3 Taylor 

(Proceed  to  10) 

■ 6 Bessel 

(Proceed  to  8) 

- 7 Cubic 

(Proceed  to  9) 

- 8 Bayllss 

(Proceed  to  10) 

The  message 

WTPED . NWTPOV,  IRAD,  WTRAD 
is  printed. 

WIPED  - The  height  of  the  pedestal 

NWTPOW  - The  power  of  the  cosine  function 
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IRAD 


- The  radius  (or  half  span  of  a linear 
array),  in  units  of  elements,  of  the 
array . 


WTRAD  - The  radius  (or  half  span  of  a linear 
array)  of  the  weighting  function. 


NOTE:  For  all  subsequent  entries,  IRAD  and  WTRAD 
have  the  same  meaning  as  above. 


Proceed  to  15 . 


6 . The  message 
IRAD,  WTRAD 


is  printed. 


Proceed  to  15 


7 . The  message 

WKASIR,  IRAD,  WTRAD 


is  printed. 


WKASIR  - The  Kaiser  variable  for  the  trade-off 
between  main  lobe  width  and  sldelobe 
ampl i tude . 


Proceed  to  15 . 


8 . The  message 

BESEDC,  BESCAL,  IRAD,  WTRAD 


is  printed. 
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BESEDG  - The  scale  factor  used  In  calculating  the 
argument  used  in  evaluating  the  Bessel 
function  from  the  actual  radial  location 
on  the  aperture . 

BESCAL  - The  maximum  amplitude  at  the  center  of 
the  aperture. 

Proceed  to  15 . 

9 . The  message 

CUBK,  IRAD,  VTRAD 

Is  printed. 

CUBK  - The  amplitude  scaling  constant 
Proceed  to  15. 

10.  The  message 
IRAD,  VTRAD 
is  printed. 

* 0 The  user  wishes  to  generate  the  complete  set 
of  tables  for  IRAD  with  NBAR  and  DB  varied. 

3 < NBAR  < 20  (19  for  the  Bayllss) 

and  20  < DB  £ 80 

in  steps  of  5 dB . 

Proceed  to  15. 

- 1 The  user  wishes  to  choose  one  dB  level  of 
interest.  Proceed  to  12. 
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The  message 
IDB 

is  printed.  The  user  should  respond  with  an 
integer  value  for  the  specified  dB  level. 

The  message 

ALL  OR  SINCLE  NBAR?  (0  or  1) 
is  printed. 

- 0 The  user  wishes  to  generate  for  the 

specified  IRAD  and  dB  level  all  possible 
NBAR  distributions  in  the  range 

3 £ NBAR  £ 20  (19  for  the  Bayliss) 

Proceed  to  15. 

- 1 The  user  wishes  to  choose  one  NBAR  of 

interest . 

Proceed  to  14. 

The  message 


m 


NBAR 

is  printed.  The  user  should  respond  with  an 
integer  number  for  the  specific  value  of  NBAR 
desired.  If  NBAR  is  too  small,  as  explained  by 
Taylor  (2)  and  Hansen  (1),  the  message 

INVALID  VALUE  FOR  NBAR 

is  printed  and  the  program  repeats  step  14. 

60 


tess&tftst 


Jirxitut'jMi  f'.~r  aua&iUttatfifoirf&aQ  outa 


i 

{ 

j 

I 

I 

I 


i 

i 


I 


I 

:& 


1 


48 

4 


PAAS-UM-TBLS 
B APR  73 


15.  The  program  now  generates  the  appropriate 

distribution  and  prints  the  tables  on  th9  TSS 
terminal  printing  device.  The  program  returns  to 
step  4 vhen  the  requested  distribution  has  been 
printed . 

IV.  SUBPROGRAMS  REQUIRED 

EXPND 

BESS 

GAM 

WEIGHT 

V.  RESTRICTIONS,  REQUIREMENTS.  AND  MISC . DATA 


1.  The  range  of  NBAR  must  be 
3 S NBAR  S 20 
for  the  Taylor  and 
3 S NBAR  S 19 
for  the  Bayl i ss . 

Reference 

1,  Hansen,  R.C.,  'Tables  of  Taylor  Distributions  for 

Circular  Aperture  Antennas, * IEEE  Trans . 

Anten.  Prop.,  AP-8,  #1,  (1/60),  pp.  23-26. 

2.  Taylor,  T.  T.,  ’Design  of  Circular  Apertures  for 

Narrow  Beamwidth  and  Low  Sidelobes,'  The 

Bell  System  Technical  Journal,  Vcl  . 47,  No. 

5,  (May-June  1968),  pp.  623~S51. 
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APPENDIX 


COMPUTER  PROGRAM  DOCUMENTATION 


In  this  appendix  the  computer  program  documentation 
(CPD)  for  the  PAAS  modules  is  presented.  The  CPD's  are 
arranged  in  alphabetical  order  and  each  is  self-contained. 
The  subprograms  BESS  and  GAM  were  obtained  from  RADSIM  for 
use  in  PAAS.  For  the  reader's  convenience  and  completeness 
of  this  document  these  have  been  incorporated  into  this 
appendix . 


Unless  otherwise  stated,  all  software  documented 
herein  is  stored  under  user  master  catalog  'BECAGD01*  in 
the  RADC  H6180  GCOS  system.  The  source  code  for  all 
programs  herein  are  stored  in  PRMFL's  having  the  same  name 
as  the  program,  e.g.  the  source  code  of  the  program  PLARY 
and  all  required  subroutines  is  stored  in  a PRMFL  having 
the  name  PLARY  under  user  master  catalog  BECAGD01 . 

The  documentation  for  each  PAAS  program  presented  in 
this  appendix  is  divided  into  six  sections.  The  order  and 
title  of  each  section  is  as  follows: 


1.  Purpose 


2.  Input  Parameters 

3.  Restrictions,  Requirements,  Miscellaneous  data 


4.  Subprograms  Required 

5.  Theory  of  Operation 

6.  FORTRAN  Listing 


The  content  of  each  section  is  explained  in  the  following 
paragraphs . 
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Purpose 


This  section  contains  a brief  description  of 
the  purpose  of  the  software  module. 

2 . Input  Parameters 

This  section  lists  all  of  the  input 
parameters  for  the  particular  software  module. 

Both  required  and  optional  input  parameters  are 
listed.  Each  parameter  entry  is  broken  into  four 
groups  of  information  and  placed  into  columns  for 
easy  reference.  The  first  column  contains  the 
parameter  name  as  it  appears  in  the  software.  The 
second  column  tells  whether  the  parameter  is 
required  or  optional  . An  R in  the  second  column 
denotes  a required  parameter  while  an  0 denotes  an 
optional  parameter.  The  variable  type,  either 
integer  or  floating  point,  is  noted  in  the  third 
column.  An  I denotes  and  Integer  type  while  an  F 
denotes  a floating  point  variable.  The  fourth 
column  contains  a brief  description  of  the 
parameter  and  how  it  is  used  in  the  program. 

3 . Restrictions,  Requirements,  Miscellaneous  Data 

In  this  section  special  notes  concerning  the 
input  parameters,  use  of  the  program,  potential 
usage  problems,  etc.  are  discussed. 

4 . Subprograms  Required 

In  this  section  the  subprograms  required  for 
the  PAAS  program  are  listed.  Both  subroutine  and 
function  subprograms  are  included. 
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Theory  of  Operation 


In  this  section  the  theory  behind  the 
programming  is  discussed  using  the  variable  names 
and  notation  as  they  appear  in  the  program.  This 
helps  the  user  in  understanding  the  operation  of 
the  program. 


6 • FORTRAN  Source  Code  Listing 


This  section  contains  a listing  of  the 
FORTRAN  source  code. 


The  documentation  for  each  PAAS  subprogram  presented 
in  this  appendix  is  divided  into  seven  sections.  The  order 
and  title  of  each  section  is  as  follows: 


1.  Purpose 


2.  Input  Parameters 


3.  Calling  Sequence 

4.  Restrictions,  Requirements,  Miscellaneous  Data 

5.  Subprograms  Required 

6.  Theory  of  Operation 

7.  FORTRAN  Source  Code  Listing 


These  are  the  same  as  those  that  were  previously  described 
for  the  PAAS  program  CPD's  except  that  the  section  entitled 
‘Calling  Sequence*  has  been  added.  A description  of  that 
section  follows: 
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3 . Calling  Sequence 

This  section  presents  an  example  of  a typical 
FORTRAN  call  for  the  module.  For  function 
subprograms  the  example  calling  sequence  is  shown 
as  an  assign  statement  but  of  course  the  function 
reference  can  be  embedded  in  a FORTRAN  arithmetic 
statement . 
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FUNCTION  BESS 


1.  MODULE  IDENTIFICATION 


RADSIM-CPD-BESS 
1 MAY  78 


Name 


BESS 


Classification  Code 


Subordinate 


Reference  Number 


Not  User  Referenced 


2.  PURPOSE 

This  function  is  used  to  compute  the  value  of  a 
Bessel  function. 

3.  INPUT  PARAMETERS 


Name 


O/R  T 


Description 

The  order  of  the  Bessel  function 

The  argument  of  the  Bessel 
function . 


4.  CALLING  SEQUENCE 
BS  - BESS(0,Z) 

Where:  0,Z  are  the  Input  arguments 

BS  contains  the  computed  value  of  the 
Bessel  function 

5 . RESTRICTIONS,  REQUIREMENTS,  MISCELLANEOUS  DATA 

a.  This  subprogram  was  obtained  from  the  Computer 
Program  Documentation  for  AF  Contract 
F30602-S7-C-0074 . 
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6 . FORTRAN  LISTING 
C 

c 

FUNCTION  BESS(0, Z) 

DOUBLE  PRECISION  A1,A2,BS,ADD,GAM,C 
C 

IF(Z.NE.0.0>  CO  TO  40 
IF(0 .EQ .0 .0)  BS-1.0 
IF(O.NE.0.0)  BS-0.0 
GO  TO  10w 
40  SKALL-1.0E-8 
IF(O) 100, 31, 32 

31  BS-1.0 
AKV-0.0 
Al-1 .0 

52  AKV-AKV+1 .0 

Al-Al *(-l . 0 ) * ( Z/2 . 0 ) * *2/( AKV  * AKV ) 
BS-BS+A1 

IF ( DABS ( Al/ES ) -SMALL  >51,52,52 
51  CO  TO  100 

32  A-0 
N-0 

13  IF(A-1 .0)10, 12 , 12 
12  A-A-1.0 
N-N+l 
CO  TO  13 
10  ARC-A 

G-GAM(ARG) 

A2-1.0 


IF(N)100,76,75 

75  DO  26  NV-l.N 
AJ-NV-1 

IF(A2 .CT.1.0E-38)  CO  TO  28 

BS-0.0 

CO  TO  100 

26  A2-A2*(Z/2 .0)/<O-AJ> 

76  A2-A2  *(Z/2 .0) **ARG/G 
BS-A2 

AKV-0 .0 
Al-1  .0 
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21  AKV-AKV+1.0 

A1-A1*<~1 .0) *(Z/2 .0) **2/(AKV*(AKV+Q) ) 

ADD-A1 *A2 

BS-BS+ADD 

IFCBS ,EQ  .0 .0)  CO  TO  21 
TEST-DABS (ADD/BS) 

IF (TEST-SM ALL >100, 21,21 
100  BESS-BS 
RETURN 
END 
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SUBROUTINE  CHOP 

1.  PURPOSE 

The  purpose  of  this  subroutine  is  to  convert  an 
integer  number  into  two  ASCII  characters.  The  ASCII 
characters  are  packed,  right-adjusted,  into  an  output 
word . 

2.  INPUT  PARAMETERS 

Name  Q/R  T Description 

IDAT  R I The  integer  word  to  be  processed 

IOUT  R I The  output  word  containing  two 

ASCII  characters. 

3.  CALLINC  SEQUENCE 
CALL  CHOP (IDAT, IOUT) 

4.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  The  input  word,  IDAT,  must  be  in  the  following 
range: 

-211  < IDAT  < 2li-l 

5.  SUBPROCRAMS  REQUIRED 
None 

S.  THEORY  OF  OPERATION 

The  input  word,  IDAT,  contains  at  most  12 
significant  bits.  The  rightmost  12  bits  are  extracted 
from  IDAT  and  separated  into  two  6 bit  characters,  II 
and  12.  The  characters  II  and  12  have  values  which 
range  from  0 to  63  and  include  the  ASCII  control 
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character  region  from  0 to  31.  In  order  to  ensure  that 
these  characters  cannot  have  values  in  the  control 
character  region,  the  number  32  is  added  to  each.  If 
this  is  not  done,  problems  arise  with  the  H6180  TSS 
processing.  These  characters  are  packed, 
right-adjusted,  into  the  output  word,  IOUT,  and  control 
returns  to  the  calling  (sub)program . 
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SUBROUTINE  CZFFT 


PURPOSE 

This  subroutine  performs  the  inverse  discrete 
Fourier  transform  of  a sequence  of  input  data  samples. 
The  mechanization  is  based  on  the  Fast  Fourier 
Transform  (FFT)  algorithm  developed  by  Langdon  and 
Sande  from  the  approach  of  J.  W.  Tukey  and  J.  Cooley. 
The  subroutine  described  herein  has  been  structured  to 
facilitate  the  efficient  computation  of  2-dimenslonal 
discrete  inverse  Fourier  transforms. 


INPUT  PARAMETERS 


Name  O/R  T 


Description 


R I 


IOFST  R I 


IHOP  R I 


Power  of  2 which  determines  the 
total  number  of  points  (NTHPOW) 
transformed  by  the  FFT 

Offset  of  the  first  sample  to  be 
transformed  from  the  front  of 
the  array  S 

Power  of  2 which  determines  the 
spacing  between  the  samples  (NHOP) 
to  be  transformed. 


CALLING  SEQUENCES 


CALL  CZFFT (S,N2, IOFST, IHOP) 


Where:  S is  a complex  array  containing  the  data  to  be 

processed.  The  output  samples  are  placed  into 
the  array  S also. 
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4 . RESTRICTIONS  > REQUIREMENTS,  MISCELLANEOUS  DATA 


a.  The  maximum  value  of  N2  is  11,  which  gives  2048 
samples . 


b.  In  order  to  minimize  the  CPU  time  requirements  of 
this  subroutine,  a complex  exponential  look  up 
table  is  used. 


c.  Source  PRMFL: 


BECAVU01/SUPORTSJR/SCZFFT 


d.  Object  PRMFL: 

BECAVUOl/SUPORTtJR/OCZFFT 


5.  SUBPROGRAMS  REQUIRED 


COS 

SIN 


6.  THEORY  OF  OPERATION 


Refer  to  RADSIM-CPD-ZFFT . 
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7.  FORTRAN  LISTING 


SUBROUTINE  CZFFT(S , N2 , IOFST, IHOP) 

C 

C S » /SUPORT/SCZFFT 

C 0 » /SUPORT/OCZFFT 

C 

DATA  PI2/6 .28318531/ 

DATA  IFLAG/0/, NHOPO/- 1/ 

COMPLEX  S ( 1 ) , C , Cl , C2 , C3 , C4 , XEXP ( 153B ) 

REAL  I,I1,I2,I3,I4,RX(2) 

INTECER  PASS,  SEQLOC,  L(15) 

EQUIVALENCE  (J,JI),<PASS,JS>,< NXTLTH , J7 ) , 

« ( LENGTH , J 8 ) , ( SEQLOC , J9 ) , ( ISCALE , J10) , 

* (IARG, Jll ) , (A1 , J12) , (RX(1) ,14 , C) , 

* <RX<2),R4) 

EQUIVALENCE  (L15 ,L<1) ) , (L14 ,L(2) ) , <L13,L<3>> , 

* (L12,L<4)),(L11,L(5)),(L10,L(6)), 

* (L9,L(7)),(L8,L(8)),(L7,L(9)), 

* (L6,L(10)),(L5,L(11)), (L4,L(12)), 

* (L3,L(13)),(L2,L(14)),(L1,L(15)) 
NHOP-2«*IHOP 
IOFST1-IOFST+1 

IF  IFLAG-0  THEN  LOAD  THE  COMPLEX 
EXPONENTIAL  TABLE,  XEXP 

IF(IFLAC.EQ.l)  CO  TO  502 
DARG-PI2/2048.0 
ARC-0.0 

DO  500  J-1,1538 
ARG-ARC+DARC 

XEXP ( J > -CMPLX (COS(ARC) ,SIN(ARG) ) 

500  CONTINUE 
IFLAG-1 

502  IF (NHOF.EQ. NHOPO)  GO  TO  503 
DO  6 J-1,15 
L(J)-NHOP 

6 IF(J.LE.N2)  L(J)-(2**(N2+1-J)) *NHOP 
NTHPOW  - 2**  N2 ; NHOPO-NHOP 
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N4P0V  - N2  /2 
503  NTTL«NTHPOtf#NHOP 

IFCN4P0W .EQ .0)  CO  TO  3 

PERFORM  RADIX  4 TRANSFORM 

DO  2 PASS-1, N4P0W 
NXTLTH-2  **(  N2  -2*PASS> 

LENCTH-4 *NXTLTH 

IDEL-2048/LENGTK 

I ADDH-NXTLTH  *NHQP 

LENCTH-LENCTH *NHOP 

DO  2 J-1,NXTLTH 

IARCl-(J-l) *IDEL 

IARG2-IARG1+IARG1 

IARG3-IARG2+IARC1 

MLOC-IOFSTl~LENCTH+( J-l > *NHOP 

DO  2 SEQLOC-LENCTH,NTTL, LENGTH 

J1  - SEQLOC+MLOC 

J2  - Jl+IADDH 

J3  - J2+IADDH 

J4  - J3+IADDH 

C1«S(J1)+S<J3> 

C2-S( Jl)-S< J3) 

C3-S(J2)+S(J4) 

C-S( J2)-S( J4) 

C4-CMPLX(-R4 ,14) 

S(  JD-C1+C3 

IF(J.EQ.l)  GO  TO  1 

S(J3)«XEXP<IARC1)*<C2+C4) 

S< J2)-XEXP(IARC2) *<C1-C3) 

S ( J4  > -XEXP  < IARC3 ) *(C2-C4> 

GO  TO  2 

1 S(J3)-C2+C4 

S(J2)-C1-C3 
S( J4 )-C2-C4 

2 CONTINUE 

PERFORM  RADIX  2 TRANSFORM  IF  REQUIRED 
C 

3 IF(  N2  .EQ .2  *N4POW)  GO  TO  5 
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NHOP2-NHOP*2 
NSTOP-NTTL+IOFST 
DO  4 J"I0FST1 , NSTOP, NHOP2 
C“S( J)+S< J+NHOP) 

S( J+NHOP) *S(J)-S( J+NHOP) 

S<J)-C 

CONTINUE 

OUTPUT  CURRENTLY  IS  ORGANIZED  WITH 
BIT  REVERSED  ADDRESSING 
THIS  SECTION  PLACES  OUTPUT  IN  THE 
CORRECT  ORDER 

IJ-1 

Jl-1 

DO  7 J2-J1,L2,L1 
DO  7 J3-J2,L3,L2 
DO  7 J4«J3,L4,L3 
DO  7 J5-J4,L5,L4 
DO  7 J6-J5,L6,L3 
DO  7 J7“JS,L7,LB 
DO  7 J8-J7,L8,L7 
DO  7 J9**J8,L9,L8 
DO  7 J10-J9,L10,L9 
DO  7 J11-J10,LU,L10 
DO  7 J12-J11 ,L12 ,L11 
DO  7 J13“J12,L13 , LI 2 
DO  7 J14  -J13,L14 ,L13 
DO  7 JI-J14,L15,L14 
IF<IJ.GE.JI)  GO  TO  7 
KJ-IJ+IOFST 
JK-JI+IOFST 
C-S(KJ) 

S(KJ)-S(JK) 

S(JK)-C 
7 IJ-IJ+NHOP 
J1-NTTL/2+IOFST 
J2-J1+1 

DO  14  J3-IOFSTl,Jl,NHOP 
C-S(J3) 

S(J3)-S(J2)  7Q 


m 
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S(J2)-C 
J2-J2+NH0P 
14  CONTINUE 
RETURN 
END 
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function  of  the  first  kind  and  order  zero  with  an 
imaginary  argument.  The  series  is  shown  below  (4.1). 


EXPND(X2)  - l‘KX2)8+(X2/'2)*+(X2/2)*+(X2/2)8+ . . . 

2 C2 1 > 2 ( 3 1 > 2 (4 1 )2 


(4.1) 
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7.  FORTRAN  LISTING 


FUNCTION  EXPNOCX2) 

XB2SQ-X2 *X2 *0 . 23 
SUM-1. 0+XB2SQ 
ADDON-XB2SQ 
DO  11®  J-2,20 
AJ-FLOAT(J) 

IF ( ABS ( ADDON ) . LT . ABS ( SUM  *1 . 0E-0 G ) ) GO  TO  200 
ADDON-ADDON *XB2SQ/(AJ*AJ) 

110  SUM-SUM-1- ADDON 
20®  EXPND-SUM 
RETURN 
END 
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SUBROUTINE  FFT2D 

1.  PURPOSE 

This  subroutine  computes  the  two-dimensional 
discrete  Fourier  transform  of  a planar  array  of 
samples . 


2.  INPUT  PARAMETERS 


Name 

O/R 

I 

Description 

N2 

R 

I 

The  power  of  2 which  defines 

the  length  of  each  side  of  the 

2-D  transform. 

LRJ 

R 

I 

The  number  of  logical  record 
blocks  in  the  x-directlon  of 
the  input  file  (horizontally). 

LRK 

R 

I 

The  number  of  logical  record 
blocks  in  the  y-direction  of 
the  input  file  (vertically). 

LRJIN 

* 

R 

I 

The  number  of  blocks  to  be 
skipped  in  the  x-direction, 
starting  on  the  left,  before 
storing  the  output. 

LRJVID 

R 

I 

The  width  in  blocks  of  the 
desired  far-field  output. 

LRKIN  R I The  number  of  blocks  to  be 

skipped  in  the  y-direction, 
starting  at  the  top,  before 
storing  the  output. 


LRKWID 


The  height  in  blocks  of  the 
desired  far-field  output. 

85 


PAAS-CPD-FFT2D 
1 MAY  78 


3 .  CALLING  SEQUENCE 

CALL  FFT2D  (N2 , LRJ , LRK , S , LRSDJ , LRSDK , SL> 

Where:  N2  - Power  of  2 that  determines  the  length 

of  each  side  of  the  2-D  transform. 

LRJ  - The  number  of  logical  record  blocks  in 
the  x-direction  of  the  input  file. 


LRK  - The  number  of  logical  record  blocks  in 
the  y-direction  of  the  input  file. 

S - A tvo-dimensio.nal  array  used  to  store 
the  blocks  of  the  input  file  for 
processing . 

LRSDJ  - Object-time  dimension  constant  for  the 
x dimension  of  the  array  S. 

LRSDK  - Object-time  dimension  constant  for  the 
y dimension  of  the  array  S. 

SL  - A one-dimensional  array  equivalenced  to 
S,  used  for  intermediate  2-D 
processing . 


4 .  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 


a.  Reference: 

Hansen,  R.  C.,  Microwave  Scanning  Antennas,  Vol . 2, 
New  York:  Academic  Press,  Inc.,  1964. 

5 .  SUBPROGRAMS  REQUIRED 

CZFFT 
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6.  THEORY  OF  OPERATION 


The  relationship  between  a finite  linear  array  of 
radiators  and  its  corresponding  far-field  is  given  by 
the  following  equation 


E(0) 


N 

- 2 c<n> 

n - 1 


e J n kdxS I n 8 


(A3 .1) 


where  k ■ 2TT/X  and  dx  is  the  interelement  spacing. 


G(n)  is  the  current  gain  of  the  nth  radiator.  This 
equation  assumes  lsotrophic  radiators.  Nov  by  letting 
a new  variable  p be  equal  to  the  following: 


p - NkdxSin0 


the  expression  in  Equation  A3.1  becomes  the  following 
equation 


N 


E(p)  - £ G(n)eJ< 
n - 1 


Pn/N) 


which  is  in  the  form  of  the  IDFT.  This  expression  can 
be  calculated  using  standard  FFT  techniques.  The 
structure  of  the  one  dimensional  FFT  algorithm  requires 
the  input  data  to  be  in  the  order  that  is  shown  in 
Figure  FFT2D-1 . If  this  reorganization  is  not 
implemented  the  output  data  will  have  a 180°  phase 
shift  from  one  point  to  the  next.  This  problem  also 
arises  when  the  2D-FFT  is  performed.  A shuffle  of 
blocks  of  data  rather  than  line  segments  must  be  done 
to  prevent  the  problem  from  occurring.  Figure  FFT2D-2 
illustrates  the  organization  of  the  aperture  shuffle. 
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In  the  case  In  which  the  antenna  aperture  has  a 
smaller  number  of  points  than  the  desired  far-field, 
there  must  be  some  ‘zero'  blocks  or  blocks  loaded  with 
zeros  to  pad  the  Input  aperture  into  the  far-field 
point  configuration.  Figure  FFT2D-3  illustrates  the 
padding  with  the  zero  blocks  and  the  shuffle. 


The  input  aperture  field  with  dimensions  LRJ  x LRK 
blocks  is  split  into  the  four  corners  of  the  transform 
field.  The  transform  field  has  2Na  points  on  a side. 
This  is  illustrated  in  Figure  FFT2D-4 . 


Each  row  of  this  matrix  is  now  transformed,  one  at 
a time,  starting  at  the  top.  Since  the  middle 
<2N2/16)-LRK  rows  are  zero,  the  transform  is  equal  to 
zero.  Therefore,  the  program  skips  these  rows  and 
begins  at  the  top  of  section  B.  This  avoids  a waste  of 
computer  time.  Now  that  the  first  transform  has  been 
executed,  if  the  whole  far-field  is  required,  the 
second  transform,  which  is  the  columns  of  the 
intermediate  result,  must  execute  a complete  2N2  x 2N2 
point  transform.  In  most  cases,  however,  the  whole 
far-field  is  not  required  and  only  a small  vertical 
section  needs  to  be  transformed.  Figure  FFT2D-5  is 
provided  to  illustrate  this  situation.  The  transform 
field  is  2N2/16  blocks  on  a side.  If  the  user  only 
wishes  to  look  at  a section  of  blocks  that  have 
dimensions  LRJWID  x LRKWID,  only  a vertical  stripe 
LRJWID  wide  needs  to  be  transformed.  Since  the 
transform  of  the  other  columns  of  the  matrix  have  no 
effect  on  the  transform  of  the  columns  in  the  stripe, 
it  would  again  be  a waste  of  computer  time.  Only  the 
double  crosshatched  area  of  Figure  FFT2D-3  is  stored  in 
the  output  PRMFL.  If  two  complete  2N2  x 2N2  transforms 
were  executed  the  total  number  of  complex  points 
processed  would  be  2N2*4.  Using  this  scheme  only 
(LRK*16) *(LRJWID*16) *2N2*2  complex  points  are 
processed.  For  the  case  of  N2  - 8,  and  LRK  - LRJWID  ■ 
4,  only  6. 25??  of  the  total  number  of  complex  points  are 
processed . 
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a)  Original  Antenna  b)  Shuffled  Antenna 

Aperture  Aperture 

Figure  FFT2D-3:-  -Organization  of  input  data  with  zero  blocks  required  by 

2-0  FFT  algorithm. 
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b)  Shuffled  Antenna 
Aperture 


Figure  FFT2D-4:  Organization  of  input  data  showing  block  measurements 
for  the  2D-FFT  algorithm. 
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7.  FORTRAN  LISTING 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


100 


220 


230 


239 


240 

290 

300 


SUBROUTINE  FFT2D ( N2 , LR J , LRK , S , LRSD J , LRSDK , S 1 ) 
COMPLEX  ALRZC296) , S(LRSDJ,*  RSDK) , Sl(l) 

CALL  PTIME(OTIME) 

NAMELIST/FILOUT/LRJIN , LRJVID , LRKIN , LRKWID 
READ ( 09 , FILOUT ) 

WRITE ( 06 , FILOUT ) 

LRSIDJ-LRSDJ/16 
ICRNR-LRSID J * < LRKIN+LRKWID ) 

LRSTP-0 
DO  100  1-1,296 
ALRZ(I)-<0. 0,0.0) 

LRNMR-0 

LR1IN-<1+LRK) *LRJ/2 

LRMKR-0 

•LR2-LRSTP+1 

LRNMR1-0 

JST-1 

JSTP-16 

CO  TO  300 

LRST-LRSTP+1 

LRSTP-LRST+ ( LRSID J-LR J > - 1 
DO  240  LR1-LRST,LRSTP 
DO  239  K-1,16 
DO  239  J-JST, JSTP 
S(J,K)-(0. 0,0.0) 

LRMKR-LRMKR+1 

JST-JSTP+1 

JSTP-JST+19 

LR1IN-LRR-LRJ/2 

LRST-LRSTP+1 

LRSTP-LRST+LRJ/2-1 

LRR-LR1IN 

DO  310  LR1-LRST, LRSTP 
LR1IN-LR1IN+1 
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500 


510 


520 


700 

710 


720 

800 


READ(01'LR1IN)<<S<J,K>,J-JST, JSTP) ,K-1 , 1S> 

LRMKR-LRKKR+1 

LRNMR1-LRNMR1+1 

LRNMR-LRNMR+1 

JST-JSTP+1 

JSTP-JST+15 

IF(LRMKR.EQ.LRSIDJ)  CO  TO  400 

IF( (LRSIDJ-LRJ) .NE .0)  CO  TO  230 

CO  TO  250 

IOFST-O 

IHOP-0 

DO  500  K-1,16 

CALL  CZFFT ( S , N2 , IOFST , IHOP ) 
IOFST-IOFST+LRSDJ 
DO  510  JST-1,LRSDJ,1S 
JSTP-JST+15 

WRITE  (03'LR2)((S<J,IO , J- JST , JSTP ) , K-l , 16 ) 

LR2-LR2+1 

CONTINUE 

LRMKR-0 

LR1IN-LR1IN+LRJ 

LRHALF- ( LR J *LRK  )/2 

IF(LRNMRl.EQ.LRHALF)  CO  TO  520 

JST-1 

JSTP-18 

CO  TO  300 

IF ( LRNMR . EQ . < LR J *LRK ) > CO  TO  800 

IFC(LRSIDJ-LRK) .EQ.O)  CO  TO  720 

LRST-LRSTP+1 

LRSTP-LRST+LRSID J- 1 

DO  710  LRl-LRST.LRSTP 

WRITE(03'LR1)  ALRZ 

LRZSTP-LRSID J * ( LRSID J-LRK/2 ) 

IF(LRSTP.CE.LRZSTP)  CO  TO  720 

LRST-LRSTP+1 

CO  TO  700 

LR1IN-LRJ/2 

CO  TO  220 

LRTTL-0 

LR2RL-0 

LR2-LRJIN  93 
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810 


830 


865 
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840 


860 


KST-1 

KSTP-256 

LRMKR-0 

LRHOP-0 

LR22-LR2 

LR2-LR2+1 

READ (03 * LR2+LRH0P)  ( SI  (KLOC)  , KLOC-KST,  KSTP) 

LRMKR-LRMKR+1 

LRHOP-LRHOP+LRSIDJ 

KST-KSTP+1 

KSTP-KST+255 

IF(LRMKR.NE.LRSIDJ)  CO  TO  830 

IOFST-0 

IHOP-4 

DO  840  K-1,16 

CALL  CZFFT ( S , N2 , IOFST , IHOP ) 

IOFST-IOFST+1 

LR2-LR22 

KST-1 

KSTP-256 

LRMKR-0 

LRHOP-0 

LR2-LR2+1 

LRBUMP-0 

LR2RL-LR2RL+1 

IFCLR2+LRHOP  .LT . ( LRKIN  *LRSIDJ ) .OR . 

* (LR2+LRHOP) .CT.ICRNR)  CO  TO  865 
WRITE ( 0 2 ' < LR2RL+LRBUMP ) ) ( S 1 ( KLOC ) , 

* KLOC-KST , KSTP ) 

LRBUMP-LRBUMP+LRJWID 

LRTTL-LRTTL+1 

LRMKR-LRMKR+1 

LRHOP-LRHOP+LRSIDJ 

KST-KSTP+1 

KSTP-KST+2I3 

IF ( LRMKR . NE . LRSID J ) CO  TO  860 
IFCLRTTL  .NE . (LRSIDJ*LRJWID) ) CO  TO  810 
CALL  PTIME(TIME) 

TIME-<TIME-OTIME>  *3600 .0 
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PROCRAM  FFT2DX 


PURPOSE 


The  program  FFT2DX  generates  a far-fleld  complex 
voltage  pattern  from  an  existing  antenna  aperture 
distribution  which  is  stored  in  a PRMFL.  FFT2DX  maps 
the  illumination  to  the  far-field  using  a 
two-dimensional  Fourier  transform.  The  far-field  is 
stored  on  a PRMFL  designated  by  the  user. 


INPUT  PARAMETERS 


Name 


N2 


O/R 

R 


I 

I 


LRJ 


LRK 


LRJIN 


LRJVID 


LRKIN 


Description 


The  power  of  2 which  defines 
the  length  of  each  side  of  the 
2-D  transform. 


The  number  of  logical  record 
blocks  in  the  x-directlon  of 
the  input  file 
(horizontally) . 


The  number  of  logical  record 
blocks  in  the  y-direction  of 
the  input  file  (vertically). 


The  number  of  blocks  to  be 
skipped  in  the  x-directlon, 
starting  on  the  left,  before 
storing  the  output. 


The  width  in  blocks  of  the 
desired  far-field  output. 


The  number  of  blocks  to  be 
skipped  in  the  y-direction, 
starting  at  the  top,  before 
storing  the  output. 
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LRKWID 


The  height  in  blocks  of  the 
desired  far-field  output. 


3. 


RESTRICTIONS,  REQUIREMENTS.  MISCELLANEOUS  DATA 
a.  N2  must  be  in  the  range  4 £ N2  S 1® 


b.  LRJ.LRK  must  be  even  and  exactly  the  same  values  as 
vere  specified  in  loading  the  aperture  distribution 
PRMFL . 


4.  SUBPROGRAMS  REQUIRED 
FFT2D 


3.  THEORY  OF  OPERATION 


FFT2DX  initializes  the  array  dimensions  required 
for  the  subroutine  FFT2D. 
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6 . FORTRAN  LISTING 


PARAMETER  LENG-8192 
COMMON  S(LENG) 

COMPLEX  S,S1<1) 

EQUIVALENCE  <S1<1),S<1>> 

CALL  RANSIZ<01,512) 

CALL  RANSIZ<02,512) 

CALL  RANSIZ(03, 512) 

NAMELIST/FFT/N2 , LRJ . LRK 
READ(05,FFT) 

WRITE ( 06, FFT) 

LRSDJ-2**N2 

LRSDK-LENG/LRSDJ 

IF < LRSDK . GT . LRSD J ) LRSDK-LRSDJ 

CALL  FFT2D ( N2 . LRJ , LRK , S , LRSDJ , LRSDK , SI ) 

CALL  EXIT 

STOP 

END 
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PROGRAM  riumo 


1 . PURPOSE 


This  program  modifies  existing  aperture  current 
distributions  that  are  stored  In  a PRHPL . The  program 
allows  the  user  to  list  and/or  change  individual 
element  values.  The  program  also  allows  the  user  to 
'punch*  holes  in  the  current  distribution  with 
specified  radius  and  center.  The  modified  file  may  be 
either  written  over  the  input  PRMFL  or  may  be  stored  on 
another  user  specified'  PRMFL. 


2.  INPUT  PARAMETERS 


Name 


LRJ 


O/R 

R 


I 

I 


Description 


Number  of  logical  record  blocks 
in  the  x-directlon 
horizontally) . 


LRK 


Number  of  logical  record  blocks 
in  the  y-directlon 
(vertically) . 


IBLK 


A pointer  to  indicate  the 
logical  record  block  to  be 
modified 


JSTRT 


The  horizontal  coordinate  to 
begin  the  element  value  listing. 


KSTRT 


The  vertical  coordinate  to 
begin  the  element  value  listing, 


JSTP 


The  horizontal  coordinate  to 
stop  the  element  value  listing. 
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KSTP 


The  vertical  coordinate  to 
stop  the  element  value  listing. 


IELJ 


The  horizontal  coordinate  of 
the  element  to  be  changed. 


IELK 


The  vertical  coordinate  of 
the  element  to  be  changed. 


VREAL 


The  real  part  of  the  new 
element  value. 


VIMG 


The  imaginary  part  of  the  new 
element  value. 


ICNTJ 


The  horizontal  coordinate  for 
the  center  of  the  hole. 


ICNTK 


The  vertical  coordinate  for  the 
center  of  the  hole. 


XHOLE 


The  radius  of  the  hole  to  be 
punched  in  the  eperture 
Illumination . 


3.  RESTRICTIONS.  REQUIREMENTS,  MISCELLANEOUS  DATA 


a . 


The  values  of  JSTRT , KSTRT , JSTP , KSTP , IELJ , IELK  are 
all  assigned  with  respect  to  the  upper  left  corner 
of  each  block  which  has  the  coordinates  (1,1). 


b.  The  values  of  ICNTJ, ICNTK  are  assigned  with  respect 
to  the  upper  left  corner  of  the  aperture  field 
having  the  coordinates  (1,1). 


4.  SUBPROCRAMS  REQUIRED 


None 
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THEORY  OF  OPERATION 


In  the  section  of  th*  program  that  lists  or  changes 
individual  element  values,  the  program  reads  the 
specified  block  into  an  array  A(J,K)  that  is  complex 
and  has  dimensions  (16,16).  The  index  J increments  the 
fastest.  When  a list  Is  requested,  the  values  of  the 
elements  start  with  the  location  ( JSTRT, KSTRT) , 

( JSTRT+1 , KSTRT) (JSTP , KSTRT) , ( JSTRT, KSTRT+1) , . . ., 

(JSTP, KSTRT+1 ) , (JSTP , KSTP ) , l.e.,  from  the  upper 

left  hand  element  down  to  the  lower  right  hand  element. 
This  process  is  repeated  for  each  block  requested. 

The  element  values  of  each  block  are  changed 
according  to  the  location  and  value  given  by  the 
equation : 


A ( IELJ , IELK ) -CMPLH ( VREAL , VIMG ) 

This  is  repeated  for  the  total  number  of  element 
changes  requested  and  for  each  block  requested. 

Holes  in  the  aperture  illumination  are  punched  with 
a radius  determined  by  KHOLE . The  center  of  the  hole 
is  located  at  the  aperture  coordinate  (ICNTJ,ICNTK) . 

The  values  of  ICNTJ  and  ICNTK  are  assigned  with  respect 
to  the  (1,1)  element  of  logical  record  block  number  one 
having  the  aperture  coordinates  (1,1).  The  upper  left 
element  in  the  aperture  field  has  the  coordinates 
(1,1).  The  program  sequentially  steps  through  the 
aperture  blocks  starting  with  block  one  and  proceeding 
to  block  LRTTL  ("LRJ*LRK).  The  distances  from  the 
elements  in  each  block  to  the  element  located  at 
(ICNTJ, ICNTK)  are  calculated.  A comparison  of  each 
distance  to  the  length  XHOLE  is  made.  If  the  distance 
is  less  than  or  equal  to  KHOLE,  the  element  value  is 
changed  to  CMPLX(0 .0 , 0 .0) . Otherwise,  the  element 
value  is  unchanged.  In  this  way  holes  with  radius 
XHOLE  are  made  in  the  aperture  current  distribution. 
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8.  FORTRAN  LISTING 


COMPLEX  A(16, 16) 

CHARACTER  FILIN  *20 , FILOUT  *20 , Y *1 , X *1 
DATA  IOK/O400000000000/, Y/' Y'/ 

100  FORMAT(V) 

200  WRITE (06, 100)  'INPUT  FILE  NAME* 

READ  100, FILIN 
LUDIN-01 

CALL  ATTACH ( 01 , FILIN ,3,1, ISTAT , ) 
IFCISTAT.EQ.IOK.OR.ISTAT.EQ.0)  CO  TO  400 
WRITE (06, 300) ISTAT 

300  FORMAT ( 'UNSUCCESSFUL  ATTACH  ISTAT- ‘,020) 
CO  TO  200 

400  CALL  RANSIZ(01 ,512) 

WRITE (06 , 100)  'OUTPUT  FILE  NAME* 

READ  100, FILOUT 
LUDOUT-02 

IF ( FILIN .NE. FILOUT)  CO  TO  450 

LUDOUT-01 

CO  TO  500 

450  CALL  ATTACH (02, FILOUT, 3,1, ISTAT,  ) 

IF ( ISTAT. EQ.IOK. OR. ISTAT. EQ.0)  CO  TO  499 
WRITE (06, 300)  ISTAT 
CO  TO  400 

499  CALL  RANSIZ(02, 512) 

500  WRITE (06, 100)  ‘LRJ,LRK * 

READ  100 , LR J , LRX 
LHTTL-LR J *LRK 
LRTTL1-LRTTL+1 

WRITE(06 , 100)  'MODIFY  OR  HOLE?  (0  OR  1)’ 

READ  100, MODFLG 

ITMP-1 

IF ( MODFLC . EQ . 1 ) CO  TO  1000 
600  WRITE(06 , 100)  'IBLK' 

READ  100, IBLK 
IF ( IBLK . EQ . ITMP ) CO  TO  650 
IF(LUDIN.EQ.LUDOUT)  CO  TO  650 
DO  620  IBK-ITMP , IBLK- 1 
READ ( LUDIN ' IBK ) A 
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WRITE  < LUDOUT ' IBK ) A 
CONTINUE 
ITMP-IBLK+1 
READ ( LUDIN ' IBLK ) A 

WRITE(06 , 100)  'ANY  ELEMENTS  LISTED?  CY  OR  N) * 

READ  100, X 

IF(K.NE.Y)  CO  TO  750 

WRITE(06 , 100)  ' JSTRT , KSTRT , JSTP , KSTP ’ 

READ  100 , JSTRT , KSTRT , JSTP , KSTP 

WRITE ( 06 , 7 30 ) < ( A ( J , K ) , J- JSTRT , JSTP ) , K-KSTRT , KSTP ) 
FORMAT ((2E12. 5)/) 

WRITERS,  100)  'ANY  ELEMENTS  CHANCED?  (Y  OR  N)  * 
READ  100, X 
IF(X.NE.Y)  GO  TO  720 

WRITE (06 , 100)  'HOW  MANY  ELEMENTS  CHANCED?* 

READ  100,NELE 
IFCNELE .LE .0)  CO  TO  770 
IF ( NELE . CT . 100 ) CO  TO  750 
DO  760  I«1,NELE 

WRITE (06,100)  * IEL J , IELK , VREAL , VIMC  * 

READ  100, I J, IK, VREAL, VIMC 
A ( I J , IK ) - CMPLX ( VREAL , VIMC ) 

CONTINUE 

WRITE (06 , 100)  'ANY  MORE  MODS  OR  LIST?  (Y  OR  N) * 
READ  100, X 
IF(X.EQ.Y)  CO  TO  700 
WRITE ( LUDOUT ' IBLK ) A 

WRITE (06 , 100)  'ANOTHER  BLOCK?  (Y  OR  N) ’ 

READ  100, X 

IF(X.EQ.Y)  CO  TO  600 

IF ( ITMP . EQ . LRTTL1 ) CO  TO  900 

IF ( LUDIN. EQ. LUDOUT)  CO  TO  900 

DO  800  IBK-ITMP , LRTTL 

READ (LUDIN 'IBK) A 

WRITE (LUDOUT 'IBK) A 

CONTINUE 

WRITE(06 , 100)  'ANY  HOLES?  (Y  OR  N) ' 

READ  100, X 

IF(X.NE.Y)  CO  TO  1700 
LUDIN-02 

IF ( FILIN. EQ.FILOUT)  LUDIN-01 
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1000  VRITEC06 , 100)  'ICNTJ, ICNTK * 

READ  100, ICNTJ, ICNTK 
WRITEC06 , 100)  ’XHOLE * 

READ  100, XHOLE 

Kl-0 

Jl-0 

DO  1500  IBLK-1 , LRTTL 
READ ( LUDIN ' IBLK ) A 
DO  1400  KK-1,1S 
DO  1300  JJ-1,16 
K-K1+KK 
J-J1+JJ 

X JD-ABS ( ICNTJ - J ) *ABS( ICNTJ- J) 

XKD-ABS(ICNTK-K) *ABS(ICNTK-K) 

DST«SQRT(XJD+XKD) 

IF (DST.LE. XHOLE)  A( JJ,KK)-(0 .0,0 .0) 

1300  CONTINUE 
1400  CONTINUE 

WRITE ( LUDOUT ' IBLK) A 
IF(HOD<IBLK,LRJ) .EQ.0)  CO  TO  1430 
Jl-Jl+16 
CO  TO  1300 

1430  K1-(IBLK/LRJ)*16 
Jl-0 

1300  CONTINUE 

WRITE (06 , 100)  'ANOTHER  HOLE?  <Y  OR  N) ' 

READ  100, X 

IF(X.NE.Y)  CO  TO  1600 
LUDIN-02 

IF ( FILIN. EQ.FILOUT)  LUDIN-01 
CO  TO  1000 

1600  WRITE (06 , 100)  'ANY  ELEMENT  CHANCES?  (Y  OR  N) ' 

READ  100, X 

IF(X.NE.Y)  CO  TO  1700 
LUDIN-02 

IF (FILIN .EQ .FILOUT)  LUDIN-01 
CO  TO  600 

1700  WRITE (06 , 100)  'ANOTHER  OUTPUT  GENERATED?  (Y  OR  N) ' 
READ  100, X 

IF(X.NE.Y)  CO  TO  1800 
WRITE(06 , 1730)  ISTAT 
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1750  FORMATC 'DETACH  OUTPUT  FILE  ISTAT-', 020) 
LUDIN-01 
CO  TO  400 

1800  WRITE(06 , 100)  'ANOTHER  FILE  MODIFIED?  (Y  OR  N) * 
READ  100.X 
READ ( LUDIN ' LRTTL ) A 
READ ( LUDOUT ' LRTTL ) A 
WRITE(06 ,1850)  ISTAT 
IF ( LUDOUT. EQ.l)  CO  TO  1899 
WRITE (06, 1750)  ISTAT 

1850  FORMAT < 'DETACH  INPUT  FILE  ISTAT-', 020) 

1899  IF(X.EQ.Y)  CO  TO  200 

1900  CONTINUE 
STOP 
END 


I 
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DOUBLE  PRECISION  FUNCTION  CAM 


1.  MODULE  IDENTIFICATION 


Name 


Classlf ication  Code 


Reference  Number 


Subordinate 


Not  User  Referenced 


2 .  PURPOSE 


This  function  is  used  to  compute  the  value  of  the 
Gamma  function. 


! M 


3.  INPUT  PARAMETERS 


Name  O/R  T 


Description 


R F 


Argument  of  the  Gamma  function 


4 .  CALLING  SEQUENCE 


C - CAM  (X) 

Where:  X is  the  Input  argument 

C contains  the  computed  value  of  the  Gamma 
function . 

3.  RESTRICTIONS,  REQUIREMENTS , MISCELLANEOUS  DATA 

a.  The  argument,  X,  must  be  within  the  following 
range:  0 < X < 1 .0 

b.  The  maximum  error  in  computing  the  Gamma  function 
is  + 3.0  x 10“7 . 

c.  Reference:  Handbook  of  Mathematical  Functions  by 
M.  Abramowltz  and  I.  A.  Slegun,  Dover,  Inc.,  p. 
237. 


> M 


! 1 

i 1. 


I 4 
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d.  External  References: 


None 


e.  Referenced  labeled  common  areas: 


None 


6.  THEORY  OF  OPERATION 


The  polynomial  expansion  for  the  Gamma  function  of 
x+1  is  given  by  the  following  expression: 


(x+1)  - £ bmxB  + e(x) 

n-0 

Where:  bQ  - 1.0 

bt  - 0.577191652 


bo  * 0.988205891 


b3  - -0 .837056937 


b4  - 0.918206857 


b5  - -0.756704078 


b6  - 0 .4^2199394 
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7 . FORTRAN  LISTING 
C 

DOUBLE  PRECISION  FUNCTION  CAM(X) 
DOUBLE  PRECISION  S,GAK 
S-+0.35868343E-1 
S«S*X-0. 183527818 
S«S*X+0. 482199394 
S“S*X-0 .756704078 
S-S*X+0. 91820S857 
s«s*X-0. 897056937 
S«S*X+0. 988205891 
S-S*X-0. 577191652 
CAM«S*X+1.0 
RETURN 
END 
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PROCRAM  PDFESTR 


1.  PURPOSE 


This  program  generates  a histogram  of  the  radiating 
elements  in  a statistically  loaded  aperture.  The  width 
of  each  radius  cell  may  be  varied  and  the  origin  of  the 
radius  is  user  specified. 


2.  INPUT  PARA.IFTSRS 


Marne 

LR.T 


Q/P  T 

f-  I 


LRK 


n 


OFSTJ' 


R F 


OFSTK 


R F 


RING 


R F 


RLIM 


R F 


n-,.  jatawwf.A.i.M«wa«^ 


Description 


Number  of  logical  record  blocks 
in  the  x-direction 
(horizontally) . 


Number  of  logical  record  bloc’rs 
in  the  y-direotion  (vertically) 


Value  added  to  the  calculated 
middle  of  the  horizontal 
aperture  field  length  tr.  give 
the  offset  origin. 


Value  added  to  the  calculated 
middle  of  the  vertical  aperture 
field  length  to  give  the  offset 
origin . 


Incremental  radius  or  radius 
cell  width  used  to  accumulate 
ihe  histogram  values. 


The  maximum  radius  value  of 
interest . 
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^r,;ym 


ICON 


R 


Mode 

- 1 


■ 0 


- 1 


- 2 


flag 

Program  halts 

Histogram  is  normalized  to 
a unit  cumulative 
distribution . 

Histogram  data  is  converted 
to  probability  density 
estimate  data.  The  rav 
data  is  divided  by  the 
product  of  the  cell  width 
and  the  total  number  of 
elements . 

Raw  histogram  data 


NOPACK 


The  number  of  incremental  radius 
histogram  cells,  RINC,  combined 
to  make  each  output  histogram 
cell . 


RESTRICTIONS.  REQUIREMENTS , MISCELLANEOUS  DATA 
None 


4.  SUBPROGRAMS  REQUIRED 


None 


5 . THEORY  OF  OPERATION 


The  program  sequentially  reads  the  logical  record 
blocks  into  the  array  A(J,K),  which  has  dimensions 
(16, IB)  and  is  complex.  The  program  starts  with  block 
number  one  and  continues  through  block  number  LRTTL 
( «LR J *LRK ) . The  elements  of  each  block  are  then 
processed.  If  the  element  value  is  CMPLXC0. 0,0.0)  then 
the  program  proceeds  to  the  next  element.  For  non-zero 
element  values  the  radius  from  the  element  to  the 
origin  is  calculated.  A counter  in  the  appropriate 
radius  cell  of  the  array  XF(IADD)  is  incremented  and 
the  program  proceeds  to  process  the  next  element. 

After  all  of  the  elements  in  the  aperture  field  have 
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been  processed,  the  rav  histogram  is  contained  in  the 
array  XF(1)  with  each  cell  corresponding  to  an  annulus 
with  radius  R - I*RINC  and  width  RINC . The  program 
then  combines  the  radius  cells  in  groups  and  stores 
them  in  the  array  DATOT( J> . NDPACK  consecutive  radius 
cells  are  put  in  each  group,  thus  reducing  the  number 
of  histogram  cells  by  a factor  of  1/NDPACK.  The 
histogram  data  is  then  modified  and  dumped  according  to 
the  value  of  ICON. 


For  ICON-0 


DATOT ( J > -DATOT ( J )/CUM 


where  CUM  is  the  total  number  of  elements.  This 
data  is  normalized  to  a unit  cumulative 
distribution . 


For  ICON-1 


DATOT ( J ) -DATOT  < J )/( CUM  *RINC  *NDPACK ) 

This  data  is  converted  to  probability  density 
estimate  data. 


For  ICON-2 


DATOT ( J ) -DATOT ( J ) 

The  raw  histogram  data  is  outputted. 
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6.  FORTRAN  LISTING 


50 


400 


350 


300 


) 


100 


COMPLEX  A(16,16) 

CHARACTER  FILIN *20 
DIMENSION  XF ( 1000 ) , DATOT ( 300 ) 

DATA  I0K/04 00000000000/ 

WRITE<06,50)  *LRJ,LRK* 

FORMAT(V) 

READ  50,LRJ,LRK 
WRITE <08, 50)  ’OFSTJ,OFSTK’ 

READ  50 , OFSTJ , OFSTK 
WRITE(06 , 50)  ’INPUT  FILE’ 

READ  50, FILIN 

CALL  ATTACH <01, FILIN, 1,1,ISTAT, 

CALL  RANSIZCOl , 512) 

IF<ISTAT  .EQ .IOK .OR .ISTAT.EQ .0)  GO  TO  300 
WRITE <06, 350)  ISTAT 

FORMAT <’  UNSUCCESSFUL* ATTACH  1ST AT- ’,020) 

CO  TO  400 

XMID J-LR J *18/2 .0+0 .5+OFSTJ 
XMIDK-LRK  *1 6/2 .0+0 .5+OFSTK 
WRITE <08, 50)  ’RINC,RLIM* 

READ  50,RINC,RLIM 
NIXF-RLIM/RINC 
DO  100  1-1 ,NIXF 
XF<I)-0 .0 
CONTINUE 
LRTTL-LR J *LRK 
Jl-0 
Kl-0 

DO  200  IBL.K-1  ,LRTTL 
READ<01'IBLK)A 
DO  500  KK-1,16 
K-KK+K1 

DO  600  JJ-1,16 
J-J1+JJ 

IF< A< JJ , KK) .EQ . <0 .0 , 0 .0) ) GO  TO  600 
XJD«<XMIDJ-J>  *<XMIDJ-J) 

XKD-<XMIDK-K) *<XMIDK-K) 

DV-SQRT  < X JD+XKD ) 
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I ADD-IFIX ( DV/RINC ) +1 
IF(IADD.CT.NIXF)  IADD-NIXF 
XF ( I ADD ) “XF ( IADD ) + 1 .0 
600  CONTINUE 
500  CONTINUE 

IF(MOD(IBLK, LRJ) .EQ.0)  CO  TO  550 

Jl-Jl+16 

CO  TO  200 

550  Kl-dBLK/LRJ) *16 
Jl-0 

200  CONTINUE 

700  WRITE(06 , 50)  *ICON,NDPACK  * 

READ  >,  r:.U*«,NCPACK 

IF ( ICON . EQ . - 1 ) CO  TO  1300 

XEND"  •'  '.XF/NDPACK 

JJ— NDPACK 

CUM-0.0 

DO  1000  J-1,KEND 
DEN-0.0 
JJ-JJ+NDPACK 
DO  900  K-l, NDPACK 
DEN-DEN+XF ( K+ J J ) 

900  CONTINUE 

CUM-CUM+DEN 
DATOT(J)-DEN 
1000  CONTINUE 

IF (ICON .EQ .2)  CO  TO  1200 
CUM-1 .0/CUM 

IF(ICON.EQ.l)  CUM-CUM/(RINC*NDFACK> 
DO  1100  K-l.KEND 
DATOT ( K ) -DATOT ( K ) *CUM 
1100  CONTINUE 

1200  WRITEC06 , 1210) (DATOTC J) , J-l , KEND) 
1210  FORMAT (FI 2 .5) 

CO  TO  700 

1300  CALL  DETACH (01, ISTAT,  ) 

CALL  EXIT 

STOP 

END 
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PROGRAM  PLARY 


1 . PURPOSE 


This  program  loads  a PRMFL  with  a user  specified 
antenna  aperture  current  distribution.  The  aperture 
parameters  include  size,  shape,  weighting,  and  several 
deterministic  phase  options,  including  beam  steering. 
Thinned  or  statistically  loaded  apertures  may  also  be 
generated . 


INPUT  PARAMETERS 


Name 


IAPTFL 


O/R  T 
R I 


HEDGE 


0 F 


XHOLE 


0 F 


NMAJOR 


0 I 


NMINOR 


NWIDTH 

NHICH 


0 I 

0 I 


Description 


Determine  the  shape  of  the 
aperture  to  be  loaded. 


- 1 Circular 

« 2 Elliptical 

- 3 Rectangular 


Radius  of  the  circular 
aperture . 


Radius  of  the  hole  in  a 
circular  or  elliptical  aperture. 


Length  of  semi-major  elliptical 
axis  . 


Length  of  semi-minor  elliptical 
axis  . 


Width  of  rectangular  aperture. 
Height  of  rectangular  aperture 
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IWTFLG 


R I 


VTRAD 


0 F 


ZJRAD 


0 F 


ZKRAO 


WTPED 


NWTPOW 


WKASIR 
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Determine  the  weighting 
function  used  in  loading  the 
aperture . 

- 0 Rectangular  weighting 
Cosine  on  a pedestal  to 
power 


- 1 


2 

3 

4 

5 

6 

7 

8 


Blackman 

Kaiser 

Bartlett  or 

Taylor 

Bessel 

Cubic 

Bayl iss 


:riangular 


Radius  of  the  specified 
weighting  function  for  circular 
apertures . 


Half  the  span  of  the  weighting 
function  in  the  x-direction  for 
elliptical  and  rectangular 
apertures . 


Half  the  span  of  the  weighting 
function  in  the  y-direction  for 
elliptical  and  rectangular 
apertures . 


The  height  of  the  pedestal  for 
cosine  on  a pedestal  to  a power 
weighting . 


The  power  of  the  cosine  function 
for  a cosine  on  a pedestal  to  a 
power  weighting. 


The  Kaiser  variable  for  ,the 
trade-off  between  mam  lobe 
width  and  side  lobe  amplitude. 
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BESCAL 


CUBK 


BESEDG 


DB 


NBAR 


DELPHJ 


DELPHK 


PHERX 


PHERY 


NBITS 


3ESERR 
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The  maximum  weighting  amplitude 
at  the  center  of  the  aperture 
for  the  Bessel  weighting. 


The  weighting  amplitude  scaling 
constant  for  the  cubic  weight. 


The  radius  scaling  constant  for 
the  Bessel  weighting. 


The  design  side  lobe  amplitude 
in  dB  for  the  Taylor  or  Bayliss 
weighting . 


The  number  of  zeros  used  to 
approximate  the  Dolph-Chebyschev 
weighting  distribution  in  the 
Taylor  or  Bayliss  weighting. 


Beam  steering  in  degrees  in  the 
x-direction . 


Beam  steering  in  degrees  in  the 
y-directlon . 


Maximum  quadratic  phase  error 
in  degrees  at  the  edge  of  the 
aperture  in  the  x-direction. 


Maximum  quadratic  phase  error 
in  degrees  at  the  edge  of  the 
aperture  in  the  y-direction. 


Number  of  bits  used  to  control 
the  digital  phase  shifters. 


The  maximum  Bessel  phase  error 
in  degrees  at  the  center' of  the 
aperture . 
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BSCAL 


The  radius  scaling  constant  for 
the  Bessel  phase  error. 


LRJ 


The  number  of  logical  record 
blocks  in  the  x-direction. 


LRK 


The  number  of  logical  record 
blocks  1*.  the  y-direction. 


XKK 


The  probability  that  an  element 
is  located  at  the  peak. 


MAD1 


The  starting  address  for 
selecting  random  numbers  from 
the  random  number  array. 

(1  £ KADI  £ 128) 


JRHD 


The  random  number  generator 
initialization  constant. 

(0  < JRND  < 236  - 1) 


3.  RESTRICTIONS.  REQUIREMENTS . MISCELLANEOUS  DATA 


a.  The  values  for  LRJ, LRK  must  be  even. 


b.  Only  circular  arrays  can  be  statistically  loaded 
4.  SUBPROCRAMS  REQUIRED 


EXPND 

BESS 

CAM 

RRAND 

WEIGHT 

FXOPT 
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5.  THEORY  OF  OPERATION 


The  program  PLARY  loads  the  generated  antenna 
aperture  into  a PRMFL.  The  aperture  is  divided  into 
'blocks'  which  are  IS  elements  on  a side  or  a total  of 
256  elements  per  block.  These  blocks  define  the  size 
and  dimension  of  the  logical  records  used  to  write  the 
aperture  into  the  PRMFL.  Thus,  the  number  of  blocks 
needed  to  load  the  aperture  also  corresponds  to  the 
total  number  of  logical  records  required.  The  aperture 
field  or  the  total  number  of  elements  available  is 
defined  by  an  even  number  of  blocks  arranged  in  columns 
(LRJ  columns)  and  an  even  number  of  blocks  in  each 
column  (LRK  rows).  Thus  the  total  number  of  logical 
record  blocks  is  LRJ*LRK. 


The  weights  at  the  edge  of  a circular  aperture  are 
assigned  values  in  a special  way  to  help  smooth  out  the 
granularity  caused  by  approximating  a circular  aperture 
with  a grid  of  rectangularly  spaced  elements.  A square 
with  sides  d/2,  where  d is  the  interelement  spacing,  is 
constructed  centered  at  each  element  of  the  array.  The 
distances  from  the  corners  of  the  square  to  the  center 
of  the  aperture  are  calculated  for  each  element.  The 
number  of  corners  contained  by  the  specified  aperture 
determines  the  weight  of  the  element.  Each  corner 
counts  a weight  of  0.25.  Therefore,  if  all  four 
corners  are  located  within  the  aperture  distribution, 
then  a value  of  1.0  is  assigned.  Similarly,  if  only 
three  of  the  corners  are  within  the  aperture 
distribution  the  value  is  0.75,  etc.  This  also  holds 
for  the  element  values  around  the  hole  of  the  aperture 
distribution.  If  all  four  corners  are  contained  by  the 
aperture  hole,  then  the  element  value  is  0.®.  The 
number  of  corners  contained  by  a radiating  part  of  the 
aperture  determines  the  value  of  each  element.  This 
technique  produces  a smoother  circular  image 
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The  method  used  to  load  the  elliptical  or  the 
rectangular  aperture  is  not  quite  so  sophisticated. 

< The  rectangular  aperture,  since  it  is  being  loaded  in  a 

; rectangular  grid,  is  simply  loaded  based  on  whether  the 

element  is  contained  or  is  outside  of  the  radiating 
portion  of  the  aperture  distribution.  Those  elements 

• contained  inside  the  specified  limits  are  assigned  a 

, value  of  1.0.  All  others  are  given  a value  of  0.0. 

* The  elliptical  aperture  is  loaded  with  the  same 
technique  as  the  rectangular  aperture.  The  value  of 

t each  element  is  determined  by  the  location  of  the 

element  with  respect  to  the  elliptical  edge.  If  the 
element  is  inside  the  ellipse,  the  value  is  1.0, 
otherwise  it  is  0.0.  This  technique  produces  a fairly 
granular  edge.  However,  if  the  ellipse  is  large  this 
effect  is  minimized. 

The  degree  of  phase  accuracy  in  beam  steering  is 
determined  by  the  number  of  controlling  bits,  NBXTS. 
The  value  of  the  least  significant  bit  of  the  beam 
steering  phase  shifter,  XLSB,  can  be  calculated  as 
follows: 

XLSB  - 360.0/2nbits  degrees 

The  beam  steering  phase  shift,  PH,  for  element  K,J  is 
computed  from  the  orthogonal  steering  angles, 

DELPHJ, DELPHIC,  as  follows 


PHI  - DELPHK  * < K-XMIDK ) + UELPHJ*(J-XMIDJ> 


PH  - FLOAT (IFIX(AMQD( PHI , 3G0 .0)/XLSB) ) *XLSB*DTR 


Where:  K,J 


The  element  location  in  the  y and  x 
coordinates  respectively 


XMXDK 


The  middle  of  the  aperture  in  the  y-span 
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XMIDJ  ■ The  middle  of  the  aperture  in  the 
x-span 


DTR  - Degrees  to  radius  conversion  constant 


Civen  the  maximum  quadratic  phase  error  at  the  x 
and  y edge  of  the  aperture,  PHERX,PHERY,  the  element  to 
element  quadratic  phase  error  is  calculated  using 


PKERR  - ( YMUK  * < K-XMIDK ) * *2  +XMU J * ( J -XMIDJ ) * *2 ) *DTR 


Where,  for  a circle 


YMUK  - PHERY/(XEDCE)**2 


KMUJ  - PHERX/(XEDGE>**2 


and  the  other  parameters  have  the  same  meaning  as 
above.  For  the  ellipse,  the  values  for  XMUJ  and  YMUK 
are 


YMUK  - PHERY/<NMIN0R>**2 


XMUJ  - PHERX/(NMAJOR) **2 


For  the  rectangular  aperture  the  values  are 
YMUK  - PHERY/(NHIGH/2) **2 


XMUJ  - PHERX/<NWIDTH/2>**2 


The  Bessel  phase  error  is  determined  by  first 
calculating  the  radius  to  each  element,  then  scaling 
the  radius  by  the  constant  BSCAL.  The  scaled  radius 
XRAD,  is  then  used  as  the  argument  for  evaluation  of 
the  Bessel  function. 


PHBSER  - BESERR*BESS <0.0, XRAD )*DTR 
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Where:  BESS(0 .0,XRAD)  « The  Bessel  function  of  the 

first  kind  and  order  zero, 
evaluated  at  KRAD 


BESERR 


DTR 


A magnitude  scaling  factor, 
determines  the  value  of 
maximum  error  at  the  center  of 
the  aperture 

Degrees  to  radians  conversion 
constant . 


The  total  deterministic  phase  error  at  each  element 
is  the  sum  of  the  three  independent  phase 
contributions . 

PHTTL  - PH  + PHERR  + PHBSER 

Where:  PHTTL  - The  total  phase  error 

PH  - Beam  steering  and  quantization 
phase  error 

PHERR  - Quadratic  phase  error 

PHBSER  ■ Bessel  phase  error. 

The  value  assigned  to  the  element  in  the  array  A(J,K) 
is 

A(J,K)  - A(J,K) *CMPLX( COS (PHTTL), SIN(PHTTL)) 
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FORTRAN  LISTING 


C 

C 28  APR  78 

C 1030 

C 

COMPLEX  A(1S , 16) 

COMMON  IVTFLC , WTPED , NWTPOW , WKASIR . F( 20) , B( 20) , ANG , 
NBAR , BESCAL , CUBK1 , PII 2 , BESS 1 , I A2 , XKX , VMAX , 
BESEDC 

COMMON/BLKRND/HAD1 . JRND , XMEAN , S IC  2 SQ , UL , UEXT 
DIMENSION  U(20) ,Z(20)  , BZEROK20) 

DATA  U/l .2196699 , 2 .2331306 , 3 .2383134 , 4 .2410621, 

& 5. 2427643, 6. 2439216, 7. 2447598, 8 .2493948 , 

& 9 .2498927.10.2462933,11.246624,12  .246900, 

& 13  .247131.14  .247334,19  .247908 , 1-  .247663, 

& 17  .247796,18.247920,19.248027,20.248129/ 

DATA  BZERO1/0 . 986067 , 1 .6970909 , 2 .7171939 , 3.726137 , 
& 4 .7312271.9.7349209,6.7368281,7.7389396, 

& 8.7398909,9.7408949,10.7417439.11,7424479, 

8.  12.7430408,13.7439477,14.7439896,19.7443679, 

& 16.7447044.17.749003,18.7492697,19.7499093/ 

DATA  IOK/O400000000000/ 

CHARACTER  OUTFIL *20, X*1 , Y*1 ,N«1 
DATA  Y,N/'Y' , N'/ 

CALL  FXOPT< 68, 1,1,0) 

PI-3 .1419926 
PII2-2 .0/<PI*PI) 

BESS1-1 .0/BESSC0. 0,0.0) 

990  WRITE(06 , 140)  'OUTPUT  FILE  NAME* 

READ  140 , OUTFIL 

CALL  ATTACH (01, OUTFIL, 3,1, 1ST AT,  ) 

IFdSTAT .EQ .IOK .OR .ISTAT .EQ .0)  CO  TO  141 
WRITE(06 , 149)  ISTAT 

149  FORMAT ( 'ATTACH  FAILED  ISTAT  - ',020) 

CO  TO  990 

141  CALL  RANSIZC01 , 912 ) 

NAMELIST/APETUR/I APTFL , HEDGE , XHOLE , NMA JOR , 

& NMINOR , NWIDTH , NHIGH 

WRITE (06 ,140)  'STATISTICAL  TAPER?' 
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610 


620 


630 

640 

660 


670 


810 


830 
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800 

710 


720 


READ  140, X 
IAZ-1 

IF(X.EQ.N)  CO  TO  600 
IAZ-0 

WRITE(06, 140)  'XKK, KADI, JRND' 

READ  140 , XKK , MAD1 , JRND 
WRITE (06, 140)  ’IAPTFL  * 

READ  140, IAPTFL 
CO  TO(610, 620, 630), IAPTFL 
WRITE (06, 140)  'XEDGE, XHOLE  * 

READ  1 40, XEDCE, XHOLE 
CO  TO  640 

WRITE (06,140)  ' NM A JOR , NMINOR , XHOLE ' 

READ  140, NMAJOR, NMINOR, XHOLE 
CO  TO  640 

WRITE(06, 140)  ’NWIDTH,NHICH* 

READ  140, NWIDTH , NHICH 

CONTINUE 

CO  TO  670 

WRITE<06 , 140)  'INVALID  IWTFLC' 

NAMELIST/W AIT/I WTFLC , WTRAD , Z JRAD , ZKRAD , VTPED , 
NWTPOW , WKASIR , BESCAL , CUBK , BESEDC 
WRITE (06, 140)  'IWTFLC* 

READ  140, IWTFLC 

IF ( IWTFLC . EQ . 8 . AND . IAPTFL . NE . 1 ) CO  TO  660 
IF ( IWTFLC . EQ . 7 . AND . IAPTFL . EQ . 1 ) CO  TO  660 
IF (IWTFLC .EQ .0)  CO  TO  700 
CO  TO  (810, 830, 830), IAPTFL 
WRITE (06, 140)  'WTRAD' 

READ  140, WTRAD 
CO  TO  800 

WRITE(06 , 140)  'ZJRAD, ZKRAD’ 

READ  140, ZJRAD, ZKRAD 
WTRAD-AKAX1 (ZJRAD, ZKRAD) 

CO  TO  (710, 700, 730, 700, 720, 740, 750, 720), IWTFLC 
WRITE(06 , 140)  'VTPED, NWTFDW ' 
ilt  AD  140,  WTPED,  NWTPOW 
CO  TO  700 

WRITE(06 , 140)  *DB,NBAR ’ 

READ  140,DB,NBAR 
NAMELIST/TAYL/DB , NBAR , SIC 
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RAT-10. 0**(DB/2®.0> 

AA-ALOC  < RAT+SQRT < RAT  *R  AT- 1 ) ) /PI 
AASQ-AA  *AA 

IF(IWTFLG .EQ.8)  CO  TO  760 

SIC-U ( NBAR  >/SQRT ( A ASQ+ (NBAR-0 . 3 ) * *2 ) 

SICSQ-SIC *SIG 
DO  232  I-1,NBAR-1 
FNUH-1 .0 
FDNM-1.0 
T-U(I) *U(I) 

XII— 0.5 

DO  234  II-1,NBAR-1 
XII-XII+1 .0 

FNUM-FNUM*(1 .0-T/(SIGSQ*( AASQ+(XII*XII> ) ) ) 
IF(II.EQ.I)  CO  TO  234 
FDNM-FDNM * ( 1 . 0-T/ (U(II)*U(II))) 

234  CONTINUE 

ARC-PI *U( I) 

F ( I ) —BESS  ( 0 . © . ARC  > *FNUK/FDNM 
F(I)-F<I)/(BESS(0 .0, ARC) **2) 

232  CONTINUE 
CO  TO  700 

730  WRITE(®6,140)  ’WKASIR* 

READ  140, WKASIR 
CO  TO  700 

740  WRITE (06 ,140)  ‘BESCAL , BESEDC ' 

READ  14®, BESCAL, BESEDC 
CO  TO  700 

730  WRITE (06 ,140)  *CUBK* 

READ  14®, CUBK 

XX-SQRT ( WTR AD  *WTR AD/3 .0) 

CUBK1-CUBK/ABS ( XX  * ( XX-WTRAD  > * ( XX+WTRAD ) ) 

CO  TO  700 

760  Z(l)-0 .9858302+0 .0333885  *DB+0 .00014064  *DB*DB 

Si  -0  . 0000019  *DB  *DB  *DB+0 . 00000001  *DB  *DB  *DB  *DB 

Z(1)-Z(1)*Z(1> 

Z( 2 >-2.00337487+0. 011 41548 *DB+0 . 0004159 *DB *DB 
Si  -0  .00000373  *DB*DB*DB+0 .00000001  *DB *DB*DB*DB 

Z(2)-Z(2)*Z(2) 

Z(3)-3 .00636321+0. 00683394 *DB+0. 00029281 *DB *DB 
i -0.000001 6 1*DB*DB*DB 
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Z(3)-2(3)*Z(3) 

Z(4)»4 .00518423+0 .00501795  *DB+0 .00021735  *DB*DB 
-0 . 00000088  *DB *DB *DB 
Z(4)-2(4)*Z(4> 

IF(NBAR.LE.4)  CO  TO  762 
DO  761  I-5,NBAR 
Z<I)-AASQ+I*I 
CONTINUE 

SIG-BZER01 (NBAR+l )/SQRT(Z (NBAR) ) 

SIGSQ-SIG *SIG 
DO  765  I«1,NBAR 
FNUM-1 .0 

T-BZER01 ( I ) *BZER01 < I ) 

FDNM- 1 . 0-T/( BZER01 ( 1 ) *BZER01 (1 ) ) 

IF(I.EQ.l)  FDNH-1.0 

DO  770  II-1,NBAR-1 

FNUM-FNUM  * < 1 . 0-T/<  SICSQ  *Z ( II ) ) ) 

IF(I-l.EQ.II)  CO  TO  770 

FDNM-FDNM  * ( 1 . 0-T/CBZERO1 C II+l ) "BZER01 ( II+l ) ) ) 
CONTINUE 

B<I)-<2 .©*T/BESS<1 .0,BZERO1(I) *PI> ) *FNUM/FDNM 
CONTINUE 

P0-0 .4797212+0 . 01 453692 *(DB)-0 .00018739*(DB*DB> 

+0 . 000002 1 8 * ( DB  *DB *DB ) -0 . 00000001  * ( DB  *DB *DB *DB ) 
P0-P0*SIG 
P0SQ-P0*P0 
PIP0SQ-P0SQ  *PI *PI 
FNUM-1 .0 

FDNM-1 . 0-P0SQ/CBZERO1 ( 1 ) *BZER01 (1 ) ) 

DO  772  1-1 ,NBAR-1 

FNUM-FNUM  *< 1 .0-P0SQ/C  SICSQ *Z(I) ) ) 

FDNM-FDNM  * ( 1 . 0-P0SQ/<  BZER01 ( 1+ 1 ) *BZER01 ( 1+1 ) ) ) 

CONTINUE 

C-PIP0SQ-1 .0 

C«C*BESS<1 .0,P0*PI) *FNUM/FDNM 
C-l .0/C 

DO  773  I-l.NBAR 
B(I)—  B(I)*C 
CONTINUE 

WRITE(06 , 140)  ‘NBITS  ' 

NAMELIST/PHASE/DELPHJ , DELPHK , PHERX , PHERY , 
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& NBITS , BESERR , BSCAL 

READ  140, NBITS 

WRITE(06 , 140)  ‘ANY  BEAM  STEERING?* 

READ  140, X 
IF(X.EQ.N)  GO  TO  900 
WRITE(0S , 140)  'DELPH J , DELPHK  * 

READ  140, DELPH J, DELPHK 

900  WRITE(06 , 140)  'QUADRATIC  ERROR?* 

READ  140, X 

IF(X.EQ.N)  GO  TO  910 

WRITE (06, 140)  *PHERX,PHERY* 

READ  140 , PHERX , PHERY 

910  WRITE (06 ,140)  ‘BESSEL  ERROR?* 

READ  140, X 
IF(X.EQ.N)  GO  TO  134 
WRITE(06 , 140)  ‘BESERR, BSCAL* 

READ  140, BESERR, BSCAL 

134  WRITE(06 , 135) 

135  FORMAT ( * LRJ , LRK ‘ ) 

NAMELIST/BLOCK/LRJ , LRK 
READ  140, LRJ, LRK 

140  FORMAT(V) 

IF(MOD(LRJ, 2) .NE.0.OR  .MOD(LRK,2) .NE.0)  GO  TO  134 
WRITE ( 6 , APETUR ) 

WRITE (6, WAIT) 

WRITE ( 6 ,TAYL) 

WRITE (6, PHASE) 

WRITE (6, BLOCK) 

XEDCE2-XEDGE  *XEDGE 
XHOLE2 -XHOLE  *XHOLE 
NSIDEJ«LRJ*16 
NSIDEK“LRK*16 
NCENTJ-NSIDEJ/2+1 
NCENTK-NSIDEK/2+1 
NCNT1J-NSIDEJ/2 
NCNT1K-NSIDEK/2 
XMIDK-NSIDEK/2+0.5 
XMIDJ-NSIDEJ/2+0.5 
LRTOTL-LR J *LRK 
DTR-0. 017453 
XLSB-360 .0/2 .0**NBITS 
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DXLSB-1 . 0/XLSB 
XLSB-XLSB*DTR 
CO  TO  (10,20,30), IAPTFL 
10  YMUK-PHERY/XEDCE2 
XMU J-PHERX/XEDGE2 
CO  TO  40 

20  YMUK-PHERY/(FLOAT(NMINOR)  > **2 
XMU J«PHERX/< FLOAT (NMAJOR) ) **2 
CO  TO  40 

30  YMUK-PHERY/< FLOAT (NHICH)/2 .0) **2 

XMU J -PHERX/ ( FLOAT ( N WIDTH ) /2 .0) **2 
40  CONTINUE 
Kl-0 
Jl-0 

IF(IAZ.EQ.l)  CO  TO  680 

WMAX-0.0 

IAZ-1 

DO  650  IRAD- 1 , IFIX ( WTRAD+ 1 ) 

RAD-FLOAT (IRAD-l) 

CALL  WEICHT ( RAD , WTR  AD , WFUNC  > 

WMAX-AMAX1 ( VM  AX , VFUNC ) 

650  CONTINUE 
IAZ-0 

680  CONTINUE 

NAMELIST/ST AT/XKK , WHAX 
WRITE ( 06. STAT) 

DO  510  LRl-l.LRTOTL 
DO  50  KK-1,16 
K-K1+KK 
DO  50  JJ-1,16 
J-J1+JJ 

PH-DELPHK  *<  K-XMIDK ) +DELPHJ * ( J-XMIDJ ) 

PH-FLOAT ( IFIX ( AMOD ( PH , 3 6 0 . 0 ) *DXLSB ) ) *XLSB 
PHERR«(YMUK*< K-XMIDK) **2+XMUJ*< J-XMIDJ) **2) *DTR 
XJ- J-XMIDJ 
XK-K-XMIDK 

XRAD-SQRT (XJ  *XJ+XK  *XK ) *BSCAL 
PHBSER-BESERR  *BESS  < 0.0, XRAD ) *DTR 
PH-PH+PHERR+PHBSER 
A ( J J , KK ) - CMPLX ( COS ( PH ) , SIN ( PH ) ) 

50  CONTINUE 
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CO  TO  (100, 300, 400) , IAPTFL 
100  DO  210  KK-1,16 
K-K1+KK 

XKSQ-(K-XMIDK)**2 
DO  200  JJ-1,16 
J-J1+JJ 

ANC-ATAN2 ( ( K-XMIDK > , < J-XMID J ) ) 

XLSQ-XKSQ+( J-XMID J)  **2 
DIST-SQRT(XLSQ) 

IF(IAZ .EQ .0)  CO  TO  205 
CNRl-( J-NCENTJ) **2+(K-NCENTK) **2 
CNR2-(J-NCNT1J) **2+(K~NCENTK> **2 
CNR3-( J-NCENTJ) **2+(K-NCNTlK) **2 
CNR4-( J-NCNT1J) **2+(K-NCNTlK) **2 
ICNFL-0 

IF ( CNR1 . CT . XEDCE2 . OR . CNR1 . LT . XHOLE2 ) ICNFL- 1 
IFCCNR2 .GT.XEDCE2 .OR .CKR2 .LT .XHOLE2)  ICNFL-ICNFL+1 
IF(CNR3 .CT.XEDCE2 .OR.CNR3 .LT.XHOLE2)  ICNFL-ICNFL+1 
IF (CNR 4 .CT.XEDCE2 .OR.CNR4 .LT.XHOLE2)  ICNFL-ICNFL+1 
IF(ICNFL .EQ .0)  CO  TO  205 
CO  TO( 201 ,202, 203, 204), ICNFL 

201  A(JJ,KK)«A(JJ,KK) *0.75 
CO  TO  205 

202  A( JJ,KK)-A( JJ,KK) *0 .5 
CO  TO  205 

203  A( JJ,KK)-A( JJ,KK) *0 .25 
CO  TO  205 

204  A(JJ,KK)-(0. 0,0.0) 

205  IF ( IWTFLC . EQ . 0 ) CO  TO  200 
CALL  WEIGHT(DIST , WTRAD , WFUNC ) 

IF ( IWTFLC. NE. 8)  CO  TO  211 
A(JJ,KK)-A(JJ,KK) *CMPLX(0 .0, WFUNC) 

CO  TO  200 

211  A(JJ, KK)-A( JJ,KK) *WFUNC 
200  CONTINUE 
210  CONTINUE 
CO  TO  500 

300  XMA JOR-NMA JQR  *NMA JOR 
XMINOR-NMINOR *NMINOR 
DO  310  KK-1,15 
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310 


315 

317 


330 

320 


400 


410 


YKSQ-< K-XMIDK)* *2 
DO  310  JJ-1,16 
J-J1+JJ 

XKSQ*( J-XMIDJ) **2 
XLSQ-XKSQ+YKSQ 

ELPSQ-YKSQ/XMINOR+XKSQ/XMAJOR 
IF ( ELPSQ . CT . 1 . 0 . OR . XLS Q . LT . XH0LE2 ) 

A(JJ,KK)-<0. 0,0.0) 

CONTINUE 

IF < IWTFLC .EQ.O)  CO  TO  500 
SCN-1 .0 
DO  320  KK-1,16 
K-K1+KK 

XKPT-ABS ( K-XMIDK ) 

IF ( IWTFLC . NE . 7 ) CO  TO  315 
WFNK-1 .0 
CO  TO  317 

CALL  WEIGHT ( XKPT , ZKR AD , WFNK > 

CONTINUE 

DO  330  JJ-1,16 

J-J1+JJ 

XJPT-ABS< J-XMIDJ) 

CALL  WEIGHT  < X JPT , Z JR AD , WFN  J ) 

IF  (IWTFLC  .EQ  .7  .AND . J .LT  .XMIDJ)  SCN— 1 .0 

A ( J J , KK ) - A ( J J , KK ) *WFNK *WFNJ *SCN 

CONTINUE 

CONTINUE 

CO  TO  500 

XWIDTH-NWIDTH/2 

XHICH-NHICH/2 

DO  410  KK-1,16 

K-K1+KK 

YK-ABS( K-XMIDK) 

DO  410  JJ-1,16 
J— Jl+JJ 

XK- ABS ( J -XMIDJ ) 

IF(YK .CT .XHIGH .OR .XK .CT .XWIDTH)  A( JJ, KK)- (0 .0,0 .0) 
CONTINUE 

IF ( IWTFLC . EQ . 0 ) CO  TO  500 
SCN-1 .0 
DO  420  KK-1,16 
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1, 


K-K1+KK 

XKPT-ABS(K-XMIDK) 

IF ( I VTFLC . NE . 7 ) CO  TO  413 

WFNK-1.0 

CO  TO  417 

413  CALL  WEIGHT(XKPT , ZKRAD , WFNK ) 

417  CONTINUE 

DO  430  JJ-1,16 
J— Jl+JJ 

X JPT - ABS  < J -XMID J > 

CALL  WEIGHT < X JPT, Z JR AD, VFNJ> 

IF(IWTFLC  .EQ  .7  .AND . J .LT  .XMIDJ)  SCN— 1 .0 
A ( JJ , KK ) - A ( J J , KK ) *WFNK *WFNJ *SCN 
430  CONTINUE 
420  CONTINUE 
300  WRITE < 01 'LR1>  A 

IF(MOD(LRl ,LRJ) .EQ.0)  CO  TO  503 

Jl-Jl+16 

CO  TO  510 

505  K1“(LR1/LRJ)  *1S 
Jl-0 

510  CONTINUE 

CALL  DETACH ( 01, ISTAT,  ) 

WRITE(06, 146)  ISTAT 
146  FORMAT < 'DETACH  ISTAT-', 020) 

CALL  EXIT 

STOP 

END 
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SUBROUTINE  PLOTD 


1 . PURPOSE 


The  purpose  of  this  subroutine  is  to  transmit  a 
data  array  to  the  DUIS  for  subsequent  plotting.  The 
data  samples  to  be  plotted  must  be  equally  spaced,, 
i.e.,  the  independent  variable  increment  between 
samples  must  be  a constant. 


2.  INPUT  PARAMETERS 


Name 


O/R  T 


Description 


DV 


R F 


The  array  containing  the 
dependent  variable  values  to 
be  plotted 


NOUT 


R I 


The  number  of  samples  to  be 
processed  in  the  array,  DV 


OR 


R F 


The  origin  of 'the  independent 
variable,  i.e.  the  value  of  the 
independent  variable  that 
corresponds  to  the  sample,  DV(1) 


DEL 


R F 


The  independent  variable 
increment  between  samples 


TH 


0 F 


Variable  used  to  accumulate  the 
maximum  dependent  variable 
modulus . 


3.  CALLING  SEQUENCE 

CALL  PLOTD (DV, NOUT, OR, DEL, TH) 
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RESTRICTIONS,  REQUIREMENTS,  MISCELLANEOUS  DATA 


a.  This  subroutine  is  structured  to  process  only  data 
which  has  a fixed  independent  variable  increment 
between  samples. 


b.  The  output  data  from  this  subroutine  has  the 
standard  DUXS  plot  data  format. 


The  output  data  from  this  subroutine  .has  an 
accuracy  of  12  bits. 


The  variable  TH  can  be  used  to  accumulate  the 
maximum  modulus  value  over  a number  of  plots.  This 
is  mainly  used  in  processing  three-dimensional 
plots . 


SUBPROCRAMS  REQUIRED 


ALOC10 

CHOP 


G.  THEORY  OF  OPERATION 


The  first  operation  performed  by  this  subroutine  is 
to  scan  the  input  array,  DV,  to  determine  the  largest, 
XMAX,  and  smallest,  XMIN,  dependent  variable  values. 

The  parameter  TH  is  then  updated. 


The  dependent  variable  range  XMAX**XMAX-XMIN  is 
computed  from  the  scanner  output  parameters.  The  mean 
of  the  dependent  variable  array,  BIAS,  is  then 
computed.  Next,  two  integers,  J anc  K,  are  found  such 
that  the  following  condition  is  satisfied: 


<J-1)*10K  < XMAX  < JU0K 


The  LSB  to  be  used  in  digitizing  the  data  is  then 
computed  as  follows: 


XLSB-FLOAK J*10k)/4000 .0 


.A,  E.  irvib  *1  M re£Z  ft/iX  lAtMt*1. 
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Nov  all  of  th8  parameters  necessary  to  characterize 
the  plot  data  have  been  determined.  Therefore,  they 
are  transmitted  to  the  DUIS  in  the  plot  header>*in©,cord . 
This  record  contains  the  following  parameters: 

NOUT , OR , DEL , BIAS , KLSB  ,3 , r-j , rjj  , y j > ; » }•./.■•;  i -A . 7 Win -'Hi 

*'*«;.  lAf/iviTiH1'.;. 

Each  plot  point  is  converted,  to  an  intege£finumb,£r 
in  the  following  manner:  • >■  , v<; -\a ;./ih 


PTL  - (DV<J)nBIAS)y<XLSB  Ci  W-U  ,V 

IDAT  - PTL+SICN(®.5,PTL)„.;<,v 


’,‘S 


« c vt a 


This  integer  number  is  converted  to  tuorjASGII  _e  G*. 
characters  and  packed  into  a plo£  idat§/jpecond^byvth'fl>. 
subroutine  CHOP.  Each  time  a record Is :^ill;ed?.C3.3„.splot 
points)  it  is  transmitted  to  the  DUIS.  l&i 

When  all  of  the  plot  points  have  been  tr$nsmljtt;erd 
to  the  DUIS  a plot  terminator  record  is  transmitted  to 
signify  that  the  pi ot ;ri'S < qompl e,te . ? yGpntrgl  > is ^e  tuijned 
to  the  calling  (sub)program . 

cvc,  HriHAJirau  oi 

V.  k HP  -Kk  tiX  M s\¥>W 

>- -GO'  (c.  ms  ■ "-.cv-.n y.vz 

i\  I ?,  /; 4 V M )( 4 a A fi 1 S 
nGAs':-  jo.-! a ■». rrci 

\ jt  •«;  Tf-o  ?j.t  r 


■"'/  r 1 ■ O '•  a . Z " < ,il  ci ) ’■>  I 

1 


D 

«v 

U 

■J 

0 

3 


*;js 


A<1  r , ’ . Vf  . r • -V'. 

; v .? ..  Oj  nc 

\ ■■  >*  04  x 

• $ T,  ' • 1 ; »•>’  l H L r <5  3 X 

f ..  i > ...  r , < .n  ' '".i 

t <'-  V , ’ '**  ' . '*  t ,\»s  { ' ; A A*  5?  S*  r 
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7 . FORTRAN  LISTING 


C 

C 

C 

C 


C 

C 

C 

C 


10® 


C 

C 

C 


C 

C 

C 


123 


17  APR  78 
142® 


SUBROUTINE  PLOTD ( DV , NOUT , OR , DEL , TH ) 
DIMENSION  DV(1) , ILINEC35) 

XMAX-DV(l) 

XMIN-DV(l) 


SCAN  DATA  ARRAY  TO  FIND  LARGEST  (XMAX) 
AND  SMALLEST  (XMIN)  VALUES 


DO  1®0  J-l /NOUT 
XMAX- A MAXI (XMAX, DV( J) > 
XMIN-AMIN1 (XMIN , DV( J) ) 
CONTINUE 


UPDATE  TH 


TH-AMAX1 < ABS ( XMAX >, ABS ( XMIN ) , TH ) 


SECTION  TO  DETERMINE  BIAS  AND  XLSB 


14® 


150 

160 


124 

125 


XMAX-XMAX-XHIN 

XF(XMAX .LT .1 .0E-10)  XMAX-1.0E-10 
BIAS-0 . 5 *XMAX+XMIN 
PTL-ALOC 1® ( XMAX ) 

ITEST-IFIX(PTL) 

IF(PTL)  150,140,14® 

X-XMAX/(10 .0**ITEST> 

CO  TO  160 

X-XMAX  * < 10 . 0 * *< I ABS  C ITEST ) +1 > ) 
N-IFIX<X+0 .98) 

IF(PTL)  123,124,124 

XMAX-FLOAT ( N >/( 10 .0**(IABS(ITEST)+1) ) 

CO  TO  125 

XMAX-FLOAT(N) *<10 .0**ITEST) 

XLSB-XM AX/4 000 .0 

140 


SSrf 


itiaiataaflacittfcfltftA* 


iSJX 


k — u V^.'> *•  <d 


w 


w 

% 


m 


i 


m 

&■■■ 


¥ 


I 


I 


m 


p. 


i 
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TRANSMIT  THE  PLOT  HEADER  RECORD  TO  DUIS 


WRITE (06 , 1000)  NOUT, OR , DEL , BIAS , KLSB 
1000  FORMAT(4HzzX  .16, 4< ' , MPE12 .5)) 


SECTION  TO  CONVERT  DATA  TO  12  BIT  FORM 
AND  TRANSMIT  TO  DUIS 


181 

200 


201 

222 

1002 


220 

180 


IFLG-0 

J-l 

K-l 

IF(J-NOUT)  200,200.201 
PTL- ( D V ( J ) -BIAS ) /XLSB 
IDAT-PTL+SICN ( 0 . S , PTL ) 

CALL  CHOP ( IDAT , ILINE ( X ) ) 
IF(K-33)  220,222,222 
IFLC-1 

WRITE(06  1002) (ILINE(L) ,L«1,33) 
FORMAT( ’>z ’ , 33R2 , *>z*) 

K-l 

IF(IFLC)  180,180,500 

K-K+l 

J-J+l 

CO  TO  181 


C 

C 

C 

C 


PLOT  TRANSMISSION  IS  COMPLETE. 
TRANSMIT  THE  TERMINATOR  RECORD 


f 

500 

WRITE <06 

1 

1003 

FORMAT ( * 

I 

RETURN 

1 

1 

& 

END 
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PROGRAM  PLTDVR 


1.  PURPOSE 

The  purpose  of  this  subroutine  is  to  process  data 
for  transmission  to  the  DUIS  for  subsequent  use  in 
preparing  3-dlmensional  plots.  The  types  of  data 
normally  processed  by  this  subroutine  are  antenna  far 
field  patterns  and  antenna  aperture  illumination 
distributions . 

2.  INPUT  PARAMETERS 


Name 


O/R  T 


Description 


AFILE  R C 


JVID 


R I 


JWIDSP 


The  name  of  the  TSS  PRMFL  which 
the  user  wishes  to  process 

The  data  array  width  in  blocks. 
This  parameter  corresponds  to 
LRJWID  in  the  program  FFT2DX  or 
LRJ  in  the  programs  PLARY, 
RNDERR , or  FILMOD. 

The  width  in  blocks  of  ths 
vertical  strip  to  be  processed 


NBMAX 


The  number  of  the  last  block  to 
be  processed 


LREC 


ISGN 


The  number  of  the  first  record 
to  be  processed 

Data  processing  mode  flag. 

- 0 modulus  data  is  processed 
■ 1 real  component  only  is 
processed . 
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3.  RESTRICTIONS.  REQUIREMENTS . MISCELLANEOUS  DATA 

a.  The  maximum  value  of  the  parameter  JSID  is  10. 

This  corresponds  to  a row  length  of  160  samples. 
This  limitation  is  arbitr  *»y  and  was  chosen  to 
minimize  the  memory  required  by  this  program. 

b.  It  should  be  recognized  that  the  run  time  of  this 
program  is  directly  proportioned  to  the  number  of 
plot  points  transmitted  to  the  DUIS.  The  number  of 
plot  points  is  given  by  the  following  expression: 

NPTS* JSID  *1 6 *NROVS 

The  time  required  to  transmit  each  plot  point  is 
16.67  ms  at  1200  Baud  and  66.67  ms  at  300  Baud. 

For  example,  a 3-D  display  4 blocks  wide  and  64 
rows  tall  would  require  68.27  seconds  at  1200  Baud. 

4*16*64*16 .67ms  - 68.27  seconds 

At  300  Baud  the  time  required  is  4.99  minutes. 

c.  This  program  is  designed  to  work  only  with  the 
DUIS . 

4.  SUBPROGRAMS  REQUIRED 

ATTACH 

RANSIZ 

PLOTD 

DETACH 

EXIT 

5.  THEORY  OF  OPERATION 


After  startup  a message  is  printed  to  the  qser 
requesting  the  name  of  the  TSS  PRMFL  that  is  to  be 
processed.  The  user  reply  is  tested  to  see  if  it  is 
the  word  ’STOP’.  If  so,  program  execution  is 
terminated.  Otherwise,  an  attempt  is  made  to  access 
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the  file  via  the  system  subroutine,  ATTACH.  If  the 
file  cannot  be  accessed  the  system  6rror  code  is 
printed  and  the  user  is  requested  to  try  again. 

Providing  an  input  file  has  been  successfully 
attached,  the  following  parameters  are  requested  from 
the  user: 


JWXD , JWIDSP , NBMAX 

The  parameter  JWIDSP  is  tested  to  ensure  that  it  does 
not  exceed  the  smaller  of  JWID  or  10.  The  following 
parameters  are  then  requested  from  the  user: 

LRECF , ISCN 


The  parameter  LRECF  is  tested  to  see  if  it  is  zero  or 
negative.  If  so,  control  is  transferred  to  Statement  # 
400  and  the  file  is  deaccessed  via  a call  to  the  system 
subroutine  DETACH.  The  program  then  returns  to  the 
procedure  for  accessing  another  file  (Statement  #110). 


The  JWIDSP  blocks,  starting  with  record  number 
LREC,  are  loaded  into  the  array  A(J,K).  Next,  the  16 
rows  of  data  are  processed  one  row  at  a time  for  output 
to  the  DUIS  via  the  subroutine  PLDTD.  If  ISCN-1  then 
the  real  components  of  the  data  samples  are  placed  in 
the  array  DV.  If  ISGN-0  then  a bipolar  modulus 
function  is  computed  from  the  output  data.  The  reason 
for  computing  a bipolar  modulus  is  that  in  the  DUIS  the 
data  to  be  plotted  is  normally  interpolated  to  provide 
smoother  curves.  The  interpolation  functions  cannot 
accurately  process  data  that  possesses  discontinuities 
which  is  the  case  for  a true  modulus  function.  The 
procedure  used  to  convert  the  modulus  is  as  follows.  It 
has  been  noted  that  the  modulus  can  become  zero  only  if 
the  real(x)  and  imaginary(y)  components  are  both  zero 
simultaneously.  Also,  it  is  known  that  both  x and  y are 
continuous  functions.  Therefore,  if  both  x and  y 
reverse  signs  between  two  sample  points  then  the 
modulus  must  have  gone  through  zero.  Therefore,  the 
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scheme  used  was  to  reverse  the  sign  of  the  modulus 

function  (IP)  each  time  x and  y reversed  signs  at  the 

same  time.  After  each  bipolar  modulus  value  has  been  t 

computed  it  is  placed  in  the  array  DV.  When  the  array  ; 

DV  is  full  (JSTOP  values)  then  the  array  is  rpocessed  1 

and  transmitted  to  the  DUIS  via  the  subroutine  PLOTD. 

After  each  sweep  is  transmitted  the  program  waits  I 

for  the  DUIS  to  return  a status  number,  ISTAT.  If  the  ] 

DUIS  replies  with  ISTAT-1  then  the  sweep  must  be  ] 

retransmitted  since  a transmission  error  was  detected.  | 

If  the  DUIS  replies  with  ISTAT-0  then  the  sweep  was  | 

received  with  no  errors  detected  and  the  processing  of  | 

the  next  sweep  is  begun.  | 

Vhen  all  16  rows  of  data  have  been  transmitted  to  , 

the  DUIS  then  another  16  rows  are  read  in  and  processed 
in  the  same  manner  as  described  above.  If  the  end  \ 

block  NBMAX  is  reached  then  the  program  waits  for  the 
user  'to  transmit  a carriage  return  and  then  the  program 
prints  the  maximum  value  of  the  data  transmitted. 

Next,  the  TSS  PRMFL  is  detached,  control  is  transferred 
to  Statement  # 110,  and  the  user  is  asked  for  another 
file  name. 

i 

c 

'.I 

i 

i 

1 

I 

i 

\ 
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6 . FORTRAN  LISTING 


17  APR  78  1343 

DIMENSION  DV(236) 

COMPLEX  A(160, 16) 

CHARACTER  AFILE *20, STOP *4 , DONE *4 
DATA  I0R/04 00000000000/, STOP/' STOP V 
DATA  DONE/0007040033014/ 

FORMAT(V) 

REQUEST  TSS  FILE  NAME  FROM  USER 

110  WRITE (06, 100)  'INPUT  FILE  NAME  ?' 

READ  100 , AFILE 

IF ( AFILE . EQ . STOP ) CO  TO  600 

CALL  ATTACH (01, AFILE, 1,1, ISTAT,  > 

IF(ISTAT .EQ .IOK .OR .ISTAT .EQ .0)  CO  TO  70 
WRITE(06, 510)  ISTAT 
CO  TO  110 

70  CALL  RANSIZ(01 , 312) 

REQUEST  PARAMETERS  FROM  USER 

WRITE (06 , 100)  'Enter  JWID . JWxDSP , NBMAX ' 

READ  100, JWID, JWIDSP, NBMAX 

IF(JWIDSP .LE .10 .AND . JWIDSP .LE .JWID)  GOTO  75 

JWIDSP-10 

WRITE(06, 100)  'JWIDSP  > 10...  SET  TO  10' 

IF ( JWIDSP. LE. JWID)  COTO  73 
JWIDSP-JWID 

WRITE (06 , 100)  'JWIDSP  > JWID...  SET  JWIDSP-JWID’ 
75  WRITE(06 , 100)  'Enter  FIRST  LREC.ISCN’ 

READ  100,LRECF,ISCN 
TH-0.0 


IF  LRECF  < 1 THEN  TERMINATE  PROCESSING  THIS  FILE 


IF(LRECF .LE .0)  CO  TO  400 
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700 
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JST0P-JVIDSP*16 
NST-1 
NSTP-1S 
CO  TO  200 


SECTION  TO  READ  IN  A STRIP  OF  DATA 


LRECF-LRECF+JWID 
IF(LRECF.CT.NBMAX)  CO  TO  83 
LR1IN-LRECF 

READ ( 01 ' LR1IN ) ( ( A ( J , K > , J-NST , NSTP ) , K-l , 1 6 ) 

NST-NSTP+1 

NSTP-NSTP+16 

LR1IN-LR1IN+1 

IF(NSTP.LE.JSTOP)  GO  TO  190 


SECTION  TO  PROCESS  A STRIP  FOR  OUTPUT 


NST-1 
NSTP-1S 
DO  310  K-l, 16 
IP-1 

DO  300  LL-1, JSTOP 
IF(ISCN.EQ.l)  CO  TO  770 
E-CABS(A(LL,K) ) 

X-REAL(A(LL,K) ) 

Y-AIMAC(A(LL,K) > 

IF(LL.NE.l)  CO  TO  112 

XO-X 

YO-Y 

P1-X0*X 

P2-Y0 *Y 

IF(E .LT .1 .0E-30)  CO  TO  111 

IF(P1  .LT  .0 .0  .AND .P2  .LT  .0.0)  IP— IP 

XO-X 

YO-Y 

DV ( LL  > -SICN (E , IP ) 

CO  TO  300 

DV(LL)-REAL(A(LL,K) ) 

CONTINUE 

CALL  PLOTD(DV, JSTOP, 0.0,1 .0,TH) 
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WAIT  FOR  DUIS  TO  REPLY  WITH  STATUS  CODE 

READ  100, ISTAT 
IF(ISTAT.EQ.l)  CO  TO  700 
310  CONTINUE 
CO  TO  210 
C 

C ALL  PLOT  RECORDS  IN  THE  FILE  HAVE  BEEN 
C SUCCESSFULLY  TRANSMITTED  TO  DUIS 
C 

85  WRITE (06 , 100)  DONE 

READ  100 , AFILE 
WRITE(06 , 100)  *TH-*,TH 
CO  TO  75 

400  CALL  DETACH < 01, ISTAT,  ) 

WRITE (06, 410)  ISTAT 
410  FORMAT ( * DETACH  ISTAT- * , 020 ) 

CO  TO  110 

510  FORMAT ( ‘UNSUCCESSFUL  ATTACH  ISTAT-’, 020) 
600  CALL  EXIT 
STOP 
END 
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PROCRAM  RNDERR 


1 . PURPOSE 


This  program  adds  random  phase  errors  to  the 
elements  of  an  existing  aperture  distribution  stored  in 
a PRMFL.  The  resulting  aperture  may  be  stored  on  the 
input  PRMFL  or  may  be  stored  on  another  user  specified 
PRMFL . 

2 . INPUT  PARAMETERS 


Name 


O/R  T 


Description 


NTYPE 


R I Determines  the  type  of 
distribution . 


MAD1 


JRND 


- 1 Uniform 
■ 3 Gaussian 

The  starting  address  for 
selecting  random  numbers  from 
the  random  number  array. 

(1  £ MAD1  £ 128) 

Random  number  generator 
initialization  constant. 

(0  £ JRND  £ 236  - 1) 


UMEAN 


Number  of  logical  record 
blocks  in  the  x-direction 
(horizontally) . 

Number  of  logical  record 
blocks  in  the  y-direction 
(vertically) , 

Mean  value  of  the  uniform 
distribution  in  degrees. 
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UUEXT 

0 

F 

Width  of  the  uniform 
distribution  in  degrees. 

XMEAN 

0 

F 

Mean  value  of  the  Gaussian 
distribution  in  degrees. 

SIGMA 

0 

F 

Standard  deviation  of  the 
Caussian  distribution  in 
degrees . 

RESTRICTIONS. 

REQUIREMENTS.  MISCELLANEOUS  DATA 

a . The 

value 

of  MAD1 

must  be  in  the  range 

1 

£ MAD1 

£ 128 

b.  The  value  of  JRND  must  be  In  the  range 
0 < JRND  < 236  - 1 

4.  SUBPROGRAMS  REQUIRED 

RRAND 

FXOPT 

5.  THEORY  OF  OPERATION 

The  program  steps  sequentially  through  the  LREND 
(«LRJ*LRK)  input  logical  record  blocks  starting  with 
block  number  1.  For  each  element  in  a block,  a random 
phase  component  is  generated  according  to  the  equation 

PHERR  - RRAND(NTYPE) *0.017453. 

RRAND  is  the  random  number  generator  function  and  the 
constant  0.017453  is  a degree-to-radian  conversion 
factor.  With  the  phase  error  expressed  in  radians  a 
complex  representation  C is  calculated  as 

C - CMPLX (COS (PHERR), SIN (PHERR) ) . 

152 


, -V  - 


PAAS-CPD-RNDERR 
1 HAY  78 


This  number  is  then  multiplied  by  the  value  of  the 
antenna  aperture  element  to  be  modified.  The 
multiplication  is  defined  by  the  equation 

A<I>  - A(I)«C. 

Note  that  this  is  a complex  multiplication  by  a unit 
vector  and  that  no  change  in  energy  is  introduced  as  a 
result  of  the  phase  errors.  This  process  is  repeated 
for  all  of  the  elements  in  each  block.  As  the  elements 
are  processed  each  block  is  written  off  to  the 
appropriate  PRMFL,  either  the  input  PRMFL  or  another 
that  is  chosen  by  the  user. 
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WRITEC08 , 550)  ISTAT 

CO  TO  580 

LREND-LRJ *LRK 

DO  700  LRIN-1,LREND 

READ(Ol'LRIN)  A 

DO  650  1-1,256 

PHERR-RRAND ( NT YPE ) *0 .017453 

C-CMPLX(COS(PHERR) , SIN(PHERR) ) 

A(I)-A(I)*C 

WRITE (LROUT'LRIN)  A 

CONTINUE 

CALL  DETACH <01, 1ST AT,  ) 

WRITE (06, 7 50)  ISTAT 

CALL  DETACH (02, ISTAT,  ) 

WRITE(06 , 750)  ISTAT 

FORMAT ( 'DETACH  ISTAT- * , 020) 

CALL  EXIT 

STOP 

END 
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1. 


FUNCTION  RRAND 


PURPOSE 

This  function  generates  random  numbers  for  use  in 
various  subprograms  of  PAAS.  Samples  from  the  uniform, 
Gaussian,  and  Rayleigh  distributions  can  be  generated. 
This  function  is  based  on  the  function  RRAND  used  in 
RADSIM . 


2.  INPUT  PARAMETERS  (Common  Area  BLKRND) 


Name  Q/R  T 

NADI  R I 

JRND  R I 

XMEAN  R F 

SIC2SQ  R F 


UL  R F 


UEXT  R F 


Description 

Random  Number  Table  pointer 

Random  integer  from  previous 
execution  of  RRAND 

Mean  value  of  the  Gaussian 
distribution 

An  intermediate  parameter  used 
in  computing  Gaussian  and  Rayleigh 
distributions 

An  intermediate  parameter  used 
in  determining  the  uniform 
distribution  mean  value 

An  intermediate  parameter  used 
in  determining  the  uniform 
distribution  width 
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3.  CALLING  SEQUENCE: 

VAR  - RRAND  (NTYPE) 

NTYPE  R I Control  Integer  which  specifies 

the  type  of  distribution  to  be 
generated . 

NTYPE  - 1 Uniform  distribution 

(floating  point  output) 
- 2 Rayleigh 
« 3 Gaussian 


VAR  contains  the  random  sample  generated  by  the 
function  from  the  NTYPE  probability  distribution. 

4.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  Before  any  call  can  be  made  to  the  function  RRAND 
the  labeled  common  area  must  be  loaded  as  follows: 

NADI  Any  integer  on  the  interval:  0 £ HAD1  £ 127 

JRND  Any  integer  on  the  interval:  0 £ JRND  £ 2s 5 

XMEAN  Mean  value  of  Gaussian  distribution 

SIC2SQ  - -2.0*SICMA*SICMA 


Where:  SIGMA  is  the  standard  deviation  of  j | 

the  Caussian  or  Rayleigh  ; fjj 

distribution  1 


UL  - UMEAN  - 0.5  * UEXT 

Where:  UMEAN  is  the  mean  value  of  the 

\ 

uniform  distribution 

UEXT  is  the  width  of  the  uniform 
distribution 
158 


M 

I 

i 

I 

M 

I 

1 

% 

I 

I 


-s 

% 

% 


"mWKWSSMBSSSRM^ 


PAAS-CPD-RRAND 
16  JUN  77 


UEKT  - UEXT/235 


b . 


For*  convenience  and  to  minimize  program  steps,  the 
array  IRAND  was  equlvalenced  to  the  arrays  NRND1 
and  NRND2,  but  displaced  by  one  location.  This 
structure  allows  an  address  of  zero  to  be  used, 
i.e.,  an  address  of  zero  will  access  IRAND  (0), 
which  overlays  NRND(l).  If  this  were  not  done,  a 
test  would  have  to  be  performed  on  MAD1  to  ensure 
that  an  address  of  zero  did  not  occur. 


c . 


The  distribution  transformations  used  herein  are 
documented  in  the  following  reference: 


Robert  E.  Hachol  (ed.).  System  Engineering 
Handbook.  McGraw-Hl 11 , N.Y.,  1963,  pp.  40-28, 
40-29 . 


d.  Referenced  labeled  common  areas: 


BLKRND 


e.  Source  PRMFL: 


/PLARY 


5-  SUBPROGRAMS  REQUIRED 


FLD 

FLOAT 

SQRT 


6 . THEORY  OF  OPERATION 


For  each  call  to  the  function  RRAND  a number  KRND 
is  selected  from  the  random  number  table,  IRAND.  The 
address  of  the  number  selected  from  the  table  is  MAD1, 
a random  number.  The  number  KRND  is  added  to  the 
product  of  JRND  and  IMULT  and  stored  in  IRND.  The 
variable  JRND  is  the  random  number  which  was  generated 
by  the  previous  execution  of  the  function  and  IHULT  is 

159 


iV  >'»^.<fV*'<'J*  ; y7,  V*  "yi, ,”  —■  vsAv^->f?^<57AVv>?-’  '^'^fu-yi*-.\''v»-;^-  '+}’“'-■**> — !f^'‘'.?'^W^^^!f}!f'*f.'(,!^ 


PAAS-CPD-RRAND 
IB  JUN  77 


a multiplier  chosen  because  it  results  in  good  bit 
scrambling.  The  multiplication  of  JRND  by  IMULT  causes 
the  bits  of  JRND  to  be  scrambled.  This  scrambled  word 
is  then  added  to  KRND,  just  retrieved  from  the  table  to 
form  the  nev  random  number , IRND.  The  sign  bit  of  XRND 
is  set  to  zero  to  ensure  a positive  number.  The  random 
number  XRND  just  generated  is  an  integer  having  a 
uniform  distribution  from  0 to  235-l.  XRND  is  placed 
in  the  random  number  table  location  previously  occupied 
by  the  KRND.  In  this  manner  the  random  number  table  is 
updated  by  generating  nev  random  numbers  and  inserting 
them  into  the  table.  From  the  random  number  XRND,  7 
bits  are  selected  to  determine  the  nev  address  MAD1  to 
be  used  in  the  next  call  to  the  function.  The  7-bit 
address  field  allovs  the  addresses  to  range  from  0 to 
127.  Once  the  random  numcer  is  generated  JRND  is  set 
equal  to  XRND  for  use  in  the  subsequent  executions  of 
the  function.  In  order  to  convert  this  number  to  a 
floating  point  number,  r,  having  a uniform  distribution 
from  0 to  1.0,  the  folioving  conversion  is  used: 


r - FLOAT(IRND)/235 


From  this  uniform  distribution  other  probability 
distributions  can  be  generated  by  using  transformations 
vhlch  map  a uniform  distribution  into  the  desired 
distribution.  The  folioving  is  a list  of  the 
transformations  used  in  this  function: 


a.  Uniform  distribution 


p(x)  ■ l/(b-a) 


, <a  £ x i b) 


xn  - (b  - a)rn  + a 


b. 


Rayleigh  distribution 

p(x)  - (x/c2)  exp(-x2/2a8)  ,<x  2 0) 


xn  - -2cra  In  r„ 
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c.  Gaussian  distribution 

p<x>  - <1/7ttc>  exp[(x-H>8/2c8]  ,<-«  < x < ») 

xn  - 2c 8 In  rn  * Cos  2nrn+l 
y„  - 2 c 8 In  rn  * Sin  2TTrn  + l 
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7 . FORTRAN  LISTING; 


FUNCTION  RRAND(NTYPE) 

C0MM0N/BLKRND/MAD1 , JRND , XMEAN , SIG2SQ , UL , UEXT 

DIMENSION  IRAND (128), NRND1 (65), NRND2 (64) 

DATA  NRND1/15134181997, 27509664464,30323512272, 

& 14051007893,16402190290,26306990212,11260717646. 

& 16801629773,11849273156,19404991345,06977712830, 

& 02883434137,33025570091,11012391622,13411365861, 

& 31267410086.13462139250,26463885902,24219774296, 

& 11557820695,30512809719,12630506319,17722780814, 

& 04722597022,16900280091.16243824041,16388044606, 

& 26212698408,13570004754,11188309528,29134237821, 

& 13164942096.29908968258,03564986686,24513426529, 

8i  25262307992,16416251777,32749370939.21116178576. 

& 19395173043,20743061171,21319359579,19074491967, 

& 19244390324,08846123356,27142309994,15825176938, 

& 16410917813,23416520791,28825638452,10800745449, 

& 01702686304,17006458873,16841482774,26473264721, 

& 17160292937,29260744156,8883554486,3669953728, 

& 16068801392,5883873859,14824731880,18081451748, 

& 8160418880,30068227383/ 

DATA  NRND2/12068158044, 06847664659, 15416782660, 

& 25052201840,13988647055,01734737408,07289355507, 

& 28120785669,32320902560,19471392797,07683759917, 

& 24386072834,29317493972,07114843643,16232718423, 

8r  29170604246,26866574818,20335880812,14861357546, 

& 25072568248,31374670078,13676667951,30463132192, 

& 20172084006,16184261842,14974210467,10283018420, 

& 18310257399,18938188207,01286074697,19662214195, 

& 01577045480,16742867695,11686848767,18174114680, 

8t  30892487160,30892487160,28360949700,33368415709, 

& 17235921632,25322444850,30007056175,13488881553, 

& 30224148581,07655423387,32626402591,13101024674, 

& 30533512969,07218771539,00229536870,29198604401, 

& 33122308420,29107616508,16534467415,3669736170, 

& 3491463822,5804776974,30256545186,10832795361, 


162 


4. 


«8  «9 


PAAS-CPD-RRAND 
IB  JUN  77 


18174114680,10556707007 , 10140208896 , 9779017119 , 
19382343178/ 

DATA  IMULT/1220703 125/, IMAX/4 294967296/ 

DATA  N2P16/65536/,CI/2 . 9103830E-11/ 

EQUIVALENCE  ( IR AND <1 ) , NRND1 ( 2 ) ) , ( IRAND (65), NRND2 ( 1 ) ) 

10  KRND-IRAND(MADl) 

IRND-KRND+ JRND  *IMULT 
IRND-FLD (1,35, IRND ) 

IRAND (MAD1) -IRND 
NADI -FLD (15,7, IRND ) 

CO  TO  ( 200 , 300 , 400 ) , NTYPE 

200  RRAND-FLOAT ( IRND ) *UEXT+UL 
JRND-IRND 
RETURN 

300  RR AND-SQRT ( SIC 2SQ  * ALOC ( FLOAT ( IRND ) *CI) ) 

JRND-IRND 

RETURN 

400  I1-FLD( 1,17, KRND)-N2P1G 
I2-FLD(18 , 17 ,KRND)-N2P16 
IS-I1 *11+12 *12 
IF(IS.LT.IMAX)  CO  TO  20 
JRND-IRND 
CO  TO  10 

20  S-l. 0/FLOAT (IS) 

VCOS-S  *FL0AT(I1 *11-12  *12 ) 

VSINE-S*2 .0*FLOAT(I1 *12) 

DUM-SQRT (SIC2SQ*ALOC( FLOAT ( JRND) *CI) ) 

RRAND-DUM *VCOS+XKEAN 

DUM-DUM  *VSINE+XMEAN 

JRND-IRND 

RETURN 

END 
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PROGRAM  RTI4 


1 . PURPOSE 


The  purpose  of  this  subroutine  is  to  produce  a 
compact  representation  of  three-dimensional  data,  e.g 
antenna  far  field  pattern,  antenna  aperture 
illumination  distributions,  and  radar  ambiguity 
diagrams.  The  procedure  used  herein  is  to  represent 
the  modulus  of  each  sample  with  an  alphanumeric 
character . 


2 . ISJPU7  PARAMETERS 


Name 


O/R 


Description 


AFILE  R C 


The  name  of  the  TSS  PRMFL  which 
the  user  wishes  to  process 


FLOOR 


R F 


The  reference  in  dB  below  which 
everything  is  represented  by 
dashes(-)  on  the  RTI  plot 


YINC 


R F 


The  increment  in  dB  between  each 
successive  alphanumeric  symbol. 
See  Table  RTI4-1 . 


JWID 


The  data  array  width  in  blocks. 
This  parameter  corresponds  to 
LRJVID  in  the  program  FFT2DX 
or  LRJ  in  the  programs  PLARY, 
RNDERR,  or  FILMOD. 


NBMAX 


The  number  of  the  last  block 
to  be  displayed 


LREC 


The  number  of  the  first  record 
to  be  processed 
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The  width  of  the  output 
character  matrix. 


3.  RESTRICTIONS.  REQUIREMENTS . MISCELLANEOUS  DATA 


a.  The  maximum  value  of  the  parameter  IWD  is  128.  In 
addition,  the  value  of  IVD  should  be  less  than  or 
equal  to  JWID*16. 


b.  Two  examples  of  the  relative  dB  values  for 

specified  values  of  FLOOR  and  YINC  are  shown  in 
Table  RTI4-1. 


SUBPROGRAMS  REQUIRED 


FPARAM 

ATTACH 

RANSIZ 

DETACH 


L&k£QtHd£&ri 
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TABLE  RTI4-1 


EXAMPLES  OF  RELATIVE  VALUE  OF  LETTERS,  NUMBERS,  AND  SYMBOLS 


FLOOR  - -20. 0 

20.0 

t 

30.0 

FLOOR  - 10.0 

YINC  - 1.0 

19  .0 

• 

29.5 

YINC  - 0.5 

18  .0 

* 

29  .0 

17.0 

+ 

28  .5 

16  .0 

0 

28  .0 

15  .0 

1 

27.5 

14.0 

2 

27  .0 

13.0 

3 

26.5 

12.0 

4 

26.0 

11.0 

5 

25.5 

10.0 

6 

25  .0 

9.0 

7 

24.5 

8.0 

8 

24.0 

7.0 

9 

23.5 

6.0 

A 

23.0 

5.0 

B 

22.5 

4 .0 

C 

22.0 

3.0 

D 

21.5 

2 .0 

E 

21.0 

1 .0 

F 

20.5 

0.0 

C 

20.0 

-1.0 

H 

19.5 

-2  .0 

I 

19.0 

-3  .0 

J 

18.5 

dB  level 

-4  .0 

K 

18  .0 

dB  level 

-5.0 

L 

17.5 

-6  .0 

M 

17.0 

-7.0 

N 

16.5 

-8  .0 

0 

16  .0 

-9  .0 

P 

15  .5 

-10.0 

Q 

15.0 

-11  .0 

R 

14  .5 

-12  .0 

S 

14  .0 

-13  .0 

T 

13  .5 

-14  .0 

U 

13  .0 

-15  .0 

V 

12  .5 
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5.  THEORY  OF  OPERATION 


After  startup  a message  is  printed  to  the  user 
requesting  the  name  of  the  TSS  PRMFL  that  Is  to  be 
processed.  The  user  reply  is  tested  to  see  if  it  is 
the  word  'STOP*.  If  so,  program  execution  is 
terminated.  Otherwise,  an  attempt  is  made  to  access 
the  file  via  the  system  subroutine,  ATTACH.  If  the 
file  cannot  be  accessed  the  system  error  code  is 
printed  and  the  user  is  requested  to  try  again. 


Assuming  the  user  has  successfully  accessed  a file, 
the  program  then  requests  the  following  parameters: 


FLOOR,  YINC,  JVID,  NBMAX 


Then  the  following  parameters  are  requested: 
LREC , IVD 


The  parameter  LREC  is  tested  to  see  if  it  is  zero  or 
negative.  If  so,  the  file  is  deaccessed  via  a call  to 
the  system  subroutine  DETACH.  The  program  then  returns 
to  the  procedure  for  accessing  another  file  (Statement 
# 99).  The  parameter  IVD  is  tested  to  ensure  that  its 
value  is  in  the  range  from  1 to  128.  The  number  of 
blocks  (NREC)  required  for  the  specified  display  width 
(IVD)  is  computed  and  compared  to  JVID.  If  NREC  > JVID 
then  NREC  is  set  equal  to  JVID  and  IVD  is  set  equal  to 
NREfc*16.  In  other  words,  IVD  is  made  as  large  as 
possible  for  the  set  of  data  to  be  processed. 


The  NREC  blocks  starting  with  LREC  are  read  in  and 
processed  in  the  following  manner.  The  modulus  of  each 
complex  valued  sample  is  computed  and  converted  to  dB. 
Then  the  dB  value  (KM)  is  mapped  to  an  integer  number 
(IADD)  on  the  interval  from  1 to  41  by  the  following 
procedure : 
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IADD-IFIXC <XM-FLOOR)/YINC+0 .5)+l 
IF(IADD  .GT .40)  IADD-41 
IF(IADD.LT.l)  IADD-1 

The  Integer  numbers  41  and  1 correspond  to  S and  -, 
respectively,  and  represent  values  either  too  large  or 
too  small  to  be  displayed  for  the  set  of  parameters, 
FLOOR  and  YINC,  specified  by  the  user.  The  Integer 
numbers  are  used  to  'pull*  the  corresponding  character 
ASCII  code  from  the  character  table,  CTABL.  These 
characters  are  then  stored  in  the  output  character 
matrix,  XRTI. 


Once  the  NREC  blocks  have  been  processed,  the 
character  matrix  is  transmitted  to  the  user.  This 
character  matrix  contains  16  rows  of  data.  The  above 
procedure  is  repeated  until  the  record  number  NBMAX  is 
encountered.  Control  then  returns  to  the  statement  (# 
75)  requesting  the  first  logical  record  to  be 
processed . 
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6 . FORTRAN  LISTING 


C 

C 26  APR  78  0843 

C 

COMPLEX  TEMP (16, 16) 

DIMENSION  XRTI(32 , 16) , CTABL(41) 

DATA  CTABL/0137, 0132, 0131, 0130, 0127, 0126, 0125, 

8.  0124,0123,0122,0121,0120,0117,0116,0113,0114,0113, 

& 0112,0111,01 10,0107 , 0106 , 0105 , 0104 , 0103 , 0102 , 0101 , 

& 071,070,067,066,065,064,063,062,061,060,053, 

& 052,056,044/ 

CHARACTER  AFILE *20, STOP *4 

DATA  IOK/O400000000000/, STOP/' STOP"/ 

C 

C REQUEST  TSS  FILE  NAME  FROM  USER 
C 

CALL  FP ARAM (1,130) 

99  WRITE(06, 105) 

105  FORMAT ( ‘INPUT  DESIRED  FILE  NAME') 

READ(05,50)  AFILE 
50  FORMAT (A20) 

IF ( AFILE. EQ. STOP)  GOTO  310 
CALL  ATTACH ( 01 , AFILE, 1 , 1,ISTAT , ) 

IFdSTAT .EQ .IOK .OR  .ISTAT .0)  GO  TO  101 
WRITE(06 , 410)  ISTAT 
GO  TO  99 
C 

C REQUEST  PARAMETERS  FROM  USER 
C 

101  WRITE(06 , 100) 

100  FORMAT ( ‘ FLOOR , YINC , JWID , NBM AX  ' ) 

READ  1 15, FLOOR. YINC, JWID, NBMAX 
CALL  RANSIZ(01,512) 

75  WRITE(06 , 110) 

110  FORMAT( ‘Enter  FIRST  LREC,  DISPLAY  WIDTH’) 

READ  115,LREC,IWD 
115  FORMAT(V) 

C 
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C 

C 


IF  LREC  < 1 THEN  TERMINATE  PROCESSING  THIS  FILE 

IF(LREC.LE.0>  CO  TO  700 

IF(IWD .LE .128 .AND .IWD .GT .0)  GOTO  150 

IWD-128 

WRITE (06 ,115)  'DISPLAY  WIDTH  TOO  BIG.  SET  TO  128* 
150  NREC- ( IWD+15 )/16 

IF ( NREC . LE . J WID ) GOTO  160 

NREC-JWID 

IWD-NREC*16 

WRITE(06 , 115)  'DISPLAY  WIDTH  > AVAIL  DATA..,' 

* ' CHANGED  TO:  IWD-JWID*16' 

BEGIN  PROCESSING  FOR  A STRIP 

160  IWD4-IWD/4 

LREND-LREC+NREC-1 
IF(LREND.CT.NBMAX)  GOTO  75 
NST--1 

DO  800  IREAD-LREC,LREND 
READ (01 ' IREAD)  TEMP 

PROCEDURE  TO  PROCESS  ONE  BLOCK  OF  DATA 

DO  200  J-1,16 
IBIT“MOD( J-l , 4 ) *9 
IWORD-l+(NST+J)/4 
DO  300  IR-1,16 
XM-CABS(TEMP(J,IR)> 

IF(XM .LT .1 .0E-10)  CO  TO  120 
XM-20.0*ALOG10(XM> 

GO  TO  121 

120  XM— 100.0 

121  IADD«IFIX((XM-FLOOR>/YINC+0.5)+l 
IF(IADD .CT .40)  I ADD- 41 
IF(IADD.LT.l)  IADD-1 

FLD(IBIT , 9 , XRTI ( IWORD , IR ) ) -CTABL ( IADD ) 

300  CONTINUE 
200  CONTINUE 
NST-NST+16 
800  CONTINUE 


I 

1 


c 

c 
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SECTION  TO  TRANSMIT  OUTPUT  TO  USER 


DO  450  IR-1,16 

WRITE ( 6 , 1000) (XRTHJ , IR) . J-l , IWD4 ) 
1000  FORMAT (1H  , 32A4) 

450  CONTINUE 

LREC-LREC+JWID 
COTO  160 


C 

C * DISCONNECT  FROM  USER  TSS  FILE 
C 


700  CALL  DETACH < 01, ISTAT,  ) 

500  WRITE (06,510)  ISTAT 
CO  TO  93 

410  FORMAT( ‘UNSUCCESSFUL  ATTACH*, * ISTAT-  *.020) 
510  FORMATC  DETACH  ISTAT  - *,020) 

310  CALL  EXIT 
STOP 
END 
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PROCRAM  TBLS 


PURPOSE 


The  program  TBLS  computes  and  tabulates  the  sampled 
values  of  selected  weighting  functions.  The  program 
also  generates  data  which  are  used  in  checking 
probability  density  functions  of  space  tapered  arrays 
estimated  by  program  PDFESTR. 


INPUT  PARAMETERS 


Name 


O/R  T 


NTYPE 


R I 


WTPED 


0 F 


NWTPQW 


0 I 


IRAD 


R I 


t£M  iid 


Description 


Determines  the  weighting 
function  tabulated. 


2 

3 

4 

5 

6 

7 

8 


Cosine  on  a pedestal  to  a 

power 

Blackman 

Kaiser 

Bartlett  or  triangular 

Taylor 

Bessel 

Cubic 

Bayllss 


The  height  of  the  pedestal  for 
cosine  on  a pedestal  to  a power 
weighting . 


The  power  of  the  cosine  function 
for  cosine  cn  a pedestal  to  a 
power  weighting. 


The  radius  (or  half  span  of  a 
linear  array),  in  units  of 
elements  of  the  array. 
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WTRAD 


The  radius  (or  half  span  of  a 
linear  array)  of  the  weighting 
function . 


WKASIR 


BESEDC 


BESCAL 


CUBK 


IDB 


NBAR 


The  Kaiser  variable  for  the 
trade-off  between  main  lobe 
width  and  side  lobe  amplitude. 

The  radius  scaling  constant  for 
the  Bessel  weighting. 

The  maximum  weight  amplitude  for 
the  Bessel  weight.  This 
corresponds  to  a radius  of  zero 

The  weighting  amplitude  scaling 
constant  for  the  cubic  weight. 

The  design  side  lobe  amplitude 
in  dB  for  the  Taylor  or  Bayliss 
weighting . 

The  number  of  zeros  used  to 
approximate  the  Dolph-Chebyschev 
weighting  distribution  in  the 
Taylor  or  Bayliss  weighting. 


3.  RESTRICTIONS,  REQUIREMENTS , MISCELLANEOUS  DATA 

a.  The  range  of  NBAR  must  be 

3 £ NBAR  £ 20 
for  the  Taylor  weighting. 

b.  The  range  of  NBAR  must  be 

3 £ NBAR  £ 19 

for  the  Bayliss  weighting. 
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SUBPROCRAMS  REQUIRED 

EXPND 

BESS 

CAM 

WEIGHT 

THEORY  OF  OPERATION 

This  program  is  a driver  for  the  subroutine  WEIGHT. 
A loop  is  set  up  that  starts  at  zero  and  goes  to  IRAD 
in  steps  of  one.  In  each  pass  of  the  loop  a call  to 
WEIGHT  is  made  and  a value  of  the  selected  weighting 
function  is  returned  for  the  particular  radius  cell. 

The  values  are  stored  and  printed  on  the  time-sharing 
terminal  in  a tabular  fashion. 
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6.  FORTRAN  LISTINC 


COMMON  NTYPE , WTPED , NWTPOW , WKASIR , F ( 2® ) , B ( 20 ) , ANC , 
NBAR , BESCAL , CUBK1 , PII 2 , BESS1 , IAZ , XKK , VMAX , 
BESEDC 
INTEGER  DB 

DIMENSION  U < 2®  > , OUT ( 205 ) , SIGG ( 5 ) , IX < 20 ) , BZEROl ( 2® ) 
DIMENSION  Z(20> 

DATA  U/l. 2196699 ,2. 2331306 ,3. 2383154, 4. 2410628, 

9. 242 7643, 6. 2439216,7. 2447598, 8.2493948, 
9.2498927,10.2462933.11.246624,12.246900, 
13.247131,14.247334,15.247508,16.247663, 
17.247796,18.247920,19.248027,20.248125/ 

DATA  BZERO1/0. 586067, 1.6970509, 2. 7171939, 3. 726137, 
4.7312271,5.7345205,6 .7368281,7.7385355, 
8.7398505,9.7408945,10.7417435,11 .7424475, 
12.7430408,13.7435477,14.7439856,15.7443679, 

16 .7447044,17 .745003,18.7452697,19 .7455093/ 
WRITEC06 , 50)  'PDFESTR  DATA  OR  TABLES?  (0  OR  1)* 
READ  50, IAZ1 
IAZ-1 
ANC-0.0 
PI-3.1415926 
PII2-2 .0/<PI*PI) 

BESS1-1 . 0/BESS (0 .0,0.0) 

WRITE (06, 50)  'NTYPE  * 

READ  50, NTYPE 
IF(NTYPE.EQ.-l)  CO  TO  1830 

CO  TO  (1100,1200,1300,1200,1500,1360,1370,1500), 
NTYPE 

1100  WRITE(06 , 50)  'WTPED, NWTPOW, IR AD, WTRAD' 

READ  5®, WTPED, NWTPOW, IRAD, WTRAD 
CO  TO  1600 

1200  WRITE(06 , 50)  'IRAD, WTRAD' 

READ  50, IRAD, WTRAD 
CO  TO  1600 

1300  WRITE(06 , 50)  'WKASIR, IRAD, WTRAD' 

READ  50, WKASIR, IRAD, WTRAD 
CO  TO  1600 

1360  WRITE(®6 , 50)  'BESEDG, BESCAL, IRAD, WTRAD* 


1850 
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READ  50 , BESEDC , BESCAL , IRAD , WTRAD 
CO  TO  1600 

1370  WRITE <06, 50)  ‘CUBK, IRAD, WTRAD * 

READ  50, CUBK, IRAD, WTRAD 
XX-SQRT< WTRAD *WTR AD/3 .0) 

CUBK 1 -CUBK/ABS ( XX  * ( XX- WTRAD ) *<XX+VTRAD) ) 

CO  TO  1600 

1500  WRITEC06 , 50)  ‘IRAD, WTRAD* 

READ  50, IRAD, WTRAD 

WRITEC06 , 50)  ‘ALL  OR  SINCLE  DB  LEVEL?  <0  OR  1)‘ 
READ  50, IA 
IF<IA.EQ.0)  CO  TO  60 
WRITE <06, 50)  ‘IDB  * 

READ  50, IDB 
IDB1-IDB 
IDB2-IDB1-1 
IDB3-1 

WRITE<06 , 50)  ‘ALL  OR  SINCLE  NBAR?  <0  OR  1)‘ 
READ  50, IAX 
IF<IAX .EQ .0)  CO  TO  70 
410  WRITE<06 , 50)  *NBAR‘ 

READ  50, INBAR 
IBAR1-INBAR 
IBAR2-IBAR1-1 
CO  TO  80 
60  IDB1-20 

IDB2-80 
IDB3-5 

70  IBAR1-3 
IBAR2-20 
80  CONTINUE 
50  FORMATCV) 

DO  20  1-1,20 
IX<I)-I 
20  CONTINUE 
Il-IRAD+1 

DO  300  III-IDB1 , IDB2 , IDB3 
DB-FLOAT<III) 

IOUT-3 

SICl-l 

XA-10.0** <111/20.0) 


IS  JUN  77 


A- ALOC ( XA+SQRT ( XA  *XA-1 ) )/PI 

AASQ«A*A 

IT-0 

IT1-0 

NSCFL-0 

DO  200  NBAR-IBAR1 , IBAR2 

IF ( NB  AR . EQ . 20 . AND . NTYPE . EQ . 8 ) CO  TO  200 

IFOiTYPE .EQ .8)  CO  TO  i900 

IFCNSGFL .NE  .0)  CO  TO  30 

SICP1 -U ( NBAR+1 >/SQRT ( AASQ+ < NBAR+0 . 5 ) * *2 ) 

30  SICC < SIG1 ) -U  < NBAR  >/SQRT ( AASQ+ ( NBAR-0 . 5 ) * *2 ) 

SC-SICC(SICl) 

SGSQ-SG*SG 

IF(SC.LE.SICP1 .AND.IAX.EQ.l)  CO  TO  400 

IF(SG.LE.SIGPl)  CO  TO  200 

NSCFL-1 

SICP1-0.0 

DO  292  I-l.NBAR-l 

FNUM-1 .0 

FDNM-1 .0 

T-U(I)*U<I) 

XII— 0.9 

DO  254  II-l , NBAR-1 
XII-XII+1 .0 

FNUM-FNUM  * <1 . 0-T/<  SCSQ  *< AASQ+ ( XII *XII ) ) ) ) 

IF (II  .EQ .1)  CO  TO  254 
FDNM-FDNM  *< 1 . 0-T/(U  < II ) *U ( II ) ) ) 

254  CONTINUE 

ARC-PI *U( I) 

F ( I ) —BESS  ( 0 . 0 , ARC ) *FNUM/FDNM 
F(I)-F(I)/(BESS(0 .0, ARC) **2) 

252  CONTINUE 
CO  TO  1910 

1900  Z(l>-0.9858302+0. 0333885 *DB+0 . 00014064 *DB*DB 
& -0 .0000019  *DB*DB*DB+0 .00000001 *DB*DB *DB*DB 

Z(1)-Z(1)*Z(1) 

Z(2)-2. 00337487+0. 01141548 *DB+0. 0004159 *DB*DB 
& -0 . 00000373  *DB  *DB  *DB+0 . 00000001 *DB *DB  *DB *DB 

Z(2)-Z(2)*Z(2) 

Z(3)-3 .00636321+0 .00683394  *DB+0 .00029281 *DB*DB 
& -0 . 000001 6 1 *DB *DB *DB 
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Z(3)-Z<3)*Z<3> 

Z( 4 )-4 .00518423+0 .00501795  *DB+0 .00021735  *DB*DB 
& -0. ©000008 8 *DB*DB*DB 

Z<4)«Z(4)*Z<4> 

DO  761  1-5, NBAR+1 
Z<I)-AASQ+I*I 
761  CONTINUE 

IF(NSCFL.NE.O)  CO  TO  1930 
SICP1 -BZEROl ( NBAR+2 >/SQRT ( Z ( NBAR+1 ) ) 

1930  SIGC( SIC 1 ) -BZEROl < NBAR+1 )/SQRT( Z (NBAR) ) 
SC-SICC(SICl) 

SCSQ-SC*SC 

IF(SC.LE.SICPl.AND.IAX.EQ.l)  CO  TO  400 

IF(SC.LE.SICPl)  CO  TO  200 

NSCFL-1 

SICP1-0.0 

DO  1940  I- 1, NBAR 

FNUM-1 .0 

T-BZEROl < I ) ‘BZEROl < I ) 

FDNM-1 .0-T/C BZEROl ( 1 > *BZER01 ( 1 ) ) 

IF(I.EQ.l)  FDNM-1 .0 
DO  1950  II«1,NBAR-1 
FNUM-FNUM*(1 .0-T/(SCSQ*Z(II> > > 

IF(I-l.EQ.II)  CO  TO  1930 

FDNM-FDNM  *( 1 .0-T/(BZEROl ( II+l ) *BZER01 ( II+l ) ) ) 

1950  CONTINUE 

B(I)»(2 .0*T/BESS<1 .0, BZEROl (I) *PI ) ) *FNUM/FDNM 
1940  CONTINUE 

P0-0. 4797212+0. 01456692 *<DB)-0. 00018739 *<DB*DB) 

& +0 . 000002 18  *(DB*DB*DB)-0 . 00000001  * ( DB  *DB *DB *DB ) 

P0-P0*SC 
P0SQ-P0 *P0 
PIP0SQ-P0SQ  *PI *PI 
FNUM-1 .0 

FDNM-1 .0-P0SQ/(BZEROlCl) *BZER01C1>) 

DO  772  I-1,NBAR~1 

FNUM-FNUM * ( 1 . 0-P0SQ/( SCSQ  *Z  < I ) ) ) 

FDNM-FDNM  *< 1 .0-P0SQ/< BZEROl  < 1+ 1 ) *BZER01  ( 1+1 ) ) •) 

772  CONTINUE 

C-PIP0SQ-1 .0 

C-C  *BESS ( 1 . 0 , P0  *PI > *FNUM/FDNM 
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C-l .0/C 

DO  773  I-1,NBAR 
B(I)--B(I) *C 
773  CONTINUE 
1910  SIG1-SIC1+1 

IF(IAZl.EQ.l)  CO  TO  5S 
WMAX-0.0 
DO  51  1-1 , IRAD+1 
RAD-FLOAT (1-1) 

CALL  WEIGHT (R AD, VTR AD, VFUNC) 

WMAX-AMAX1 ( WMAX , WFUNC  > 

51  CONTINUE 

55  IT-IT+1 

DO  100  K-l, IRAD+1 
RAD-FLOAT (K-l > -0 . 5 * ( I AZ1-1 ) 

CALL  VEICHT  < RAD , WTR AD , WFUNC ) 

IFCIAZ1 .EQ.0)  WFUNC-WFUNC*2 . 0 *PI *RAD/WHAX 

OUT (IOUT) -WFUNC 

IOUT-IOUT+1 

IF(IAX.EQ.l. AND .K.EQ. IRAD+1)  CO  TO  500 
IF<NTYPE.NE.8)  CO  TO  56 

IF ( ( MOD ( IT , 5 ) . EQ . 0 . OR . NB AR . EQ . 1 9 ) . AND . K . EQ . IRAD+ 1 ) 
& CO  TO  120 

CO  TO  100 

56  IF( (MOD(IT, 5) .EQ .0 .OR  .NBAR .EQ .20) . AND. K .EQ .IRAD+1) 

& CO  TO  120 

CO  TO  100 

500  WRITE(06, 995)  III, A 
WRITE(06 , 997)  NBAR 
997  FORMATC  NBAR-*, 4X, 12,//) 

WRITE<06 , 996)  SICC(l) 

996  FORMATC  SICMA- * ,/,  3X,F12 .10,//) 

DO  975  LL-1, IRAD+1 
975  WRITE(06 , 994)  LL-l,OUT(LL) 

994  FORMAT(I3,F12 .10) 

CO  TO  100 

120  m-ITl+l 

IF(ITl.CT.l)  CO  TO  130 
WRITE(06 , 995)  III, A 

995  FORMAT <//,  *DB- M2  ,//'A- ' ,F12  .10,//) 

130  WRITE (06 , 998 ) (IX ( I > , I-NBAR-IT+1 , NBAR) 
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998  FORMAT ( 'NBAR- ‘ , 4H , 12 , 4 (10X, 12 ) ) 

WRITE (06 , 1998 ) 

WRITE (06 , 993 ) (SIGG( JJ) , JJ-1 , IT) 

993  FORMAT ( 'SIGMA- * 3X, 5<F12 .10) ) 

WRITE(06 , 1998) 

1998  FORMAT (/) 

ITT-5 

IF (NTYPE .NE  .8)  CO  TO  57 

IF  < NBAR . EQ . 1 9 . AND . IT . NE . 5 ) ITT-MDD ( IT , 5 ) 

CO  TO  58 

57  IF ( NBAR . EQ . 20 . AND . IT . NE . 5 ) ITT-MOD(IT, 5) 

58  DO  75  LL-1 , IRAD+1 

75  WRITEC06 , 999)  tL-1. (OUT(LL+I-Il) , 1-11,11 *ITT, II) 

999  FORMAT (13, 5 (F12. 10)) 

WRITE(06, 1999) 

1999  FORMAT (//) 

IOUT-1 

IT-0 

SIC1-1 

100  CONTINUE 
200  CONTINUE 
300  CONTINUE 
GO  TO  420 

400  WRITE(06 , 50)  ’INVALID  VALUE  FOR  NBAR* 

CO  TO  410 
1600  WMAX-1.0 

IF(NTYPE.EQ.B)  WMAX-BESCAL 
IF(NTYPE .EQ .7)  VMAX-CUBK 
DO  1700  K-l, IRAD+1 
RAD-FLOAT ( K- 1 ) -0 . 5 * ( I AZ 1 - 1 ) 

CALL  WEIGHT ( R AD , WTR AD , WFUNC ) 

IF(IAZ1 .EQ .0)  WFUNC-WFUNC *2 . 0 *PI *RAD/WMAX 
OUT(K) -WFUNC 
1700  CONTINUE 

WRITE(06, 1705)  NTYPE 
1705  FORMAT (//, ' NTYPE- ', 13 ,//) 

CO  TO  (1710, 1800, 1730, 1800, 1830, 1740, 1750), NTYPE 
1710  WRITE(06 , 1715)  WTPED,NWTPOW 
1715  FORMATC  WTPED-  " , F12 . 5 , ’ NWTPOW- *,I3 ,//) 

CO  TO  1800 

1730  WRITE(06,1735)  WKASIR 
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1735  FORMAT ( * WKASIR- ‘ , F12 .5 ,//) 

CO  TO  1800 

1740  WRITE(06 , 1745)  BESCAL, BESEDG 

1745  FORMATC  BESCAL-* , F12.5,  ‘ BESEDG- \ FI 2 .5,//) 
CO  TO  1800 

1750  WRITE(0S,1755)  CUBK 
1755  FORMATC  CUBK- * ,F12 .5,//) 

1800  DO  1820  LL-1, IRAD+1 

VRITE(06 , 1810)  LL-l,OUT(LL) 

1810  FORMAT(I3 ,F12 .5) 

1820  CONTINUE 
420  CONTINUE 
CO  TO  1850 
1830  CALL  EXIT 
STOP 
END 
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SUBROUTINE  WEIGHT 


PURPOSE 


This  program  is  used  to  compute  the  values  of 
various  weighting  functions.  The  weighting  functions 
Include  the  cosine  on  a pedestal  to  a power,  Blackman, 
Kaiser,  Bartlett  or  triangular,  Taylor,  cubic,  and 
Bayliss . 


2 . INPUT  PARAMETERS 


Name  O/R  T 

RAD  R F 

WTRAD  R F 

IAZ  R I 


Description 

Independent  variable  for  the 
weighting  function  evaluation. 

Radius  of  the  specified 
weighting  function.  (For  a 
linear  array,  WTRAD  is  the  half 
span  of  the  weighting  function) 

Flag  that  determines  whether 
the  subroutine  generates 
amplitude  or  statistical 
weighting  data. 

- ® Statistical  weighting 
specified  by  the 
probability  density 
function  defined  by  the 
chosen  weighting  function. 

• 1 Normal  amplitude  weighting 
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IWTFLC  R 


For  IWTFLC 

- 1 

WTPED 

0 

NWTPOW 

0 

For  IWTFLC 

- 3 

WKASIR 

For  IWTFLC 

- 5 

F<2®> 

0 

DB 

0 

NBAR 

0 
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Determines  the  weighting 
function  evaluated. 

- 1 Cosine  on  a pedestal  to  a 

power 

- 2 Blackman 

- 3 Kaiser 

- 4 Bartlett  or  triangular 
“ 5 Taylor 

- 6 Bessel 

- 7 Cubic 

- 8 Bayliss 


Height  of  the  pedestal  for 
cosine  on  a pedestal  to  a power 
weighting. 


I Power  of  the  cosine  for  cosine 
on  a pedestal  to  a power 
weighting . 


Kaiser  variable  for  the 
trade-off  between  main  lobe 
width  and  side  lobe  amplitude. 


F A set  of  constants  used  in 

evaluating  the  Taylor  weighting 
(See  Section  4)  . 

F The  design  side  lobe  amplitude 
in  dB  for  the  Taylor  weighting. 

I Number  of  zeros  used  to  approach 
the  ideal  pattern  function  for 
the  Taylor  weighting. 
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PII2  0 

F 

A constant,  PII2  ■ 2.0/TT2. 

BESS1  0 

F 

A constant,  BESS1  • 1.0/JO(0.0). 

For  IWTFLG  - G 

BESCAL  0 

F 

Maximum  weighting  amplitude  for 
the  Bessel  weighting  (at  RAD  - 
0.0). 

BESEDC  0 

F 

Radius  scaling  constant  for  the 
Bessel  weighting. 

For  IWTFLG  - 7 

CUBK  0 

F 

Amplitude  scaling  constant  for 
the  cubic  weighty 

For  IWTFLG  • 8 

B(20)  0 

F 

A set  of  constants  used  in 
evaluating  the  Bayliss 
weighting  (See  Section  4). 

ANC  0 

F 

Azimuth  angle  independent 
variable  for  evaluation  of  the 
Bayliss  weighting. 

DB  0 

F 

The  design  side  lobe  design 
amplitude  in  dB  for  the 

Bayliss  weighting. 

NBAR  0 

I 

Number  of  zeros  used  to 
approximate  the  ideal  pattern 
friction  fcr  the  Bayliss 
weighting . 
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iaz  - a 

XKK  0 F 


UMAX  0 F 


MAD1  0 I 


JRND  0 X 


UL  OF 


UEXT  0 F 


A thinning  factor  used  in  the 
statistical  loading.  Equals 
the  probability  of  an  element 
occurring  at  the  normalized 
peak  of  the  chosen  weighting 
function . 

Peak  of  the  chosen  weighting 
function  used  for  weight 
normalization  in  statistical 
loading . 

Starting  address  for  selecting 
random  numbers  from  the  random 
number  array  in  the  call  to 
Function  RRAND  (3) 

(1  < MAD1  < 128). 

Random  number  generator 
initialization  constant  used  in 
the  call  to  Function  RRAND 
(0  < JRND  < 236  - 1). 

Constant  used  in  Function  RRAND 
to  set  up  uniform  random  number 
distribution  (UL  - 0.0). 

Constant  used  ir  Function  RRAND 
to  set  up  uniform  random  number 
distribution 

(UEXT  - 2 .910383046E-11) . 
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3.  CALLING  SEQUENCE 

CALL  WEIGHT  ( RAD , WTR AD , VFUNC  > 

Where:  RAD  - Independent  variable  for  the  weighting 

function  evaluation. 

WTRAD  - Radius  (half  span)  of  the  weighting 
function . 

WFUNC  - Returns  the  value  of  the  weighting 
function  evaluated  at  RAD. 

4.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  Two  common  statements  are  required  for  the 
subroutine  WEIGHT.  These  contain  the  input 
variables  to  the  subroutine.  The  statements  must 
be  in  the  form  shown 


COMMON  IWTFLC , WTPED , NWTPOW , WKASIR , F ( 20 ) , B ( 20 ) , 
ANG , NBAR , BESCAL , CUBK , PII2 „ BESS1 , IAZ , XKK , 
WMAX , BESEDG 


The  labeled  common  block.  BLKRND,  is  used  in  the 
call  to  the  Function  RRAND. 

All  the  variable  names  are  the  same  as  those 
described  above  in  Section  2. 

b.  The  constants,  F(20),  are  used  in  evaluating  the 
Taylor  weighting  function.  To  decrease  execution 
time  the  constants  should  be  calculated  once  in  the 
calling  program  for  each  weighting  design  and  the 
values  stored  in  the  array  F(20).  The  equations 
for  these  constants  are  the  following: 
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NBAR-1 


m . v*  \ 

n*l  1 ~7(A2+(n+%)2)  ’ 


1 NBAR-1 

^%7  TT" 


I I ,i  . A,2  } , 

n«1  u I »ni*l»...l 

ndm  L Hn  J 


,m=NBAR+l  ..».- 


Where:  A ■ cosh"1*! 


n - 10. 0DB/a° 


«x  " HHBAR/<Aa+(NBAR-l/2)2)l/a 


The  zeros  of  the  Bessel  function 


- 0,  n - 1,2.  . . . 


c.  The  constants,  B<20),  are  used  in  evaluating  the 
Bayllss  weighting  function.  To  decrease  execution 
time  the  constants  should  be  calculated  once  in  the 
calling  program  for  each  weighting  design  and  the 
values  stored  In  the  array  B(20) . The  equations 
for  the  constants  are  given  below. 


iW|i_-  NBAR-1  r „ 

^ 17  {’  • («„> ) 


NBAR-1  . ^ 

IT  {'  • *■ 

t.  ^ 


,D)=0,1,...  NBAR-1 


,m*NBAR,NBAR+1 . . . 


Where:  - The  zeros  of  the  Bessel  function 


J^CH.n)  - 0,  m - 0,  1,  . 


Zte+JamJj tf. — 1-*»  9- -»v-  Z-iL. 
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Zn  *< 


±5, 


» n * 0 

. a * 1,  2,  3,  4 


(i  (A^n2)*5  , n - 5,  6,  ... 

* 1 - 0.9858302  + 0 .0333885 *DB  + 0 .0001'40‘DB2 
-0 .0000019 •DB3  + 0 .00000001 »DB 4 
?2  - 2.00337487  + 0.1141548‘DB  + 0 .0004159 *DB2 

- 0. 00000373 »DB3  + 0. 00000001 ‘DB 4 

*3  - 3.00836321  + 0 .00883394 *DB  + 0 .00029281 *DB 2 

- 0. 00000181  vDB3 


e*  - 


4.00518423  + 0.00501795 *DB  + 0 .0021735 »DB2 
- 0 .00003088 -DB 3 


A - cosh~lTi 

tt 


71  - 10.0DB/2° 


- wnBAR 


znBAR 


NBAR-1 


1/C  * ((P0cnr)2-l)J1  (p0c77r)- 


n«l 


£>2' 

n 


NBAR-1 


n=0 


1 - 


P„cr  9 

(^-)2 
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p0  - 0.4797212  + 0.1456692-DB  - 0 .0018739 -DB8 


+ 0. 00000218 *DB3  - 0.00000001* OB4 


The  normalization  co'nstant,  C,  is  selected  such 
that  the  weighting  function  will  produce  a peak  of 
unit  height  in  the  far-field. 


d.  The  value  of  CUBK  must  be  normalized  to  the  peak  of 
the  cubic  weighting.  The  equation  for  this  is 
shown  below. 


CUBK-CUBK/ABS ( XX  * ( XX .VTRAD) *(XX+ VTRAD) ) 


Where:  XX  - <WTRAD)8/(3 .0) 
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5 . SUBPROGRAMS  REQUIRED 

RRAND 

EXPND 

BESS 

CAM 

6 . THEORY  OF  OPERATION 

The  cosine  on  a pedestal  to  a power,  Blackman, 
Kaiser,  and  Bartlett  weighting  are  described  in 
Oppenheim  and  Schafer  <4.e>.  Details  of  the  Taylor 
weighting  function  may  be  soen  in  the  articles  by 
Taylor  and  Hansen  <4.e).  Details  of  the  Bayliss 
weighting  function  may  be  seen  in  the  article  by 
Bayliss  (4.e).  The  value  of  the  weighting  function, 
WFUNC,  for  a cosine  on  a pedestal  to  a power  is 
described  in  the  equation  below. 


WFUNC"WTPED+  < 1 -WTPED ) *<COS(RAD*PI/WTRAD*2> ) ) **NWTPOW 
For  the  Blackman  window  the  equation  is  given  below. 
WFUNC-0.42  - 0 .5  *CQS(ARG)+0 .08 *COS(ARG+ARG) 


Where: 


ARG- ( ( R AD/VTR AD ) + 1 ) *PI 


The  equation  for  the  Kaiser  window  is  given  below. 

WFUNC-EXPND ( CONK  *SQRT  < SQN-RAD  *RAD ) ) *DENOK 

Where:  DENOM  - 1 .0/EXPND(WKASIR) 

CONK  - WKASIR/VTRAD 
SQN  « WTRAD*WTRAD 

The  equation  for  the  triangular  weighting  is  given 
below . 

WFUNC-1  - RAD/WTRAD 
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The  equation  for  the  Bessel  weighting  is  given  belov. 

WFtWC«BESCAL*J0(RAD*BESEDG) 

The  equation  for  the  cubic  weighting  is  given  below. 
VFUNOCUBK *RAD * ( RAD  - WTRAD)  *<RAD+WTRAD) 
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7.  FORTRAN  LISTING 


**»«*xst«««**<*«x**mi»t*i*****m*tx* 


SUBROUTINE  WEIGHT  C RAD , WTR AD , WFUNC ) 

COMMON  IWTFLC , WTPED , NWTPOW , WKASIR , F ( 20 ) , B < 20) , ANG , 
NBAR . BESCAL , CUBK , PII2 , BESS1 , IAZ , XKK , WHAX , 
BESEDC 

COMMON/BLKRND/MAD1 , JRND , XMEAN , SIG2SQ , UL , UEXT 
DIMENSION  U ( 2® ) , BZEROl ( 20 ) 

DATA  U/l .2196699,2 .2331306,3 .2383154,4 .2410628 , 

5. 2427643, 6 .2439216,7. 2447398 ,8. 2433948, 

9. 2438927,1®. 2462933, 11. 246624, 12. 246900, 

13  .247131,14  .247334,13 .247508,16.247663, 

17  .247796,18.247920,19 .248027,20.248123/ 

DATA  BZERO1/0 .386067 , 1 .6970309 , 2 . 7171939 , 3.726137 , 

4.7312271,3.7345203,6.7368281,7.7383356, 

8.7398503,9.7408943,10.7417433,11.7424473, 

12 .7430408,13 .7433477,14 .7439856,15 .7443679, 
16.7447044,17.743003,18.7432697,19.7433093/ 

UL-0.0 

UEXT-2 . 91038304 6E-11 
IF (RAD .LE .WTR AD)  GO  TO  28® 

WFUNC-0.0 
CO  TO  200 
280  PI-3.1415926 

CON-PI/(WTRAD*2) 

GO  TO(21® , 220, 230, 24®, 250, 260, 270, 281) , IWTFLC 
21®  TEMP-1 .0-WTPED 

WFUNC-WTPED+TEMP *COS ( RAD  *CON ) * *NWTPOW 
GO  TO  20® 

220  ARC- ( RAD+WTR AD  > *2  *CON 

WFUNC-0 . 4 2-0 . 5 *COS ( ARG  > +0 . 08  *COS ( ARG+ ARC > 

GO  TO  200 

23®  DENOM- 1 . 0/EXPND (WKASIR ) 

CONK-WKASIR/WTRAD 
SQN-WTRAD  *WTRAD 

WFUNC-EXPND ( CONK  *SQRT < SQN-R AD  *RAD ) ) *DENQM 
CO  TO  20® 
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240  WFUNC-1  . 0-RAD/WTR AD 
CO  TO  200 

250  P-PI*RAD/WTRAD 
CSTRT-0 .0 
DO  256  I-1,NBAR-1 

GSTRT-GSTRT+ ( BESS < 0 . 0 , U < I > *P  > *F  < I > ) 

256  CONTINUE 

WFUNC-PII2  *<BESS1+GSTRT > 

CO  TO  200 

260  RAD-RAD *BESEDC 

WFUNC-BESS (0 . 0 , RAD ) *BESCAL 
CO  TO  200 
270  X-CUBK*RAD 

WFUNC-X  * ( R AD+WTR AD ) * ( RAD-VTR AD ) 

CO  TO  200 

281  P-PI *RAD/WTR£D 
CSTRT-0 .0 
DO  285  1-1 ,NBAR 

CSTRT-CSTRT+ < -B  < 1-1 > *BESS  < 1 . 0 , BZEROl (1-1) *P>) 
285  CONTINUE 

VFUNC-COS(ANG) *CSTRT 
200  IF(IAZ.EQ.l)  CO  TO  300 
WFUNC-XKK *VFUNC/VMAX 
WTMP-SICN <1.0, VFUNC ) 

RRND-RRAND< 1 ) 

IF < RRND .GT .ABS< VFUNC) ) WTMP-0  . 0 
WFUNC-WTMP 
300  RETURN 
END 
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MISSION 

of 

Rome  Air  Development  Center 


RADC  plans  and  conducts  research,  exploratory  and  advanced 
development  programs  in  command,  control , and  communications 
(C3)  activities,  and  in  the  C3  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  communications,  electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 
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